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Heat and Mass Transfer From a 
Rotating Disk 


The analogy between heat, mass, and momentum transfer is applied to a rotating 


disk. 


Experimentally measured mass-transfer rates from a disk rotating in an infinite 


environment under laminar and turbulent conditions are related to the corresponding 


heat-transfer process by means of an analogy method. 


The experimental analog is 


shown to eliminate difficulties associated with accurate measurements of heat-transfer 


coefficients. 


Experimental data on the effect of an adiabatic surface placed at various distances 


parallel to the disk on the transfer rate from the disk are presented. 


Observations of 


some unusual flow patterns resulting from Goertler type vortexes in the transition regime 
and from some as yet unexplained turbulent vortex phenomena are also reported. 


1 Introduction 


H. AT TRANSFER BY CONVECTION from a rotating body 
is of importance in the thermal analysis of rotating components of 
various types of machinery. The rotating disk is one of a number 
of geometrical configurations of interest because many practical 
svstems can be idealized in terms of a disk rotating in an infinite 
environment or in a housing. 

The hydrodynamic phenomena associated with a disk rotating 
in an infinite environment have been investigated theoretically 
[3], and Stuart [4]; 
, Smith [6], and 
Heat transfer by convection has been in- 
and Millsaps and Pohl- 


by von Karman [1], Cochran [2], Goldstein 
and experimentally by Theodorsen and Regier [5 
Gregory and Walker [4]. 
vestigated theoretically by Wagner [7 

1 Professor of Aeronautical Engineering, University of Colorado, 
Boulder, Colo., formerly, Associate Professor of Mechanical Engineer- 
ing, Lehigh University, Bethlehem, Pa. Mem. ASME. 

?Supervisor of Dispersion Research, National Printing Ink Re- 
search Institute, Lehigh University, Bethlehem, Pa 

> Research Assistant in Mechanical Engineering, Lehigh Univer- 
sity, Bethlehem, Pa. 

* Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of Tae AMERICAN 
Society OF MECHANICAL ENGINEERS and presented at the ASME 
AIChE Heat Transfer Conference, Chicago, Ill., August 18-21, 1958. 

Norte: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, Novem- 
ber 7, 1958. 


hausen [8] Kreith and Taylor [9] for 


turbulent flow, and experimental data for cooling in air have been 


for laminar flow and by 


obtained by Young [10] in laminar flow, and by Cobb and 
Saunders [11] in the laminar as well as in the turbulent flow 
regime 

For a disk rotating in an enclosure Pantell [12] and Jimbo [13 
have measured skin friction and velocity profiles experimentally, 
and Soo [14] and Schultz-Grunow [15] have investigated the 
boundary layer and flow characteristics theoretically; neither 
experimental nor theoretical work on heat transfer are available 
to date 

This paper contains the results of one phase of a long range re- 
search program designed to investigate transfer phenomena in 
rotating systems. In this phase of the work mass-transfer rates 
from a disk rotating at various speeds in an infinite environment 
in laminar and turbulent flow were measured and compared with 
the results of an analysis which is also applicable to heat transfer 
The effect of placing a stationary surface at various distances 
parallel to the rotating disk surface was investigated experi- 
mentally in the laminar flow regime. Some unusual observations 
of flow patterns resulting from transitional ‘‘Goertler vortexes’’ 
and from some as yet unexplained turbulent vortex phenomena 
are also reported. 


1! Experimental Procedure 


During the initial stages of developing equipment suitable for 
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thermal conductivity, 
Btu/hr ft F 

local mass-transfer coefficient at r, 
lb/hr 8q ft (b/ecu ft) 

average mass-transfer coefficient 
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lb/hr sq ft (ib/eu ft 

mass-flux parameter (m/A) 
lb/hr sq ft 

rate of mass transfer, lb/hr 

partial pressure of vapor com- 
ponent A, #/sq ft 

total pressure, #/sq ft 

rate of heat transfer by convec- 
tion, Btu/hr 


parameter (q 


radial distance, ft 

critical radius at which transition 
from laminar to turbulent flow 
begins, ft 

outside radius of disk, ft 

ideal gas constant for naphthalene 
vapor, ft #/Ib R 

V 9.7, p, shear velocity, ft/hr 

™ temperature, F 

net velocity at r relative to disk 
velocity at r, ft/hr 

radial velocity component, ft/hr 

peripheral velocity 
ft/hr 

distance from disk surface, ft 
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component, 


thermal diffusivity, 
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studying transfer processes from a rotating disk, an electrically 
heated test section was considered. The design of this test section 
involved the use of a guard heater at the rim of the rotating disk 
and sliprings for thermocouple or thermistor leads as well as for 
the electrical heating element. The initial cost of constructing 
this type of test section was estimated to be quite high and the 
reliability of measured heat-transfer coefficients appeared to be 
incertain. In order to circumvent the experimental difficulties 
involved in taking heat-transfer and surface temperature measure- 
ments, it was decided to use a mass-transfer analogy, similar to 
one employed by Sogin [16]. This decision shifted the experi- 
mental difficulties to finding a suitable material which could 
withstand rotative speeds up to 10,000 rpm, but eliminated the 
need for taking heat flux and temperature measurements in a 
rotating body, and also obviated the necessity of correcting for 
heat leakages. 

The arrangement of the equipment used in the experimental 
part of this study is shown in Fig. 1. An 8-in-diameter aluminum 
disk was mounted horizontally on a Serval “Superspeed” centri- 
fuge which could be adjusted to and maintained at speeds ranging 
from 400 to 10,000 rpm 


strobos« ope 


The speed was measured by means of a 
The temperature of the air above the disk was de- 


PREC/S/ON 


termined with a precision thermometer. The tests were carried 
out in a large closed room to eliminate stray air currents. 

Before each test run crystal grade naphthalene (C;oHs, molecu- 
lar weight 128.2, melting point 79-81 C, residue after ignition less 
than 0.001 per cent) was melted in a covered flask in an electrically 
heated oven at about 120 F. The molten naphthalene was then 
cast into the undercut and recessed center of the aluminum disk 
on whose rim a split retaining ring was mounted. The disk with 
the molten naphthalene was then cooled in air at room tempera- 
ture. After the naphthalene had solidified the retainer ring was 
removed and the surface of the naphthalene was machined on a 
lathe to obtain a smooth surface raised about '/3 in. above the rim 
of the disk 

Prior to each test the surface of the naphthalene was carefully 
cleaned of loose particles and the aluminum disk with the cast 
naphthalene coating was weighed on a precision balance to within 
0.002 gm. The disk was then mounted on the centrifuge and 
brought up to speed. The time was recorded at initial weighing, 
at centrifuge start, at attainment of the desired speed, at centri- 
fuge power off, at centrifuge stop, and at reweighing. 

Since the measured weight difference between initial and final 
weighings included losses between weighing and mounting, be- 
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and stationary shroud surface, 
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tween dismounting and weighing, and during the start-up and 
slow-down periods, appropriate corrections had to be applied to 
obtain the difference in weight due to mass transfer at the test 
speed. The losses before mounting and after dismounting were 
due to natural convection transfer. Mass-transfer rates during 
natural convection were therefore measured independently, mul- 
tiplied by the total time elapsed between the initial weighing and 
mounting and between dismounting and the final weighing, and 
the resulting quantity subtracted from the measured weight dif- 
ference. To correct similarly for losses during speeding up and 
slowing down of the disk, it was assumed that the mass-transfer 
rate varied linearly with time and that the rates of mass trans- 
fer were independent of the acceleration. These assumptions were 
checked by repeating test runs at a given rotational speed, but 
varying their total duration. Since the results were reproducible 
within the estimated accuracy of the experimental data, the 
method of correction was considered satisfactory. All test runs 
were of sufficient duration that the total of the mass loss correc- 
tions necessary for any of them was less than 1.5 per cent of the 
measured weight difference. 

The average mass-transfer coefficient for the system k, is de- 
fined by the equation 


k. = mR,T,/p,.,A (1)5 


where m is the total rate of mass transfer in lb/hi 


R, is the ideal gas constant of the vapor in ft #/lb R 
T, is the surface temperature of the solid naphthalene in R 
P,, is the vapor pressure of the naphthalene at 7’, in #/sq ft 
A is the surface area of the naphthalene in sq ft. and 


k. is the average mass-transfer coefficient in lb/hr sq ft (lb/ft?) 
In Equation (1) the vapor pressure P,,, is the saturation pressure 
corresponding to the temperature 7’, which was not actually 
measured. However, as shown in Reference [17], the tempera- 
ture measured by the thermometer differs from the disk surface 
temperature only by two small corrective terms. One accounts 
for the surface temperature depression by cooling due to sublima- 
tion. Its numerical value is about 0.1 F. The other term corrects 
for the difference between the recovery temperature indicated by 
the thermometer and that of the cylinder surface. For the rotat- 
ing disk this term is [(wr)?/2Jg,c,](1 — 7,) where 7, is the recovery 
factor. Assuming that the recovery factor for a rotating disk is 
of the same order of magnitude as for a flat plate at the same 
velocity, this correction term never exceeds 1 deg F. The omis- 
sion of these correction terms does therefore not introduce an ap- 
preciable error in the results. An analysis of the accuracy of the 
experimental results in terms of uncertainty intervals based:on 20 
to 1 odds [18], indicates that the value of any of the mass-transfer 
coefficients is accurate to within +6 per cent without applying a 
surface temperature correction. 

The physical properties used in the reduction of the experimen- 
tal data, including the diffusion coefficient of naphthalene vapor in 
air, the specific gravity of solid naphthalene, and the vapor pres- 
sure of solid naphthalene, were taken from Table 2 of Reference 
[19] which is based on the best available data 

In the tests designed to study the effect of a shroud on the 
transfer process, the experimental procedure described here was 
followed, but after the disk was mounted a 2-ft-square glass plate 
was placed parallel to and at a predetermined distance above the 
disk surface on an adjustable stand. The same data reduction 
method as for the free disk was used. 

5 The concentration of naphthalene vapor in the air above the plate 
is so small that no correction need to be made to distinguish the 
process from equimolar counterdiffusion. 
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li Evaluation of Average Mass and Heat-Transfer Co- 
efficients for a Disk Rotating in an Infinite Environment 


The experimental technique used in this investigation as well as 
that in the heat-transfer measurements of Cobb and Saunders 
{11] yielded average coefficients of mass and heat transfer. To 
compare the measured values with those calculated analytically 
it is therefore necessary to evaluate average values of the co- 
efficients. 

As shown in Appendix I as well as in References [7 and 8], the 
theory predicts that in the laminar flow regime mass and heat- 
transfer coefficients are uniform over the entire surface of a rotat- 
ing disk. For the system used in the experiments, i.e., for a 
Schmidt number of 2.4, the approximate analysis presented in 
Appendix I yields the equation 


Sh = k,r/D, = 0.67Re®-5 (2) 


Equation (2) applies also to heat transfer for a fluid having a 
Prandtl number of 2.4 if the Sherwood number is replaced by the 
Nusselt number h,r/k. 

An exact solution of the complete Navier-Stokes and energy 
equations for the rotating disk was obtained by Millsaps and 
Pohlhausen [8]. For cases where viscous dissipation can be 
neglected, the usual situation in laminar flow, their results can 
be reduced to the equation 


Nu = Sh = CRe?®-5 (3) 


where the constant C depends on the Schmidt or Prandtl number 
as shown in Table 1 
Table 1 


Pr or Se 0.74 1.0 2.5 
Cin Eq. (3) 0.33 0.39 0.60 


5.0 
0.80 


As shown in Reference [4], laminar flow exists at Reynolds 
numbers wr?/y below about 2 K 10%. At about this value of 
Reynolds number disturbances in the laminar boundary layer on 
the disk are amplified and rapidly cause transition to turbulent 
flow. At a Reynolds number wr,?/v larger than about 2.8 * 105 
the flow is completely turbulent. On a rotating disk having a 
Reynolds number larger than 2.8 10° there exists therefore a 
laminar inner core surrounded by a ring shaped transition region 
and a ring shaped turbulent regime. 

In the turbulent regime mass and heat-transfer coefficients in- 
crease with increasing radius. As shown in Appendix II, a 
modified mass, heat, and momentum-transfer analogy predicts 
local values of Sherwood or Nusselt numbers given by the equa- 
tions 

/ 
te Se B,, V Cp,/2 


Sh = 


58,Se + 5 In (58, Se + 1) + W2/Cp, — 14 


Re Pr By, WC >p,/2 


+ 5 In(58,Pr + 1) 4+ 


Nu = 7 


/ 


56 ,P1 Vv 2, Up, 14 


respectively, where the symbols are defined in the nomenclature. 


To evaluate average values of the mass coefficients the local 
values must be integrated over the entire disk surface according to 


the equation 
To 
k..2mr dr + kar? 


e 


aT," 


° 
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Table 2 Summary of results for mass transfer from disk rotating in an infinite envir 





Temperature Rate of mass Average 
of air, transfer, Sherwood Reynolds number 
Test run no Speed, rpm deg F gm /min number based on disk radius 
2900 79.7 0 02546 293 1.98 X 10° 
4500 0.04074 545 3.10 
5500 0.05776 773 3.78 
2000 0.01762 231 37 
7000 0.08817 1154 81 
7200 76 0.08520 1216 96 
1000 5.{ 0.01137 162 691 
2000 72.5 0.01421 240 40 
3000 2.4 0.01812 305 09 
6520 ¢ 0.06238 937 52 
500 9.8 0.00568 114 352 
1200 . 0.0101 179 839 
3600 79.6 0.03364 388 45 
3300 0.02931 373 26 
8300 7 0.10949 1397 
t 10300 , 0.15111 2005 
17 700 0.00852 144 
18 1500 0.01065 202 
19 1730 0.01169 220 
20 2250 0.01896 258 
21 2600 0.02144 275 


Joisjci=i 
CAIUINNCOCOKRN KEK Of ee 


——s 7) 
cr hy 


* In this test a small piece of solid naphthalene chipped off near the rim and it was estimated that 
the actual mass-transfer rate is about 6 per cent less than the measured value 


where r, is the radius corresponding to the critical Reynolds The average Sherwood number corresponding to the average 
number value of the mass-transfer coefficient was then evaluated accord- 
The same procedure must be followed to evaluate the average ing to 
and h., 


derived in Appendix II cannot be integrated analytically, the 


heat-transfer coefficient Since the expressions tor Bes Sh = E.r,/D, (g 
averaging Was done numerically according to A similar procedure was followed in the evaluation of the average 
Nusselt numbers. 


Ar 
IV Presentation and Discussion of Results 


The experimental results for mass transfer from the disk rotat- 
ing in an infinite environment are summarized in Table 2. The 
speed was varied from about 500 to 10,000 rpm, corresponding to 
a rotational Reynolds number variation from about 3.5 K 10‘ to 

In the averaging procedure the transition region was neglected 


- 7 < 105. At the highest Reynolds number about 70 per cent of 
and it was assumed that turbulent flow sets in at the radial dis- ‘ . sinh 


, the total disk surface was covered by a turbulent boundary layer. 
tance corresponding to the critical Reynolds number (see Fig. 2 - hs one 

rhe test results are presented graphically in Fig. 3 where the 
average Sherwood number is plotted as a function of the rota- 
tional Reynolds number. Also shown in Fig. 3 are the results of 
theoretical analyses for the laminar flow regime and the results 
of the analogy calculations according to Equations (7) and (8). 
In addition to the mass-transfer results, averaged heat-transfer 
data obtained by Cobb and Saunders [11] with a disk rotating 
in air are compared with the results of the analogy calculations in 
Fig. 3. 

In the laminar regime the experimental] data fall slightly below 
the curve predicted by the analysis in Appendix I, but slightly) 
above Millsaps and Pohlhausen’s theory [8]. When the Reynolds 
number exceeds the critical value and part of the disk surface is 
covered by a laminar boundary and part by a turbulent boundary 
layer, the experimental mass-transfer data fall slightly below the 
results predicted by the analogy if a critical Reynolds number of 
2 X 10° is used, but above the calculated curve if 2.5 * 105 is 

taken to be the critical Reynolds number. In view of the numer- 
ous assumptions made in the analysis the agreement between 
theory and experiment is quite satisfactory for a Schmidt number 
of 2. 4. 
ceaainienibin ; The experimental heat-transfer data, that is the average Nusselt 
) TURBLLENT numbers for the disk, agree in the laminar regime within 5 per 
cent with the results predicted theoretically in References [7 and 
} 8]. At turbulent Reynolds numbers the heat-transfer data fall 
Fig. 2 (A) Velocity profile with observer riding on disk; (8) boundary —** much as 10 per cent below the result predicted on the basis of 
layer profile the analogy calculations for a Prandtl] number of 0.74 and a 6, of 
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Fig. 3 Average Sherwood number and average Nusselt number as a 
function of rotational Reynolds number for a disk rotating in an infinit 
environment 





unity. One possible reason for this deviation is that the analogy 
calculations take no account of variations in physical properties 
which are insignificant in mass transfer, but should be taken into 
account in heat transfer. Moreover, since the disk with which 
the heat-transfer data were obtained was electrically heated, the 
surface temperature may not have been uniform. It is of course 
also possible that the assumptions made in the analogy calcula- 
tions could be responsible for the deviation between theory and 
experiment. 

Some observations incidental to the primary objective of the 
research are worth mentioning. At the end of several tests in the 
turbulent Reynolds number range the imprints of certain features 
These 
and Walker 
[4] by means of the china-clay evaporation tec hnique in the transi- 
tion 


of the flow could be observed on the naphthalene surface. 


traces corroborate the observations made by Gre gory 


regime According to these authors a vortex system is 


generated in the region of flow upstream of transition, and between 
two critical radii vortexes occur which are 


stationary relative to 


the moving surfac Also the imprints left on the naphthalene 
surface indicate the presence ol fairly level scale vortex pattern 
locked to the surface in the transition regime. However, whereas 
in the British experiments the vortex 


were visible in the turbulent regime, the imprints on the naph- 


{| no traces oi pattern 
thalene suggest that also in the fully developed turbulent regime 
vortexes remam attached to the disk surtace and pe netrate the so- 
The 
ippears to be much 


called laminar sublaver in a more or less regular manner. 
scale of the vortexes in the turbulent regime 
the existence of 


smaller than in the transition regime. To verify 


stationary vortexes in the turbulent zone the tests should be re- 
peated with a disk whose surface is covered with china-clay, but 
instead of looking for traces merely after completion of a test, 
the streak pattern should be recorded photographically through- 
out the test. It is planned to perform such tests in the future. 
The effect of a shroud on the mass-transfer process from a 
rotating disk was studied experimentally in the laminar flow re- 
gime at rotational speeds of 1000 and 2000 rpm, corresponding 
to rotational Reynolds numbers of about 7 104 and 1.4 X 105, 


The test results are summarized in Table 3 and presented graphi- 
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cally in Fig. 4 where the ratio of the average Sherwood number 
for the shrouded disk to the Sherwood number for the free disk 
at the same rotational Reynolds number Shy/Sh,, is plotted 
versus the ratio of the distance between the disk and the shroud 
surface and the disk diameter, Y /d. 

An inspection of the results shows that within the range 
of Reynolds number covered in the experiments the effect of the 
Reynolds number on the Sherwood number ratio is less than the 
scatter of the data. The presence of a shrouding surface de- 
creases the rate of mass transfer. Ata clearance to diameter ratio 
of 0.25 the mass transfer for the shrouded disk is only about 10 
per cent less than that for a comparable free disk, but drops to 
80 per cent of the free disk mass-transfer rate at a clearance ratio 
of 0.05. As the clearance is further decreased the mass-transfer 
rate decreases rapidly. 

These observations are in qualitative agreement with analyses 
of the boundary characteristics in this sort of system [14, 15]. 
When the fixed boundary is sufficiently far away from the surface 
of the rotating disk there will be two distinct boundary lavers, one 
on the stationary, the other one on the rotating surface, leaving 
between them a “core’’ of finite width which rotates at approxi- 
mately one half the velocity of the disk. In this regime the 
presence of the stationary surface would not be expected to affect 
the heat or mass transfer appreciably. However, as the clearance 
is decreased a point is reached where the two boundary layers 
meet. Any further decrease in clearance will impede the circula- 
tion appreciably and eventually cause choking of the central core. 
In this regime the rate of heat or mass transfer would be expected 


Qo é 


2 o 
Wa 


Fig.4 Ratio of Sherwood number for shrouded disk to Sherwood number 
for free disk as a function of clearance ratio 


Table 3 Summary of results for mass transfer from a disk rotating in a 
shroud 


Dimensionless 
shrouding ? 
Test run listance __Tatio 
no rpm (Y /d) (Shy/Sh. ) 
a L000 0.143 0.871 
b 1000 0.059 0.823 
c 1000 0.129 0.878 
d 1000 0.031 0.717 
1000 0.094 0.872 
1000 250 0.896 


Sherwood no. 


Speed, 


2000 079 0.898 
2000 031 0.754 
2000 035 0.776 
2000 023 0.741 
2000 O16 0.651 
2000 012 0.585 
2000 250 0.926 
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to drop off rapidly with decreasing clearance. The experimental 
results corroborate these phenomena qualitatively, but additional 
studies of the flow and boundary layer characteristics are neces- 
sary before the various flow regimes can be delineated quantita- 
tively 


V Conclusions 


1 The analogy between heat, mass, and momentum transfer 
can be applied to a disk rotating in an infinite environment to cal- 
culate heat and mass-transfer coefficient from friction coefficients 
in the range of Prandt!] and Schmidt numbers between 0.7 and 


orc 


a. 

2 In the laminar regime of rotational Reynolds numbers, a 
stationary surface placed parallel to the surface of a rotating disk 
reduces the rate of mass transfer, but the reduction is less than 
20 per cent of the transfer rate for a free disk at clearance to 
diameter ratios larger than 0.05. 

3 Traces left by the air flowing over a disk rotating in an in- 
finite environment suggest that stationary vortexes occur in the 


turbulent as well as in the transition regime 
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APPENDIX I 


Approximate Calculations of Rates of Heat and Mass 
Transfer From a Disk Rotating in an Infinite Environment 
in Laminar Flow 


In cylindrical co-ordinates the steady-state equation for mass 
transfer of component A by molecular diffusion through an inert 
gas in laminar motion is 


0%c 1 oc Oc 0c 
dD = + - — Uu, 
or* r or oy" or 
where D, is the diffusivity of component A in ft?/hr and c is the 
concentration of component A in the inert gas in lb/cu ft. 


The corresponding equation for heat transfer to or from a rotat- 
Ing plate is 


(= 1 or ofr oT oT 
a + — + - uU, - - - = 
or? r Or oy? or ~ OY 
if viscous dissipation is neglected and the fluid properties are 
uniform. Assuming that the temperature or the concentration 
at the surface of the rotating plate is uniform and that the tem- 


perature or the concentration in the fluid far away from the plate 
is zero, the boundary conditions pertaining to the system are 


0 10 


at y = 0: orc = ¢, 


aty = ©: orc 0 

Von Karman [1{ and Millsaps and Pohlhausen [8] have shown 
that the boundary layer thickness 6 is constant over the entire 
plate, provided the small influence of flow toward the plate is 
neglected; and, that the radial flow velocity u, at a given distance 
y from the disk surface is proportional to r. Consequently, the 
total mass rate of fluid flowing outward in the radial direction at 
any distance r is proportional to the area mr*?. Thus a concentra- 
tion or temperature distribution dependent on y, but independent 
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of r may be assumed [7]. Equations (9) and (10) reduce then to 
the ordinary differential equations 


d*c 
($5): 
dy? 
aT 
(a) dy? 


Equations (11) and (12) yield upon integration, 


de 1 y ie 
- a -—6x ual c 
dy P D, 0 ve 
dT 1 v j 
= —ex u,dy +c 
dy Pla Jo _— 


Another integration gives 


c= f exp D f u,dy + c) ay 
y ° 0 
oo 1 y 

T -f exp f u,dy + C )-dy 
y a 0 


The integration constant C can be evaluated from the boundary 
condition at y = 0 while the boundary condition at y = @© is 
automatically satisfied by the upper limit of the first integral in 
Equations (15) and (16). 


(11) 


(12) 
respectively. 


(13) 


(16) 


The concentration or the temperature 
at the plate surface are according to Equations (15) and (16), 
respectively 


(17) 


l u 
T,= f exp ( f u,dy + ¢ ) ay 
0 a 0 j 


On division of Equation (13) by Equation (17) or Equation (14) 
by Equation (18) the integration constant C drops out and 


1 (¥ 
exp ( vs) 
; D, Jo _ 
a cD l (ey , F 
ex u,ay al 


SS «& 
exp 7 | udu) 
i 0 


—T.—— 


‘ wo 1 ay 
f exp ( : ws) dy 
0 a Jo 


For naphthalene vapor diffusing from the solid naphthalene-air 
interface at the disk surface into the surrounding air the concen- 


—¢ 


tration in the gas phase is related to the partial pressure p, by the 
equation 
c = p,/R,T (21) 


where p, is the partial pressure of the naphthalene vapor in the 
air in psf, R, is the ideal gas constant in ft #/lb R, and T is the 
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absolute temperature in deg R. The mass-transfer coefficient k, 
and the heat-transfer coefficient h, are given by 


1 
ke 3 | tins, ( =) ] : 7 [- 7 ( 2 ) | 
é. Oy / yo Py Oy / y-0 


/ 


/ re l y 
= D, f exp ( - f u ay) ay (22) 
i Je WNiTe 
- _ 1 y 
[ _, (2 ) | =k exp ( - f ud )-d 
Oy | yo 0 a Jo 


(23) 


and 


The corresponding dimensionless numbers, the Sherwood 
number and the Nusselt number, are 


k. - 1 v 

c r / exp ( - f vst) (24) 
if 0 D, Jo 

; {re 1 y 

Nu / fe exp f u,dy )-dy (25) 
k 0 @ Jo 


The integrals of Equations (24) and (25) can be evaluated 
graphically with the aid of von Karman’s solution [1] of the hy- 
drodynamica] problem according to which 


uv 2 
fory $6, wu, = 20 [1026 (F) (1 - §) (1424) 
| dy = —2Qwy [os ( : ) 


Sh = 


‘a \3 2 
- 052( 4) - 017 (1) +051(4) 


and for 


y = —0.708(vw)'/? 


l/s 
- v Ss 
= 2.58 ( 
Ww 


The mean temperature during the experiments was 74 F. At 
atmospheric pressure and 74 F the kinematic viscosity of air v is 
0.60 ft?/hr [19] and the vapor diffusivity D, is 0.25 ft*/hr [19], 
so that Se = 2.4. With this value for D, Equation (24) becomes 
upon evaluation of the integral in the denominator 


wr \'/2 vy 
Sh = 0.67 . = 0.67 Re’? 
v 


Similarly, for Pr = 


where 


. (29) 
2.4 the Nusselt number is 


/ . (wr? / : 
Nu = 0.67 ( ) = 0.67 Re’/? 
v 


[7] evaluated the heat-transfer coefficient in air (Pr 
and found 


2\°/3 
Nu = 034 (@") = 0.34 Re'/* 
Vv 


It should be noted that both heat and mass-transfer coefficients 
are uniform over the entire surface in laminar flow. 


(30) 


Wagner 
= ().74 


(31) 
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APPENDIX II 


Approximate Calculation of Rates of Heat and Mass 
Transfer From a Disk Rotating in an Infinite Environment 
in Turbulent Flow 


The analogy between heat, mass, and momentum transfer has 


been used successfully to predict heat and mass-transfer co- 


efficients friction coefficients 


Typical analogy solutions for flow inside of tubes or ducts and for 


from experimentally measured 
flow over plane surfaces are described in References [20, 21, 22 
ad ; 
23, and 24 


analogy 


Kays and Bjorklund [25] obtained recently a1 
infinite 
environment and found that the solution agreed satisfactorily 


solution for a circular cylinder rotating in an 
with experimental data for air 
The analogy method is based on the assumption that the trans- 


fer mechanisms for heat, mass, and momentum are physically 


similar, or can at least be described by similar equations of the 


d\ 
g.7,/p = (€y + Vv) 
dy 


momentum transfer and resulting shearing stress 


( q : ’ , ( dT ) 
f “ = | ~ ¢ 
—" dy 


for the rate of heat transfer, and by 


( m (< 
) = (e, + D, (34 
A}, dy 


for the rate of mass transfer, where €y, €,, and €,, are dependent 


form 


for the 


on the amount and kind of turbulent mixing at a point. To com- 
plete the analogy one assumes that €, = SByéy and €,, = B,,€y 
where 8, and 8, must be determined theoretically or experi- 
mentally inasmuch as the equation of motions for a fluid contains 
In the 


the surface the shear stress, the heat flux, and the mass 


terms which are not present in the diffusion equation. 
vieinit 
flux are nearly constant and if the fluid properties are uniform, 
33), and (34 


conditions at the surface as 


Mquations (32 can be written with reference to 


vy = (€y + 


v > § 
M (8. + - 


se L 


The application of the analogy method to heat and mass trans- 
fer from a rotating disk is greatly simplified if one assumes that 
the flow in the immediate vicinity of the disk surface is similar 
to the flow near the wall of a pipe. This assumption can be partly 
justified on the basis of available data on the drag coefficients in 
the turbulent regime [5] which are well correlated by the equa- 
tion 


l 
V Op, 


/ . 
= constant + 4.07 logw (Re V Cp,) 


Equation (38) is of the same form as the well known von Karman- 
Nikuradse resistance equation for pipe flow and the coefficient 
4.07 corresponds to a magnitude of 0.4 for the universal nondi- 
mensional profile constant for turbulent flow near a surface. 
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az 
r) 9040 so 720 200 
gy 
Fig. 5 Evaluation of buffer layer thickness by comparison of (T,+ — 
T..)/(T, — T..) versus Y * calculated from analogy theory with experimen- 
tal temperature profile measured by Cobb and Saunders, Reference [11] 


Moreover, measurements by Gregory and Walker [4] and 
Theodorsen and Regier [5] indicate that the velocity profile near 
a disk is logarithmic outside the so-called laminar sublaver and 
that the radial velocity component is of a smaller order of magni- 
tude than the tangential velocity component. 

Although boundary layers on rotating disks are three dimen- 
sional and the differential equations for energy and momentum as 
well as the boundary conditions are dissimilar, for the purpose of 
the following approximate calculations a conventional analogy 
The flow field was divided into a so- 
called laminar sublayer extending from y* = 


approximation was made. 
0 to y = 5, a 
buffer layer from y+ = 5 to y* = 30, and a turbulent region 
= 30. The upper limit of the buffer layer thickness 
was estimated by comparing the temperature distribution pre- 


bevond v ’ 


dicted on the basis of strict analogy calculations [25] with the 
temperature profile measured by Cobb and Saunders [11]. The 
solid line in Fig. 5 represents the calculated temperature profile 
while the dashed line was faired through experimental data taken 
from Reference [11]. The experimental data tollow the pre- 
tlicted curve quite well up to a value of ¥* of about 30, but beyond 
this point the experimental temperature gradient becomes less 
than the calculated value. On the basis of similar observations in 
experiments on heat transfer from a rotating cylinder Kays and 
Bjorklund [25] suggested that in the region beyond y* = 20 the 
transfer mechanism could be approximated by the simple Reynolds 
analogy Although the two systems are 
respects, the 


dissimilar in many 
teynolds analogy was found to be a satisfactory 
approximation for heat and mass transfer from a rotating disk 
beyond y* = 30. Although the following derivation is for mass 
transfer, it applies equally well to heat transfer since Equations 
(33) and (34) are analogous. 

In the laminar sublayer V*+ = y*, €, = O and integration of 
Equation (37) gives therefore 

Ac, = 5M Sc/S (39) 

In the buffer laver between y* = 

Karman-Nikuradse profile [21] is 


5 and y+ = 30 the von 
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Vt = —3.05 (40) 


+ 5.0 In y* 


which gives a dimensionless velocity gradient dV * 
substituting 5/y* 


/dy* = 5/y* 


for the velocity gradient in Equation (35) 
and solving for €,, yields 


so that 


7 ’ 
€,, = Bay ( 
0 


.) 


Equation (37) can then be written as 


dc, = . 
1/D,5e 


Integration between y* = 5 and 30 yields the concentration po- 
tential drop in the buffer laver Ac,, or 


5M 
Ac, = In 58 s 
6,8 


In the flow regime outside of y* = 30 it will be assumed that the 
intensity of transfer processes by eddy diffusion is much greater 
than by molecular diffusion and that the shear stress and mass 
flux are the same function of y so that 


Integration between y* = 30 and vields Ac, 


potential at the outer edge of the buffer | 


the concentration 


Ac 


Adding the concentration potentials 


V 20, 


noting that wr 


gives 


Ac 


> 


7) 


+ 5 In(5B,S« 


For the calculations in this paper, the 
Dhanak [26 


experimental results of 
were used to evaluate 8 ind the loeal drag co- 
«ficients were taken from Reference [5 


Similarly, the analogy gives for heat transfer 


Nu By VC 
58 ,,Pr 


+ 5 In(56,,Pr 


The selection of an appropriate value of By for the numerical 
evaluation of Equation (4) was difficult because even for a pipe or 


a duct available experimental results are inconclusive [27, 28, 29, 
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30]. Experimental measurements of 8, taken with air flowing in 
rectangular channels are summarized in Reference [27]. They 
indicate that the space average value of 8, increases with 
Reynolds number from about 1.1 at a Reynolds number of 2 X 107 
to 1.4 at a Reynolds number of 9 * 10’ but varies also with the 
distance from the surface at any given Reynolds number [27, 28]. 
Data obtained with liquid metals [29] give values of 8, in the 
range from 0.7 to 0.9. 

In the past a value of unity for 8B, [20, 22, 25] has been used 
Also 
in this study a value of 8, equal to unity was found to yield satis- 


with good success for heat-transfer analogy calculations. 


factory results, giving agreement within 10 per cent between the 
Nusselt numbers calculated from Equation (14) and the experi- 
mental results of Reference [11]. For comparison the calcula- 
tions were repeated with 8, equal to 1.4, the other extreme. 
This caused the calculated Nusselt numbers to deviate about 27 
per cent from the measured values. 

In passing it is interesting to note that semiempirical methods 
deduced originally from experimental data on the basis of a 
thought model implying a strictly two-dimensional structure of 
the turbulent boundary layer have in this work and in Reference 
[25] been applied to systems in which the boundary layer struc- 
ture must be three dimensional in nature. Since these methods 
are apparently applicable to cases in which the structure of the 
boundary layer can by no stretch of the imagination be considered 
two dimensional, it seems appropriate to question the hypothesis 
of a layer in the other 
The penetration of the so-called laminar sublayer by 
vortexes should supply additional food for thought. 


two-dimensional turbulent boundary 


cases. 


DISCUSSION 
S$. L. Soo° : 


Professor Kreith and his associates have presented an interest- 
Professor 
Kreith did not elaborate on the matter involving a stationary sur- 
face placed parallel to the disk. 


ing study with a rather unique experimental approach. 


It might be worth pointing out 
that as far as heat-transfer analogy is concerned, his system pro- 
vides an analog of a constant disk surface temperature, an in- 
sulated stationary surface, and the fluid drawn in and out from a 
constant temperature reservoir following the flow depicted in 
Fig. 9(c 


of his reference [14]. Such a system is probably not of 


practical interest from disk cooling point of view. The discusser 
would like to suggest the following two different experiments for 


the constant disk temperature application: 


1) Insulated stationary surface with radial outflow cooling 
with mass-transfer analogy provided by a nonevaporating station- 
ary surface but carry a central hole to pump air through the gap: 
Based on Fig. 10 of his reference [14], 1-in-diam hole, nonrecir- 
culating flow can be maintained for a system at '/»-in. gap, 
8-in. radius, and 1000 rpm. 

Constant stationary surface with recirculation cooling: 
The mass-transfer analogy can be provided by an intensely 
cooled stationary surface such that condensation of the evapo- 


rating medium can oecur at the stationary surface. 


Case (1 
enclosed in an insulated housing; 


in practical application is source-flow cooling of a disk 
”») 


a disk cooléd by 
water-jacketed housing. Ixperimental data on such cases will 
not only furnish useful design reference, but also provide the ex- 
perimental approximation for analytical study with infinite disk 
and infinite stationary plate. These data will help improve our 


understanding of the heat transfer of these significant systems. 


€ Associate Professor of Mechanical Engineering, 
versity, Princeton, N. J. 


Princeton Uni- 
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Myron Tribus’ 


The use of the analogy between heat transfer and mass trans- 
fer in the case of the rotating disk is a very attractive method for 
getting around the problem of measuring temperatures and pro- 
viding power through slip rings. If one is to rely upon the anal- 
ogy, however, it is necessary that the boundary conditions for 
the mass-transfer system and the heat-transfer system be identi- 
cal, The writer therefore is concerned with whether or not the 
surface of the disk was isothermal. If there occurred a variation 
in temperature over the disk surface, an even larger variation in 
the sensitivity 


vapor pressure is to be expected (because of 


of vapor pressure to temperature). These variations could very 
well change the mass transfer and therefore the corresponding 
dimensionless heat and mass-transfer moduli by 50 per cent or 
more. Does the author have any information concerning the 
temperature distribution which might have occurred over his 
disks? 


the effects of aerodynamic heating were taken into account. 


In addition, the writer would be interested to learn how 


If we can have some assurance that this use of mass transfer to 
deduce heat transfer can be applied to the rotating disk, a number 
of important problems can be solved, for example, the disk rotat- 
ing within housing, the disk with heat or mass transfer at selected 


points on the radius, and so on 


Authors’ Closure 


The authors wish to express their appreciation to Dr. Soo and 
Dr. Tribus for their interesting and constructive discussions 

The experiments in which a stationary surface was placed 
parallel to that of the rotating disk were of an exploratory nature 
und are being continued. The heat-transfer analog of the system 
investigated in the tests reported in the paper corresponds, as Dr 
Soo has pointed out, to an isothermally heated rotating disk with 
an insulated surface placed parallel to the surface of the disk 
with fluid drawn in from, and discharged into, a constant tem- 
perature reservoir. However, subsequent studies of the fluid flow 


(31 ]* have *hown that the flow pattern in the space between the 


’ Associate Professor, Engineering Department, University of 
California, Los Angeles, Calif. Mem. ASME. 
* Numbers from 31-33 in brackets designate Additional References 


ut end of this closure 


— 


| 
am 


rotating disk and the stationary surface does not correspond to 
that predicted theoretically and depicted in Fig. 9(c) of Ref. 
[14]. As reported at the Symposium on Stall at the 1958 ASME 
Annual Meeting in New York [32], when the stationary surface is 
larger than the rotating disk the flow is periodic in nature with 
distinct inflow and outflow regions rotating in the space be- 
On the other hand, 
if the stationary surface is of the same size or smaller than the 
rotating disk, the flow is steady as assumed in Ref. [14]. The 
results of the investigation of the fluid flow mechanism in the 


tween the disk and the stationary surface. 


space between a rotating disk and a stationary parallel surface 
will be presented in detail in a forthcoming paper [33] 

The results of some preliminary experiments on the influence 
of the size of the stationary surface on the rate of mass transfer 
are shown in Fig. 6. Rates of mass transfer in a system con- 
sisting of an &-in. diameter rotating disk with an 8in. diameter 
stationary disk placed parallel to the rotating surface were 
measured with various clearances at speeds of 1000 and 2000 
rpm. The results of these tests are presented graphically in Fig. 6 
where the ratio of the average Sherwood number for the shrouded 
disk to the Sherwood number for the free disk at 

teynolds number Shy/Sh, is plotted versus the clearance ratio 
Y /d. 


diameter stationary surface while the solid line, reproduced from 


the same 
The dotted line represents a faired average for the 8-in- 
Fig. 4, represents the average for the case of a very large station- 


Although the number of tests is not sufficient to al- 
low any generalized conclusions, it is apparent that the size of the 


ary surface 


stationary surface affects the rate of mass transfer appreciably. 
For a given clearance ratio and speed, the rate of mass transfer is 
larger with the 8-in. stationary surface, i.e., in steady flow, than 
with a large stationary surface where the flow is periodic 
Following Dr. Soo’s suggestions, a few test runs were also 
made with a radial outflow system. An 8-in-diameter stationary 
disk with a l-in. diameter hold was placed at various distances 
from the rotating disk and rates of mass transfer were measured 
at 1000 and 2000 rpm. The results of these tests 
Fig. 7. An inspection of the data shows that even at clearance 


are plotted in 


ratios as small as 0.02 the rate of mass transfer in the radial out- 
flow system is equal to that for a free disk. At 2000 rpm and 
clearance ratios between 0.01 and 0.02 the Sherwood number in 
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Fig. 7 Mass transfer with radial outflow—ratio of Sherwood number for shrouded disk to Sherwood 


number for free disk as a function of clearance ratio 


the radial outflow system with naturally induced flow was found 
to be actually larger than the corresponding Sherwood number for 
a free disk. These results indicate that radial outflow systems 
may indeed be of practical interest in turbines and other types of 
rotating machines where disk cooling is a problem. A test pro- 
gram designed to investigate the effect of forced flow through the 
central hole on the flow and transfer mechanisms will be under- 
taken at the I niversity of Colorado in the near future. 

Dr. Tribus’ question regarding the influence of aerodynamic 
heating on the temperature distribution over the disk surface can 
be answered by 


the following analysis Although the surface 


temperature 7’, was not measured directly, it can be calculated 
by means of a heat balance. Assuming that conduction within 
the solid and radiation to and from the surroundings are negligible 
compared to convection and sublimation, the rate of heat trans- 
fer by convection (including aerodynamic heating) to the surface 
For 
an observer riding on thedisk the heat balance becomes there- 
fore 


equals the rate of heat transfer as a result of sublimation. 


h.[T. + (nV? 29Jc, k.p,A/R,T, (45) 
where 

7 is the recovery factor 

J is 778 ft #/Btu 


d is the latent heat of sublimation in Btu/Ib 


Applying the conventional heat-mass 


P/RsirT 


transfer analog and ap- 


proximating Pasir by we get 


k, = (h./paire,)(Pr/Se /s = (h.Reirl’,/Pc,)(Pr/Se)* (46) 


Replacing kin I quation (45 by 


T, = T. + (nV? 
where K; = (Ryir/R,)(Pr/Se)"X/c, = 


thalene system. 


this expression yields 
K:p,/P (47 


1860 | 


2q Jc, 


for the air-naph- 


Journal of Heat Transfer 


Introducing now the stagnation temperature 7’, « 


T, = To, (V2, 29S, (1 n Kip,/P (48 


The temperature recorded by a thermometer placed in the air 
reservoir above the disk is T,,,. (The surface temperature of the 
disk 7, will therefore be less than 7’ 
Equation (48). 
dissipation and aerodynamic heating, the second for the cooling 
due to 


as can be seen from 


The first correction term accounts for viscous 
sublimation. The temperature depression caused by 
sublimation is less than 0.2 F and can be neglected in most 
cases. The correction term which accounts for the viscous forces 
in the boundary layer depends on the local disk velocity and may 
therefore produce a radial temperature variation over the disk 
surface 
rotating disk is, as for turbulent flow over a flat plate, approxi- 
mately equal to the third root of the Prandtl number, even at the 
highest speed used in the tests, the maximum temperature varia 
tion was less than 0.5 F. 


However, assuming that the recovery factor for a 


The effects of aerodynamic heating 
could therefore be neglected for the range of variables covered 
in this study, but they may become significant at higher speeds 
The preceding analysis could then be used to calculate the tem- 
perature distribution over the surface and to make the corrections 
necessary for the proper reduction of the data and application of 
the mass-transfer analog to a heat-transfer problem. 
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swaodo commaven text Transient Heat-Conduction Problems 
reo Lanois F Involving Melting or Freezing 


Associate Professor of 
Mechanical Engineering, New 


York Universty, New Yor Mf Part |—Method of Analysis and Sample Solutions 
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Previously published methods for the solution of one-dimensional heat-conduction 
problems with melting or freezing are briefly reviewed and weaknesses of previous 
analytic and numericai methods are outlined. Two new and more generally applicable 
numerical methods, applicable to digital and analog computation, are developed in the 
paper. Sample problems using these new techniques are solved and compared with 
the results of a widely used conventional method. An evaluation of the various methods 
and recommendations on areas of application are included 


1 Introduction their range of applicability and introduces two novel numerical 
ipproaches which extend the range of initial and boundary condi- 
tions that can be covered. They will also permit a continuous 


M... PRACTICAL HEAT-TRANSFER PROBLEMS involve a 


, dete it the tus travel : “ate ) rni 
change of phase of the material due to either melting or freezing letermination of the fusion front travel and indicate the internal 


temperature distribution with a greater degree of accuracy than 
Applications include the solidification of castings, freezing or i» : iste ts 8 


. the previously published numerical methods. The utility of the 
thawing of soil, ice formation iblation”’’ of missile skins under een published nume nethod -™ sadlovcca 


new ‘thods will be discussed in terms of both digital and analog 
verodynamic heating, and others. In general, all these problems rew method ill be discussed in terms of both digital and analog 


n ber nablv well described by unidimensional heat flux mechanizations Details of the actual implem«e ntation on digital 
can be reasonably . desc vere ) i “ns f 1e@! c 


Despite the importance of the topic, a literature search reveals ind analog equipment and questions ol stabilitv, accuracy, and 


re | atively few solutions which mav he extended to pr actical prob- convergence will be presented in au subsequent companion paper. 
lems Analytical solutions are limited due to the mathematical P . : 
complexity, while most numerical solutions have not been de- 2 Review of Methods and Solutions Presented in the 
veloped in a convenient and general form which combines good — Literature 
weuracy with computational ease 

The present paper briefly reviews the existing solutions and The applicable differential equations for one-dimensional heat 


flow with constant thermal properties in both the solid and the 
Abstracted from the thesis submitted by Mr. Murray in partial 
fulfillment of the requirements of the Eng ScD degree at New York . 
Univerait law and Jaeger {1])?, are: 


liquid region (as found in most standard references, such as Cars- 


Contributed by the Heat Transfer Division of THe AMERICAN So- 
creTy OF MECHANICAL ENGINEERS and presented at the ASMI ov, oO”, 
AIChE Heat Transfer Conference, Chicago, Ill., August 18-21, 1958 ol " Ox? 

Nore: Statements and opinions advanced in papers are to be 
inderstood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, June 4 
1U58 Paper No. 58--HT-16 Numbers in brackets designate References 


in the solid region, and: 


Nomenclature 





specific heat, Btu-lb~'-deg I fusion front (Method II), di- fusion temperature T'; 


half thickness of a finite slab, ft mensionless deg F 


thermal conductivit eat flux (positive in direction of = temperature excess at position z, 
u deg se increasing 4 stu-sec ft-? ime ¢, deg 
Btu-ft | | Bt t t t, deg | 
latent heat of fusion, Btu-lb ! number of space increments in = temperature excess 
latent to specific heat effect ratic the solid region (Method 1), di- point n, time point » 
4 mensioniess 9 = density o ie material 
l | f | t f tl t 
dimensionless 
c(7 T’, f time variable or instant, sec 
increment of space, or “lump” temperature, deg F Subscripts 
number, dimensionless ambient temperature (=7 = referring to fusion condition 
number of space increments, or Stefan problem), deg F = referring to liquid material 
lumps, in £, dimensionless - space variable or location, ft = evaluated at time mAf 
space increment undergoing . thermal diffusivity k/ep, evaluated at x = nAx 
sion (conventional method - It™sec evaluated at z = E 


mensionless position of the fusion front, ft evaluated at z = 0 


space increment containing temperature excess above the 8 referring to solid material 
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= a, (1b) 


ot Ox* 


in the liquid region, if internal natural convection can be neg- 
lected. If convection within the liquid is not negligible, an effec- 
tive liquid conductivity must be defined, which accounts for both 
convection and conduction. In soil thawing or freezing problems, 
k, and k, are bulk-average values to account for both the soil and 
the ice or water, respectively. 

At the fusion point, these equations are coupled due to the 
latent heat of fusion by: 


where de/di is the rate of travel on the fusion front, and L the la- 


tent heat of fusion, taken as positive tor a freezing process. To 
avoid changes in the physical dimensions as the fusion front pro- 
gresses, p, = p, = p will be specified 

Methods of solution previously published can be grouped into 
analytical solutions, complete difference methods (where both 
space and time variables are written in difference form), and into 
“lumped parameter” methods, wherein the region is grouped into a 
finite number of discrete space lumps, but time is retained as 
a continuous variable. In the lumped parameter method the par- 
tial differential equations are reduced to a set of coupled ordinary 
differential equations, which can then be continuously integrated 
For solutions 
on a digital computer, the complete difference method results in 


in time by active or inactive analog computers. 


the required set of algebraic equations 


The earliest analytical solutions, given by 


Neumann [1, 2], 
were limited to a semi-infinite region, initially either all solid or 
liquid, and subject to a temperature step-input at the boundary. 
Stefan [2, 3] obtained direct solutions for the semi-infinite solid 
initially at fusion temperature with a step-input temperature at 
the boundary or for the specific heat flux distribution resulting in 
a constant fusion front velocity. Other Neumann or Stefan type 

t through 12], which also in- 
clude Lightfoot and Lachmann’s analyses for the region bounded 
by parallel planes 


solutions are given in References 


A more recent extension for the case of con- 


stant fusion velocity was given by Kreith and Romie [13] for 
evlinders and spheres 

Weiner [14 treated the case of constant freezing 
front velocity in a finite slab of metal bounded by semi-infinite 


molds. 


analytically 


In practice, his generalization to consider successive 
different velocities will require numerical means to obtain a solu- 
tion. An approximate solution for the fusion front travel’ in 
semi-infinite slabs, evlinders, and spheres was obtained by London 
and Seban [15], who assumed negligible subcooling effects in the 


' : L 
solid, limiting the solutions to cases where L* = —— — 
eT, —T; 


[16] compared the London and Seban results to the 


> 1.0 


Cochran 
Stefan solution and also indicated correction factors as a function 
of L* to permit a qualitative extension of the method to other 
solidification problems. Landau [17] 
equations in a very general form and gave analytic solutions for 
the limiting cases of infinite or negligible latent heat effects. 


formulated the governing 


The severe limitations imposed on the available analytical 
solutions can be lifted by either fully numerical or by lumped 
parameter techniques. Depending on the method selected, the 
problem can be extended to include nonuniform initial tempera- 
tures, more arbitrary boundary conditions, and, in certain in- 
stances, temperature dependent material properties 

The most complete analysis of numerical methods was pre- 
sented by Eyres and his associates [18] in a pioneer paper. Eyres 
proposed three separate methods and computed examples on an 
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electromechanical differential analyzer. The most powerful 
and most accurate of the three methods is a conventional finite 
difference procedure which was also used by a number of other 
authors. Since it forms the principal basis for comparison with 
the suggested new methods, it will be reviewed briefly. 

In this analysis the total region is divided into N equally 
spaced intervals of length Az. 
expressed at its center or nodal point, is considered to be uniform 
over the lump at each instant of time. For N —1 internal nodal 
points, the temperature corresponds to a lump extending + Az/2 
from the nodal point, while the boundary temperatures at points 
n = 0, N are characteristic of a half lump of size Ar/2. 

The lumped parameter equations for this case become [from 


The temperature of each lump, as 


equation (1) ] for the nth internal lump: 


d0, A+ ~-Rh+ha 
= @ 


dt Az? 


with a, or a, used, depending on whether the point is in the sclid 
or liquid region. 
The complete difference form of this equation is obtained w hen 
On.m+1 — 8 
At 


coupled by the common boundary condition 6, = 6, = 0, where 


is substituted for d@,/dt. Both regions are 


p is the lump containing the fusion front. This boundary condi- 
tion is maintained for a period such that the heat lost by the 
lump p equals the latent heat capacity Lp. 

Using either an R-C network analog or an electronic differen- 
tial analyzer, lumped parameter solutions similar to the method 
given were presented by Martin and Bromberg [19], Kreith and 
Romie [13], and Aldrich and Paynter [20]. Dusinberre [21], 
Price and Slack [22], and Liebmann [23] used the fully numerical 
difference technique. Liebmann further developed a simple re- 
sistance analog to implement his solution. 

To evaluate the advantages and shortcomings of the conven- 
tional method, which forms the basis of most published examples, 
These 
will be discussed in detail in section 4, but it should be noted at 
this time that the method fails to describe the travel of the fusion 


the authors checked it for two critical sample problems. 


front continuously and to give reliable indications of material 
temperatures near the fusion point. This becomes of practical 
importance in the solidification of many metals where mechanical] 
and thermal problems may be significantly affected by the rate of 
solidification 

These objections are partially overcome by the remaining 
Cochran [16] simplified the 


analysis by using a single lump of variable size to represent the 


methods available in the literature. 


complete solid region (between the input boundary and the fusion 
front) for semi-infinite slabs, cylinders, and spheres all initially 
at the fusion temperature. While the technique leads to a con- 
tinuous freezing front solution, the inaccuracies due to single 
lumping and the limited results available severely restrict its 
utility 

Otis [24] adapted the concept of a moving heat source to ac- 
count for the latent heat as initially suggested by the work of 
Lightfoot and Landau. He subdivided the region into a finite 
number of space increments, which vary in size with time, but 
where the temperature in each lump is again considered to be uni- 
form at any instant. The method requires a co-ordinate trans- 
formation in terms of an artificial time variable, which limits the 
analysis to materials initially or finally at the fusion temperature. 

Forster [25] extended the conventional complete difference 
equations for fixed network spacing for the case of materials either 
initially or finally at the (3) }. 
The fusion front travel was computed continuously by using a 


fusion temperature [equation 


quadratic through the interface and at two adjoining points to ob- 
tain the slope of temperature at the fusion point. 
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Although the method yields accurate results and allows ex- 
tension to include a variety of boundary conditions, initial (or 
final) conditions are limited to the fusion point temperature, 
and computation with the particular choice of the quadratic ex- 
pansion is cumbersom« 

Douglas and Gallie (26, 27] presented a stable, unique numeri- 
cal analysis of Stefan type problems using a three-point approxi- 
mation for the time derivative, which requires an implicit itera- 
tion of the temperature distribution (see also Reference [28)). 
However, the method imposes considerable mathematical com- 
plexity and requires that Ae be evaluated in large finite steps 


rather than continuously 


3 Proposed New Methods 


The two generalized methods suggested by the authors will 
vield considerably greater accuracy than the conventional method 
especially near the fusion point, while requiring only moderate 
increases in computational complexity. One approach utilizes a 
traveling space network where the fusion point travel can be ac- 
curately determined, but local temperature-time histories within 
the medium must be interpolated. The second method extends 

it leads to a more accurate 
prediction of the transient temperatures at any internal point, 
continuously computes the fusion front travel, but is subject to a 
small and rapidly decaying starting error. 


ind generalizes the work of Forster; 


Both methods will be applicable to a region with arbitrary ini- 
tial temperatures and will permit arbitrary boundary conditions 
for specified surface temperature, convective or radiation heat ex- 
change, or for a prescribed heat flux. The only restrictions are 
that only one liquid and one solid region be present, and that the 
material thickness is finite (except in the case where 0,5 = 0, = 
0, when the Stefan problem may again be analyzed). 

(c) Method |. Variable Space Network. Presenting first the variable 
space network method, it is assumed that the solid region (x < 
€) is divided into r equally sized space increments of thickness 
Ar, = é€/,, increasing as the freezing front progresses. Similarly, 
the liquid region (z > €) remains also divided into N —r equally 
spaced intervals of thickness Az, = (EF — €)/(N 


in time 


r) which shrink 
This is illustrated in Fig. 1 for the special case of N = 
2r = 8 network intervals 

The method focuses attention on the substantial temperature- 
time derivative of each internal point, given by: 


d0 06| dz . 206 
dt |, or, dt Ot |, 
where the rate of travel of each point is related to the fusion front 
velocity in the solid by 


dz/dt 


de dt 


Zin € 


for a uniform solid network spacing 


Combination of equations (1), (4), and (5) leads to: 


d6 r 00. de 


dt Or 


6, 
a, 
dt ox* 


in the solid region for the network points n = 1, 2....(r — 1). 
The corresponding equation in the liquid region will be: 


do E—z, 00, de 


Peet] 


ay . (6b 


dt |,, E-—e ox dt oz? 


forn=r+i,r+2 N—1. At the fusion point € these 


equations are again coupled by equation (2). 


If these equations be written in a difference form, all terms on 
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Fig. 1 Description of variables in Method | 


the right-hand side should be expressed with a consistent approxi- 


mation. A consistent three-point approximation will result in: 


d6 n (On+1 _ 6, ») de (6, i- 20, — 6.41) 


s = 
e 


= on : as? 
dt € 2 dt tay 


{ 7a ) 


d6,, (N—n) (Oy41 


- 6,, 1) de 


at E — ¢) : dt 
(0,1 — 20 


+ a,(N —r)? 


» “v Bn + 


(7b 
(E — ¢€)? i 


for the solid and liquid regions in the lumped parameter form with 
the coupling equation: 
+2 — 46,4:)) 


de l 5) k 
i 2E—e) } 


dt pL | 2e 


(6, Page! 46, ) 6, 
r ky! 


+ 


The equivalent statement for the full difference form is obtained 
by substituting 


—_— a = 6 | Emti — Em, de 
for and - =——— for 


At dt |, At dt” 


These same substitutions will also apply to all subsequent equa- 
tions which are given for the lumped parameter case only. 

In addition to the differential equations, two boundary condi- 
tions must be prescribed. These may include temperature at the 
boundary or heat flux at the boundary as determined from :* 


dd _ 2a, J -—% | Qutl 


dt ™ Az; 1 Az, 836k ff 
y _ 2a, fOvi- Oy Qyitl 


dt Az, ( Az, k, J 


10, 
ran 9b) 


The boundary heat flux can be a function of time, surface tempera- 
ture or external temperature, as would be the case with convec- 
tive or radiation boundary conditions. 
When € = 0, some terms of equations (7 
infinite or Therefore, the problem must be 
started with a small assumed initial value €, and a starting tem- 
perature must be assigned to the points in the solid region. By the 
same argument, the solution must be stopped before (EF — €) = 
0 and the final time period be extrapolated. Thereafter, the 


and (8) become either 
indeterminate. 


* Further details of these derivations may be found in Reference 


[29]. 


Transactions of the ASME 





subsequent time history of the solid can be 
analysis of the conventional diffusion equation 

In the fully numerical solution the choice of the initial fusion 
front distance €) will have considerable effect on the solution 
time required; due to stability requirements a small initial At, 
corresponding to the initial values of Az,, materially increases the 
solution time for the early (and late) stages. This problem will 
not arise in the lumped parameter solution, where stability con- 
siderations do not enter and the starting value, €, is determined 
by the equipment capacity in R-C network analyzers or by the 
number of scale changes desired in electronic differential ana- 
lyzers. 

(b) Method Il. Fixed Space Network. This alternate method which 
permits the same range of boundary and initial conditions as the 
previous technique utilizes a fixed lump size as illustrated on Fig. 
2 for the specific case of an eight space network. 

At some intermediate time the fusion front € will be in the 


; . Ar Az 
qth lump (: e., at a location gAx + oz, where — i 6, < ar 


continued by an 


For all points except 0, g, and N the conventional fixed space net- 
work equations will still apply: 


dé, 6, ; = 26,, + 6, +1 3) 
= é ( 
dt : Az? 


In the conventional method, the whole fusion point lump was 
assumed to remain at the fixed fusion point temperature 0p = 0; 
here the fusion front will be traced as it travels through the qth 
lump and the temperature of the lump node will be continuously 
adjusted while the fusion front moves. 

In the lump q (that containing the fusion front) a discontinuity 
in the temperature curve occurs at z = €. Two temperatures 
are calculated for the lump gq, one by interpolation from tempera- 
tures in the solid region and the other by interpolation from tem- 
peratures in the liquid. In each case the fusion temperature and 
the location of the fusion front, € = gAzr + 4z, 
form the interpolation: 


are used to per- 


262r/Axr 


6r/Ax 
“Ges 6, 1 : 
+ 6x/Azr 


2+ 6r/Az 


) 


: 10a) 
+ 62/Az)(1 4 


62/Azx 


262/Azx 


6x/Az 


62/Azx 


- (10b) 
2 — 6x/Az 


where again 6, = 0. Although 6z may be either positive or 

negative and both equations (10a) and (106) are mathematically 

defined over the full range of 52, 6g is represented by 0,, when 0 
dz and by 8,, when 0 > dz. 

The temperatures 8,, and 6,, are substituted into the uniform 
difference equations (3) to permit calculation of 0,-. and 0,41. 
This procedure facilitates the computations as compared with the 
alternate and equivalent procedure of using 6, in divided differ- 
ence calculations. (By this latter means the numerical deriva- 
tives are obtained over a nonuniform network requiring considera- 
ble modification of the equations.) 

In order to obtain the expression for fusion front motion, the 
slopes at the fusion front are developed from the first two terms 
of the difference form of the Taylor expansion: 


o0. 0-92, ns 59 
~ . : 4. 7 4 fe) 
2Azr ; . 


7 


OL z=. 
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Fig. 2 Description of variables in Method Il 
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=~ 1 (Ar — 8x) %! 
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Use of 0,, and 8,, in equations (11) and further substitution into 
equation (2) results in the fusion front velocity: 


de 1 j 0,, — 0,2 ‘ — 26 1 + 0 
= k qs "a-¢ 4 (A. - 7 «| 
dt Lp | .[ 2Az ‘ Az? 


Cu-@ 0.5 — Ber + Oess 
ca gat ante te 


(12) 
For the starting and terminating points g = 0,1, N — 1, and N, 
three internal adjacent network points are not available and the 
equations must be reduced to alternate slope and curvature ap- 
proximations. The writing of the starting and terminating 
fusion front equations are dependent on the boundary conditions 
and both must be resolved together. 

As in Method I, conditions at the extreme boundaries can be 
defined in terms of temperature and/or heat flux. When 6o(t) is 


scribed [see Fig. 2(c)], and when q = - 8) 9 
prescribed [see Fig. 2(c)], and when g = O(€< >], 0 = 


6o(t), by definition. Only two points are available 


06 


to compute 


and the two-point formula is used: 
OF \r=« 


since 6, : , ‘ are calculated by the 


and (11 
: Ax 3Ar 
When g = 1 : ae and when @(t) is prescribed, 


methods of equations (10 


») 


6. is computed from the two-point interpolation: 


Ar Ale — Ax 


0,, = 0 + 6. — &) = 


and again is calculated from equation (13 


ze 


Similar procedures are used forg = N — 1, N, 
prescribed: 


when 9,(0) is 
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08, 


Or 


For the boundary conditions just given, a starting or termination 
error results in the calculation of de/dt due to the two-point slope 
approximation When Q(t) is prescribed, or computed 4s a part 
of the problem, three-point approximations for #, may be ob- 
tained. While the fusion front is in the interior region g = 2, 3, 
\ —2, boundary conditions are identical to those listed for 
Method I, except that both Az, andAz, are now fixed rather than 
varying network spaces 
Again the fusion front must be started at some small, finite dis- 
tance from the boundary with the starting value, €&. Selection of 
the starting value will depend on the computer scaling and will be 


discussed in the second paper 


4 Sample Solutions 


To permit an evaluation of both new methods and the conven- 
tional technique, solutions of two sample problems were worked 
out on both a digital computer and an electronic differential ana- 
lyzer. The first problem is a Stefan analysis, which is subject to 
analytical confirmation the second sample problem covers 4a 
finite slab initially at a uniform temperature above the fusion 
temperature to be solidified with large subcool ng 

The physical variables and the boundary conditions selected 
The L 
as typical for practical metal solidification problems, 
which range from a low value of 297 F for lead to a high 950 F for 


are summarized in Table | c ratio selected may be con- 


sidered 


nickel and include aluminum, copper, magnesium, steel, and tin 
Water has an L/c value of 144 F; sodium has 165 F 

Right network spaces, corresponding to seven internal and two 
boundary nodal This 


number represents a compromise between satisfactory accuracy 


points, were selected for all solutions 
and computational complexity, and it is near the upper limit 
that can be effectively programmed on an electronic differential 
analyzer, The excellent agreement found between the numerical 
solution and an analytical evaluation of the fusion front travel 
for the Stefan problem ippears to justify this selection 

The digital solutions were carried out on the Burroughs Corp., 
E-101 computer 


Electronic 


while analog simulation was performed on the 
PACE 16-31-R, 
digital computations were carried out with a time interval deter- 


Associates, Ine., augmented. All 


METHOD I 
wETHOO 0 
ANALYTIC 
CONVENTIONAL 
2°T APPROX - METHOO I 


COCHRAN - 1 LUMP 





500 
t~ SECONDS 


Fig. 3 
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Fusion front motion with time, Stefan problem 
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Table 1 Sample freezing problems 


Stefan problem 
Semi-infinite slab 


General problem 

Geometry Symmetrical finite 
slab, thickness 
2E = 1.6 (half 
slab computed) 

62.0 = O(initially at 02. = 500 F (all in 
fusion tempera- liquid region 
ture) 

&.2 = —1000 F 
step-input 
Physical constants = 0.2778 < 10 

= 500 F 
= 0.5 
=() 01389 Btu-sec ~'-ft~! F = kz 
= 0.1 Btu-lb~ deg F 
. = 500 lb-ft = Py 


Initial conditions 


Boundary conditions 6.. = —1000 F 
step-input 


> ft? sec = al 


mined from @ , which exceeds the stability criterion 


Ax? 
and provides improved accuracy. 
Fig. 3 shows a plot of fusion front position against time for 
Method I, Method II, 


Points resulting from a conventional 


digital computations using and the 
analytic solution of Stefan. 
method are shown for comparison, as are points resulting from 
Cochran's one lump parameter technique 


eli 
Irom & two-point approximation ol 


The error resulting 


similar toan approni- 
ze 
mation used by Otis) is indicated by the plot with the approxima- 
tion used in Method II. Analog solutions for Methods I and II 
showed excellent agreement with the digital calculations through- 
out and are not plotted on this or any subsequent figure except 
as noted. A comparison of errors at various times in the solution 
Table 2. In Fig. 4, temperature-time plots at the 
nodal points are presented for the digital solutions by Method II 
and the conventional method. 


18 show n in 


Temperatures at various times 
calculated by interpolation from the moving space scale resulting 
from the digital solution by Method I are shown as are selected 
points resulting from the analytic solution 

Results of the fusion front motion with time for the second 
problem are presented in Fig. 5. Here no analytic solution is 
available for comparison, and the Cochran solution does not ap- 
ply. Digital solution by Methods I and II, and points deter- 
mined by the conventional method are shown. Since there was 
some divergence between the analog and digital results with 
Method I, particularly in the later stages of computation (pre- 
sumably due to integration of noise signals and to drift in integra- 
The 
analog solution of Method II was again in excellent agreement 
with the digital] results. 


tors), the analog solution by this method is also indicated 


Fig. 6 shows temperature histories at nodal points resulting from 
digital solutions by Method II, and the conventional technique. 
Temperatures obtained by interpolation from the Method I 


SECONDS 
400 $90 





aNaurT 
METHOO I 


METHOO Ir 
CONVENTIONAL 


Fig. 4 Temperature histories, Stefan problem 
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Table 2 Errors resulting from numerical solutions of a typical Stefan 


problem 
Time, sec 13.82 55.3 221.5 498 Notes 
e, Analytic 

solution ft 0.1 0.2 0.4 


-er cent error 


0.6 


Conventional lose 1/2 
latent heat is 
more accurate 
see Ref. [22]) 


7.5 5 Time to 


Cochranllump- 6 


Method I 


Very small starting 
(digital) 


error. Started 
at « = 0.01 
Method I 


Starting error 
(analog) 


duced fast. 
Started at 
0.05 


6 = 


Method II 


Starting error due 
digital) 


two-point ap- 
prox. until «€ = 
0.15. Started at 
€ 0.01 

Method II 


Starting error due 
(analog) 


two-point ap- 
prox. until e = 
0.15. Started at 
€ 0.05 
Method ITI 

2 pt approx. ) 


These 
serve to demonstrate the continuous evaluation of the fusion 
Method I 
computation. 
Method IT offers accuracy in front travel and temperature his- 


digital solutions are also shown. sample calculations 


front motion made possible in both new methods 


vields accurate temperatures throughout the 


tory after a starting error due to the two-point approximation for 

0, - 
. This error would be reduced in the case where Q(t) is 

OF ir=me 

the boundary condition, since a three-point approximation for 

og 

would then be available. 

e 


A comparison of the number of program steps 


quired for a digital solution of the general problem to € = 0.6 ft 


OZ 


ind the time re- 


suggests that the conventional method is fastest and most con- 
venient if there is no demand for an accurate knowledge either of 
the fusion front travel or of the temperature distribution, espe- 


cially in the neighborhood of the fusion or, alternately, if 


pomt 
the latent heat effects are small. 

If comparisons were made on the basis or equ il digital com- 
puting time, the conventional method with twice the number of 
network points may well lead to the same over-all accuracy as the 


proposed new approaches, and this operation migl t be possible on 


the same size computer. No such check was performed; the re- 


maining objections to the poor representation of temperatures 
near the fusion point and the lack of continuous calculation of 
front motion may still make the new methods preferable 

R-C 
the nonlinear 
fusion front coupling equations generally suggests that the mini- 


For an inactive analog computer, such as an network 


analyzer, the mechanization required to allow for 


mum number of network points be used and the new methods 
appear preferable. The same applies to an electronic differential 
analyzer; here more than eight to ten network points are likely to 
lead to poor results due to internally induced drift errors. For a 
given number of network spaces the new methods require sig- 
nificantly more equipment than does the conventional technique. 

Further details on the choice of analog or digital computing 
equipment, on the problems inherent in a mechanization of all 


three methods on both types of computers, and questions re- 
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garding stability and accuracy will be discussed in the subsequent 
companion paper. 


5 Conclusions 


The principal results described in this paper may be summa- 
rized as follows: 


1 The inherent difficulties in formulating analytical methods 
for freezing problems with practically important boundary con- 
ditions suggest the use of analog or numerical solutions in most 
instances. Even if analytical formulations of the results can be 
achieved, the actual evaluation generally requires a major com- 
puting program, so that their practical advantage over a difference 
method is doubtful 

2 Most previously published difference methods are limited 
in their applicability and will lead to inaccurate solutions for ma- 
terials with marked latent heat effects, unless the space network 
is highly refined. 


2 


3 The new methods presented in this paper are generally 


applicable and offer the highest over-all accuracy for a given 
network spacing. In terms of computing time or analog equip- 
ment requirements, they are somewhat less advantageous than 
the conventional numerical method for the same space network, 

4 The variable space network is preferable if a continuous 
evaluation of the fusion front travel is the primary objective, 
while the fixed space network will lead to a more convenient 
Both methods are 


the conventional method in describing fusion front travel, al- 


representation ot temperatures superior to 
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though Method II introduces a starting error when temperatures 

are prescribed at the boundary 
5 The both 

equipment has been established 


feasibility of methods on digital and analog 
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A Theory of Rotating Condensation 


An analysis is made for film condensation on a rotating disk situated in a large body 
of pure saturated vapor. 
condensate outward along the disk surface, and gravity forces need not be involved. 
The problem ts formulated as an exact solution of the complete Navier-Stokes and energy 
Numerical solutions are obtained for Prandtl numbers between 0.003 and 
100 and for c,AT/h,, in the range 0.0001 to 1.0. 


The centrifugal field associated with the rotation sweeps the 


Results are given for the heat 


transfer, as well as for the condensate layer thickness, torque moment, and temperature 


and velocity profiles. 


Introduction 


Ris E the pioneering work of Nusselt in 1916 [1],! 
studies of condensation heat transfer have been largely concerned 
natural 
condensation processes, it is the gravity force alone which brings 


with the situation of gravity-induced flow. In such 
about the removal of the condensed liquid from the cooled surface 


and, in this way, controls the condensation rate. Unfortunately, 
the magnitude of the gravity force is beyond our control, and so it 
appears that there are definite limits to any natural condensation 
process. This limitation becomes especially severe when appli- 
cations at high altitudes are considered, e.g.., space vehicle power 
plants. In these situations, the gravity force diminishes to a 
negligible value and natural condensation is impossible. 

To overcome the limitations inherent in natural condensation, 
a conventional alternative such as pumping or blowing might be 


But, a 


used more intriguing idea is to create an artificial 


“gravity” by use of a centrifugal field, and this is the problem on 
The system to be studied, 


as shown schematically in the sketch, Fig. 1, 


which attention will be focused here 
is a cooled rotating 
disk at uniform temperature 7’, situated in a large quiescent body 
It will be assumed that the condensed 
Fluid in this film 
will be moved radially outward along the disk under the action 


ot pure saturated vapor. 
liquid forms a continuous film on the disk 


1 Numbers in brackets designate References at end of paper 

Contributed by the Heat Transfer Division and presented at the 
Annual Meeting, New York, N. Y., November 30—December 5, 1958, 
of Tue AMERICAN Society oF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received by ASME Headquarters, July 21, 
1958. Paper No. 58—A-70. 


field (in 
qualitatively the same way that 


of the centrifugal force 
gravity causes condensate to flow 
downward in the natural condensa- 
tion process ). 

The prime results of this investi- 
gation are the heat-transfer charac- 
Heat- 
transfer coefficients are presented 
for fluids having Prandtl num- 
bers in the range 0.003 to 100 
include the film thick- 
ness, temperature profiles, and torque moment requirements 


teristics of the system. 


Fig. 1 


Other results which are to be given 
Readers who are primarily interested in results are invited to 
pass over the section on Analysis. 

Aside from the practical aspects of the problem, there is a 
strong theoretical motivation, namely, that the problem can be 
formulated as an exact solution of the Navier-Stokes and energy 
equations. Very few cases are known which permit such rigorous 
determination of the temperature and velocity fields. 


Analysis 


General Considerations. To achieve our ultimate goal of find- 
ing heat-transfer results, it is first necessary to analyze, in de- 
tail, the velocity and temperature distributions in the condensate 
layer. The problem is, of course, governed by the basic conserva- 
tion laws: Mass, momentum, and energy; and it is these which 
constitute the starting point for our study. Since it will be sup- 
posed that the condensate is incompressible and that properties 
are constant, there will be a total of five unknowns to be con- 





Nomenclature 


specific heat of condensate 

dimensionless velocity varia- 
bles defined 
(6b) 

local heat-transfer coefficient 
q/(T sat — T,,) 


latent heat of condensation 


by equation 


thermal conductivity of con- 
densate 

torque moment 

face 


dimensionless pressure, 


fined by equation (64 
Prandtl number, c,u// 
static pressure 


local 


unit area 


heat-transfer rate per 
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radial co-ordinate 
temperature; 
temperature; 
tion temperature of vapor; 
AT = 
radial velocity component 
tangential velocity component 
axial velocity component 
co-ordinate 
tance normal 
thickness of condensate laver 
dimensionless 
variable, (w/v) /*z 
dimensionless 
laver thickness, 


dimensionless temperature, 
(Tar — T)/(T oat at 
absolute viscosity of conden- 


T,,, disk surface 

T sat, Satura- 
sate » 
Trt — T. kinematic viscosity of con- 


densate 
condensate density 


radial shear stress, 


(2%: ov, 

ab T 

or 02 
tangential shear stress, 


(ers 1 OV, ) 
- O02 r OO 


angular co-ordinate 


measuring dis- 


to disk 


sur- 


independent 


condensat« 


(w/v) (% angular velocity 
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sidered: The 


temperature 


three velocity components, the pressure, and the 
Correspondingly, there are five equations at our 
disposal: Three from the Navier-Stokes equations (momentum 
conservation ), and one aprece from mass and energy conserva- 


tion. Solving such a set of five nonlinear partial differential 
equations would appear, at first glance, a too formidable task. 
Fortunately, we can draw on the experience of von Karman [2], 
who successfully solved the simpler problem of a disk, rotating 
in « large quiescent body of single phase fluid 

The Conservation Equations and Their Transformation. The equa- 
tions expressing conservation Ol mass, momentum, and energy 
for an incompressible, constant-property fluid may be written in 


evlindrical co-ordinates as 


12] 
v 


O2* 


The viscous dissipation terms have been deleted from the energy 
equation on the assumption that their effect will be negligible 
rhe retention of these terms would still permit the formulation of 
in exact solution. The has been omitted, but it 


gravity tores 


could be included in equation (4) as part of the pressure without 
altering the analysis 
Rather than deal with this formidable array of partial differen 


tial equations, We may transiorm them to a set of ordinary dif- 


ferential equations, which are easier to solve. The new variables 
which comprise the transformation are 


a) new independent variable 


AT 


Under the transformation, equations (1) through (5) become 
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H’ = 
F” = HF’ + 


G" = HG' + 2FG 


P’ = H’ — HH’ fa) 
6” = (Pr)H@’ 5a) 


The primes denote differentiation with respect to 4 and Pr 
represents the Prandtl! number. 

Inherent in the transformation, there are two important sup- 
positions about the nature of the velocity and temperature fields 
which are worthy of discussion. First, there is the property of 
d/do = 0 
is that, except for a simple stretching of V, and V4, the 


angular symmetry, i.e., Second (and more im- 
portant 
velocity and the temperature profiles do not change shape at 
values of r. 


different In particular, the assumption that the 


temperature distribution depends only on 2 (since 7 ~ z) implies 
the existence of a condensate layer whose thickness is uniform 
over the disk 


ordinary differential equations shows that the transformation (6 


The fact that equations (la) through (5a) are 
is consistent with the basic conservation laws 

Our prime interest is in the temperature distribution and the 
heat transfer, and hence our objective is the solution of equation 
5a Unfortunately, the integration of equation (5a) cannot be 
carried out without a prior (or simultaneous) knowledge of the 
velocity function H(n). But H is intimately connected with the 
other velocity functions F and G through equations (la), (2a), 
and (3a There is no complete escape from the formidable com- 
putational undertaking, but the task may be somewhat lightened 


by combining equations (la), (2a), and (3a 


H’"’ = HH" 


G’ = 


to give 
(H’)?/2 + 2G 
HG’ H'G 


Simultaneous solution furnishes G and H, but only the latter is 


needed for integration of equation (5a) for the temperature dis- 
tribution 
Boundary Conditions. To complete the statement of the prob- 


At the 
surface of the disk (z = 0), the standard viscous flow conditions 


lem, it is necessary to specify the boundary conditions, 


are Impose d, namely, that there is no motion relative to the solid 
surface. Moreover, the liquid immediately adjacent to the disk 
1s required to take on the surface temperature - 

It has already been observed that the condensate forms a layer 
of uniform thickness, denoted by 6, over the surface of the rotating 
disk 


usual condition of negligible shear will be imposed. 


At the interface between the condensate and the vapor, the 
Further, the 
condensate and vapor will share the common te mperature T sat at 
the interface 


The formal statement of these boundary conditions is 


The thickness 6 of the condensate layer remains to be determined 
from the analysis. In terms of the new independent and depend- 


ent variables pertinent to equations (5a 7), and (8), these con- 


’ 


ditions become 


where 7, is the value of 7 corresponding to z = é. 
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With these boundary conditions, we then may proceed to 
solve equations (5a), (7), and (8), provided that numerical values 
are first chosen for two parameters: Prandtl number and 1. 
From the mathematical standpoint, this is a completely satisfac- 
tory situation. But, from the practical view, the problem is in- 
complete, since the dimensionless condensate thickness 73 would 
not be known a priori. So, there still remains the task of relating 
ns to another parameter which contains quantities which are all 
known. Such a relationship is found later. 

The Parameter ¢,AT/h,.. 
physical quantities, we 
follows 


To relate ns (and hence 6) to known 


invoke an over-all energy balance as 


6 
+ f p2trV ic. T sat 
0 


(2) 
k - mr? (10) 
OZ | mo 

The first term on the left is the energy released as latent heat, 
while the second term is the energy liberated by subcooling of the 
condensate. The right hand represents the heat transferred from 
the condensate to the disk over a span from r = Otor = r. In 
writing equation (10), the assumption of negligible heat conduc- 
tion across the vapor-liquid interface, which is standard in con- 
densation theory, has been used. In terms of the variables of the 
analysis as defined by equations (6), the over-all heat balance 
becomes 


H(n5 
Pr : (10a) 
O'(ns 
where H(ngs) and 0’(n5 
7 = 5 
From a solution of equations (5a), (7 


represent the values of H and d@/dn at 


, and (8 corresponding to 
specified values of Pr and 75, the values of H(m; 
available, and the 


and @’ ns) are 
right side of equation (10a 
c,AT/h,, is determined. Now, if the 
, and (8) are 


may be com- 
puted. So, Prandtl num- 
ber is fixed and equations (5a), (7 solved for several 
different values of ns, We are able to compute 
of values for c,AT h In other words, for a fixed Prandtl num 


ber, there is a unique relationship between ng and c,AT/h,,. So, 


i corre sponding set 


if we wish, we can think of our solutions as depending upon Pr and 
c AT /h 


Solutions. 


, rather than on Pr and 5 

. Numerical solutions of equations (5a), (7), and (8) 

subject to the boundary conditions (9b) have been carried out on 

an IBM 653 electronic computer using the numerical integration 

were ob- 

1.0, 10, and 100; 
while the last 


solutions 
tained for Prandtl] numbers of 0.003, 0.008, 0.03 


procedures outlined in reference [3 The 


the first three span the range of liquid metals, 
three correspond to ordinary liquids. Values of the parameter 


c, AT /h,, ranged from 0.0001 to 0.1 for the liquid metals and from 
i . 


0.001 to 1.0 for the ordinary liquids. Based on these solutions, 
results for the heat transfer, condensate layer thickness, torque 
moment, and temperature and velocity profiles will be given. 


Results 


Heat-Transfer Coefficients. The most important results of 
practical interest are the heat-transfer characteristics of the sys- 


tem. The !ocal heat flux fo the disk may be computed by Fourier’s 


Law: 
( ; ) 
- k 
Oz pai) 


In terms of the transformed variables of the analysis, the expres- 


sion for q becomes 
»\'/2 
MT ~ T,) ( = ) (d0/dn) nao 
v 


The derivative (d@/dy)y, obtained from solutions of equations 
(5a), (7), and (8), depends on the Prandtl number and on c,AT7'/ 
h,;,. For a particular liquid and a fixed temperature difference, 
equation (12) shows that 


(11) 


(12) 


q ~ q'/? 


It is also worth remarking that q is uniform over the surface of the 
disk , 

Introducing the definition of the local heat-transfer coefficient 
as follows 


the heat-transfer results may be rephrased in the form 


—(d0/d]) nme (14) 


The quotient (v/w)’* has the dimensions of a length, and so the 
left side of equation may be regarded as a Nusselt number. Inas- 
much as (d@/d7)y-o depends on Pr and ¢,AT/h,,, so then do the 
Nusselt number results. In presenting these results, it is con- 
venient to look first at the higher Prandtl numbers (1, 10, 100) 
and then at the lower Prandtl numbers (0.003, 0.008, 0.03). 

By evaluating equations (14) and (10a) from the numerical 
solutions, heat-transfer results for the higher Prandtl numbers 
have been plotted on Fig. 2 for values of ¢, AT’/h,, between 0.001 
and 1.0. The choice of noteworthy. 
Rather than plot the Nusselt number alone, the group 


the ordinate variable is 
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Fig. 2 Heat-transfer results for high Prandtl number fivids 
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Fig. 3 
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has bee n used 


By inspection of the figure, it is seen that for 


small values of c,A7'/h,,, all the results are represented by 


.( 
( Pr ) ‘ » 
= 0.904 15 
c,AT/hy, 


But, as demonstrated in the Appendix, equation (15) corresponds 
precisely to the heat-transfer result for the situation where 
energy convection and acceleration terms are negligible. So, de- 
partures from this relationship represent the effects of convection 
and inertia, the former tending to increase the heat transfer and 
the latter tending to diminish it. With increasing Prandtl num- 
ber, the effect of energy convection becomes relatively more im- 
portant than the inertia effects. As a consequence, the Pr = 10 
und Pr = 100 curves of Fig. 1 tend to exceed (or, in the limit, equal) 
the limiting value of 0.904. On the other hand, the inertia effects 
become relatively more important as the Prandtl number de- 
creases, and so the Pr = 1 curve tends to drop slightly below the 
limiting value of 0.904. It is worth while noting that for a sig- 
nificant part of the technically important range of c,AT/h,,, 
equation (15) is a good representation of the heat-transfer results 
of Fig. 2 

Now, we turn to the low Prandtl number heat-transfer results 
which are plotted on Fig. 3. The ordinate and abscissa variables 
are the same as on the previous figure, but here we are concerned 
with smaller values of c,A7’/h,, to correspond to the relatively 
small values of c,/h,, for liquid metals. For the smallest values of 
c AT /h 


proach (but never quite reach 


» Which are plotted on the figure, the curves tend to ap- 
the limiting value of 0.904. The 
inertia terms play an important role for these low Prandtl num 
ber liquids, tending to substantially decrease the heat transfer as 
the condensate layer thickens (increasing values of c,AT hy, 
Condensate Layer Thickness. The analysis predicts that the 
dimensionless condensate layer thickness 6(w/v)'/? will be uni- 
form over the disk and will vary with Prandtl number and ¢,AT 
h,,. This dependence may be determined from equation (10a) 
by utilizing the numerical solutions of equations (5a), (7), and (8) 
For the limiting situation of negligible inertia and heat convec- 
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Heat-transfer results for low Prandtl number fluids 


tion effects, the following limiting relationship has been derived 


(see Appendix). 
" 1/ 
* = 1.107 (a fy 
Pr 


Fig. 4 shows the condensate layer thickness results for the high 
Prandtl numbers. The limiting result of equation (16) is achieved 
for small values of c, AT’/h,,, while the deviations from this limit 
never appear too great (in common with the high Prandtl num- 
ber heat transfer). The low Prandtl number results are plotted on 
Fig. 5. There, substantial upward deviations from the limit of 
equation (16) may be observed as a consequence of the inertia 
effects. 

By comparing Figs. 4 and 5, it can be seen that for a given value 
of c, AT /h,,, the magnitude of 6(w/v)'/* increases with decreasing 
Prandtl number. Finally, we note that for a given fluid and a 
fixed temperature difference 


(16) 


6~w 
Temperature Profiles. The temperature distribution across the 
We will content 
ourselves with displaying representative results, focusing atten- 
tion on Pr = 10 for the high Prandtl number group and Pr = 
0.008 for the low Prandtl number group. 

Dimensionless temperature profiles for Pr = 10 are plotted on 
Fig. 6. The profiles for c, AT/h,, < 0.1 cannot be distinguished 
from a straight line on the scale of this figure. The deviations 
from the straight line profile increase with increasing c,AT/h,, 
(i.e., increasing layer thickness), but do not become decisive even 
for the highest value of c,AT’/h,, plotted here. The effect upon 
the heat transfer of these deviations has already shown on Fig. 2 


condensate layer may also be of some interest 


to be small. Numerous other cases were available for plotting in 
the range 0.1 < c,AT h,, S 1, but they have been omitted for 
the suke of clear reproduction. 

Fig. 7 shows temperature profiles for Pr = 0.008. For values 
of c,At/h,, < 0.02, the curves could not be distinguished from a 
straight line. Even the extreme case on Fig. 7, c, AT h,. = 0.083, 


so closely parallels the straight line profile that it was physically 
impossible to draw other intermediate curves 

It would appear from both Figs. 6 and 7 that, for most tech- 
nically interesting situations, the temperature profile across the 
condensate layer closely approximates a straight line. 
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Condensate layer thickness for high Prandtl number fluids 
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Fig. 5 Condensate layer thickness for low Prandtl number fluids 


Fig.6 Representative temperature distributions across condensate layer 
for Pr = 10 


Velocity Profiles. A view of some of the hydrodynamic aspects 
of the problem may be gained by inspection of the velocity dis- 
tributions. We confine ourselves to representative situations, 
selecting results for a layer thickness of 6(w/v)'/? = 5 to typify 
relatively thick films and those ford(w/v 
tively thin films. 

Turning first to the relatively thicker films, we present on Fig. 8 
the distribution of each of the velocity components across the 


* = (0.5 to represent rela- 


condensate film. It is seen that the tangential velocity decreases 
significantly across the layer, and that there is a substantial 
axial velocity carrying mass inward toward the disk surface. 


The general shape of the curves is not significantly different from 
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z/8 
Fig. 7 Representative temperature distributions across condensate layer 
for Pr = 0.008 


those for von Karman’s problem of a rotating disk in an infinite 
domain of single phase fluid (figure 5.10, reference [2]).? 

Velocity profiles corresponding to 6(w/v)'/? = 0.5 (relatively 
thin films Here it may be noted that the 
tangential velocity changes very little across the layer, while the 


are shown in Fig. 9 


axial velocity is appreciably smaller than that for the thicker 
layers. This reduction in axial velocity in turn leads to a decrease 
in the effect of energy convection. 

2 For increasingly thicker films, the velocity profiles should ap- 
proach more and more closely to those of von Karman's problem; 
while the temperature profiles should approach those of Millsaps and 
Pohlhausen [4], who solved the heat-transfer portion of von Karman’s 
problem. 
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Fig. 9 Velocity distributions for 5(w/v) 


rotation of the 
shaft 


i steady 


it the 


Torque Moment. [nn order to sustain 


disk, it is necessary to provide torque The torque 
moment required to overcome the shear forces due to the con- 
densate The tangential stress 


film may be computed as follows: 


on the surface of the disk T,,, is given by 


In terms of the variables of rial the expression for the 


stress becomes 
T od pr (vw dG /dn l7a 


The torque moment Vv. wsociated with condensation on one 


of a rotating disk of radius ro is 


ue = ff, 


Integrating equation (18) gives 


2M, 


: , dG /d) 9m 19) 
Tro‘p(vw*) /* 


The derivative (dG/dn)». depends upon the Prandtl number 


and ¢,AT/h,,. 
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To shorten the presentation of the results, the 
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dimensionless torque moment? is plotted on Fig. 10 as a function 
of 6(w/p Inasmuch as Figs. 6 and 7 give 6(w/v 
tion of Pr and ¢ AT h 


Irom Fig 


7 as a tune- 
, the torque moment may be computed 
10 in conjunction with one of these two figures. 

It can be easily verified that the dimensionless torque moment 
increases with increasing c, AT’/h,, at a fixed Prandtl number and 
with decreasing Prandtl number at a fixed value of « PAT hj. 

Laminar-Turbulent the present 
inalysis were to be applied, the question would immediately arise 
is to the conditions under which the laminar flow 
would be valid. 


Transition. If results of the 
assumption 
In the absence of experimental data on rotating 
condensation, we must rely on the findings for the somewhat 
similar system of a rotating disk in an infinite domain of single 
phase fluid. For that situation, it is observed (reference [2], 
p. 79) that the laminar results are usable up to a value of 


= 3X 10 


Comparison With Natural Condensation 


It may be of interest to compare the heat-transfer results for 
rotating cor densation with those for natural gravity condensa- 


* Left side of equation (19). 
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Fig. 10 Torque * for cond ti 
(radius r.)—to be used with Figs. 4 and 5 





on one side of rotating disk 


tion. Unfortunately, analytical results for the latter situation are 
available only for Pr > 1, and so the comparison must be limited 
tothat range. The local heat-transfer result for natural condensa- 
tion on a vertical surface has been analytically determined in 


references [5 and 6] as 


hz tpi, “/4 
i gc,(p — p,)x* 


Re. v6 
= | 0.68 + _ : Pr > 20) 
c, AT 
where x measures distance downward from the top ol the vertical 


} 


surtace, g is the acceleration of gravity and p, 18 vapor density. 
In general, the heat-transfer coefficients from equation (20 
would have to be compared with those read from Fig 2. How- 


for c,AT/h,, < 0.1, the comparison can be made ana- 
lytically, since equation 


ever, 
15) becomes an adequate description of 
observe that 


- . 
the results for rotating condensation. Further, we 


equation (20) approximately reduces to 
hx ty} 4 ( h , ) ‘ 
} gc, (p — p,)x' c,AT . 


Then, taking the ratio of equations (15) and (20a 


(20a 


, and attach- 
ing subscripts rot and gr to respectively indicate rotating and 


gravity !'onal condensation, we find 


h 
h 


Z.or 


rot 


21 
where we have set (1 — p, - 1. So, the local coefficient for 
the rotating disk exceeds that for the vertical surface when 


og 
Sr 
where z is the distance from the top of the vertical surface. If 
1 foot, w must exceed 3.48 radians/sec 


r= 33.2 rpm) to satisfy 


equation (22). 
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APPENDIX 
Limiting Results for Thin Condensate Layers 


For thin condensate layers, the energy convection and inertia 
tend to become negligible,‘ and the mathematical problem formu- 
lated in the Analysis section simplifies considerably. After con- 
vection and inertia terms have been deleted, equations (5a), (7), 
and (8) become 


23) 
(24) 
(25) 


where only the centrifugal acceleration represented by G? has 
been retained in equation (24). The solution of these equations 
subject to the boundary conditions (9b) is 
— 1/ns (26) 
H = »°/3 — nn 
G=1 


In particular, 


(29) 


‘ 1 
H(ns) = — a, 6'(ns) = — 
3 ns 
Introducing the expressions for H(ns) and 6’(n5) into equation 
(10a), the relationship between 6, Pr, and c,A7'/h,, may be de- 
termined as follows 
c,AT 


- Pros‘ = 
h 3 


5(*) - 3 \'% c,AT hy, v4 
v 2 Pr 
This equation coincides with equation (16) of the text. 


The local 
equations (14 


heat-transfer coefficient is easily evaluated from 


and (26), giving, 


Pr Ve 
c, AT /h,, 


of the text 


which is identical to equation (15 
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‘It may be of interest to note that these simplifications were pre- 
cisely those used by Nusselt in his original derivation for natural 
condensation on vertical surfaces. 
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DISCUSSION 
W. E. Welsh, Jr.’ 


The authors have treated a unique and original problem. It 
would be very interesting to compare experimental results for 
this situation with the analytical assumptions and results. In 
particular, the point of assuming the temperature to be inde- 
The condensate thick- 
ness and heat flux become independent of radius when tempera- 


pendent of radius seems quite arbitrary. 


ture is independent, and the results of the analysis are heavily 
biased as a consequence. Further, the von Karman assumption 
of radial and circumferential velocity components being simply 
proportional to radius is tacitly accepted as valid for this case of 
two-phase flow. These temperature and velocity assumptions 
indoubtedly simplify the analysis considerably, but it is difficult 
to emphasize the meaning of the results without some physical 
istification for this analytical model. It would enhance the 
value of this analysis considerably if the authors could cite 
physical data of any nature which would tend to support their 
issumptlions. 


Authors’ Closure 


The authors extend their thanks to Mr. Welsh for his interesting 
omments. Before proceeding to describe some recently pub- 
lished experimental data, we wish first to explore the meaning of 
the mathematical solutions which have been obtained 

To facilitate this discussion, we turn for a moment to the much- 
studied and apparently well-understood problem of boundary layer 
flow over a flat plate The flow is described by 

and by the 


both of which are partial differential equations 


a boundary 


layer momentum equation continuity equation; 
In seeking a 
solution, Blasius made the a priori assumption that the shape ol 
the velocity profile is similar at al! stations along the plate He 
found that this assumption was compatible with the governing 


equations ind boundary conditions and was able to find a solu- 
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tion. Now, the solution found by Blasius is a solution to the 
boundary layer equations for flow over a flat plate, but (to the 
knowledge of the authors) it may not be the only solution. That 
is, uniqueness has not been proved. But, since other solutions 
have not been found, we are quite content to use Blasius’ results 

Now, returning to the rotating disk, von Karman’s analysis 
for the flow in a single phase fluid is again based on the a priori 
supposition of similar velocity profiles. This assumption is com- 
patible with the conservation equations and boundary con- 
ditions; and his solution is a valid solution of the problem. No 
proof of uniqueness exists; but because other solutions have not 
been found, we are ready, without reservation, to apply von Kar- 
man’s analysis to the laminar flow prceblem. 

Now, for rotating condensation, the same point of view has 
been adopted. The a priori assumption of similarity is com- 
patible with the governing equations and the prescribed boundary 
conditions and leads to a valid solution. No other solutions 
currently exist to compete with this similarity solution. The one 
aspect of the analysis which one may question is the boundary 
conditions at the edge of the condensate layer. If the liquid- 
vapor interface is taken to be smooth (and this is a standard 
assumption in condensation theory ), then the authors see no more 
reason for challenging the similarity assumptions in the two- 
phase problem than in von Karman’s single phase study 

In the period after this report was submitted for publication, 
an experimental study of condensation on a rotating disk was re- 
ported in a PhD thesis at North Carolina State College. By 
the use of dye injection, the flow was shown to be azially sym- 
metric, Heat-transfer experi- 
ments were carried out using methanol, ethanol, and Freon-113 
at atmospheric pressure. 


as was assumed in the analysis. 


The experimentally-determined heat- 
Phis 
rhe 
fact that the theorical prediction is higher than experiment Is due, 


transfer coefficients fell about 30 per cent below theory 


level of agreement is quite common in condensation studies 


it least in part, to the neglect of the shear forces at the liquid- 
The authors 
inalvsis to include this effect 


vapor interface are currently reformulating the 


Nandapurkar, “Condensation on Rotating Surface 
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Laminar Heat Transfer in 
Rectangular Channels 


Three cases of laminar heat transfer with linear heat input in long rectangular channels 


have been treated by the method of orthogonal trigonometric series. 


The boundary con- 


ditions of the first two problems are those of laminar and slug flows with two opposite 
faces as secondary extended surfaces. A new fin parameter K defined as (wk,,/bk,) 


has been shown to be 


the important factor in governing the Nusselt number. The third 


case discussed is the combined effects of free and forced-convection in vertical rectangular 


tubes. 


Introduction 


™ PURPOSE Of this paper is to present a theoretical 

treatment of the heat-transfer problem with fully established 
The study of the 
The first phase con- 
siders two parallel plates with interconnecting longitudinal ribs 


velocity profiles in long rectangular channels 
problem is carried out in three main phases 
forming a series of rectangular channels. The fluid flowing in 
these channels is assumed to be steady and incompressible and 
to have a fully established laminar velocity profile. The physical 
properties of the fluid are assumed constant. The inner surfaces 
of the parallel plates are assumed to have a constant temperature 
in any plane transverse tothe flow. Heat is transferred to the fluid 
along the parallel-plate suriaces at a unlorm ! 
the flow direction 


ate with respect to 
For channels of sufficient lk ngth, this results 
in a linearly increasing temperature of the enclosing walls of the 
flow passages. On the surface formed by the interconnecting ribs 


or the secondary surfaces, heat is transferred from the parallel 
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Nomenclature 


plates on the primary surface. On those fins, the temperature is 
therefore not uniform in a transverse plane. In the longitudinal 
direction, however, the temperature is linearly increasing at the 
same rate as the primary surface. 

The problem considered in this paper is in reality a generaliza- 
tion of a previous solution presented in reference [1]! wherein 
limiting values of the Nusselt number are given for the case of 
constant temperature on all surfaces in any transverse plane. 
The present analysis includes the effectiveness of the transverse 
fins and includes the results in reference [1] as special cases. The 
system described here might be found, for example, in a counter- 
flow heat exchanger with equal thermal capacity for both fluids, 
resulting in linear temperature variations on both sides 

The second phase of this analysis treats the case of slug flow 
through the same physical arrangement with identical thermal 
boundary conditions. The results obtained herein are of use in 
the case of heat transfer in fully established turbulent flow of 
liquid metals, since it has been demonstrated in reference [2] that 
a correlation of the actual Nusselt number, the limiting Nusselt 
number for slug flow, and the Peclet number agrees quite well with 
experimental data. 

The third problem deals with the superimposed effects of 
natural and forced-convection heat transfer for fully established 


laminar flow in vertically mounted rectangular channels. The 


1 Numbers in brackets designate References at end of paper. 





area, sq ft 
side dimensions, ft 


fin 


specific heat at constant pres- 7 ke /bk 
: WK C 
sure. Btu per (slug-deg F) a” 


temperature gradient in flow 
direction, deg F per ft 
pressure gradient in flow direc- 
tion, lbp per cu ft 
hydraulic diameter, ft 
correction factor, see equation 
28), dimensionless 
hody force, lbp per cu ft 
dimensionless heat generation 
39 
gravitational acceleration, ft per 
hr? 


; : = Rayleigh 
function, see equation 
less 


heat-transfer coefficient, Btu per 
(deg F-sq ft-hr) 

wall thermal conductivity, Btu 

deg F-ft-hr) 


mean axial 


mn 


per ° 
function 
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fluid thermal conductivity, Btu 
per (deg F-ft-hr 
parameter, 
pressure drop parameter, equa- 
tion (39) 
dummy indexes, integers 
Nusselt number, (hD/k,) 
pressure, psf 
heat generation rate, Btu per 
(hr-cu ft) 
function as defined in equations 
(22) and (26) 


= temperature, deg F 
= axial velocity, fph 


dimensionless velocity (u/u,,) 
Fourier coefficient for velocity 


one half of fin width, ft; refers 
to wall when used as a sub- 
dimensionless, script 
co-ordinates, ft 
dimensionless co-ordinates, 
(r/D), (y/D 


axial co-ordinate, ft 


, dimensionless 


aspect ratio, (b/a 
of 
tion, per deg F 


coefficient density contrac- 

Laplace operator 

temperature-difference function, 
deg F 

number, dimension- viscosity, lby-hr/ft? 

density, slug per cu ft 

dimensionless temperature func- 
tion 

mixed mean function of @ 

Fourier coefficients for @-func- 
tion 


velocity, fph 


thermal diffusivity, (k,/pC,) 
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effects of superposition have been discussed in references [3 and 
4]. Reference [3] considers a circular tube and presents a solu- 
Ref- 


however, treats the case of rectangular tubes but 


tion for circular symmetry in terms of the Bessel functions. 
erence [4] 
neglects the influence of the heat-generation function within the 
fluid, as is necessary in a number of engineering applications of 
the heat-transfer data. The same problem of combined effects is 
also discussed in references [5 and 6]. The problem discussed 


by Ostrach in reference [5] is concerned with a two-infinite- 
parallel-plate system, considering the effects of heat generation 
and frictional heat dissipation. The work of reference [6] is ex- 
perimental and is valid for a circular tube with uniform wall 


temperature 


Influence of Secondary Surfaces on Limiting Laminar 
Nusselt Number for Constant Axial Temperature Gradient 


Consider the flow configuration shown in Fig. 1. The problem 


can be simply stated as follows: Since the interconnecting ribs 
do not have infinite thermal conductivity, passage of heat from 
the primary surface to the extension surface and subsequently to 
the fluid would necessitate a temperature drop along the fin 
length The question is What is the parameter which governs 
the nonuniformity in the temperature distribution and what is 


the influence on the heat-transfer coefficient? 
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Fig. 1 Flow channel configuration 

As a first approximation, one may estimate the heat-transfer co- 
efficient on the basis of the uniform wall-temperature data in 
reference [1]. Subsequent correction can be made by allowing 
the effectiveness of the fins as described in any standard heat- 
transfer book such as reference {7 | 


{ Such a procedure is based on 
the concept that the surface coefficient is uniformly distributed 
and is unaffected by nonuniform distribution in temperature. 
Analysis in this section shows that the values and the distribution 
of the heat-transfer coefficients on all surfaces, primary and 
secondary, are strongly influenced by the temperature distribution 
on the secondary surface. 


Anaiysis. 


The assumptions employed here are summarized 


18 follows: 


a) Fluid properties are constant with respect to temperature. 

b) The velocity and temperature profiles are fully developed? 
and the motion is laminar 

(c) The axial temperature gradient is constant and independ- 
ent of position in the channel cross section. 

? The assumption of fully developed flow for laminar motion implies 
vanishing velocity components transverse to the axis. A rigorous 
proof ean be found in reference [8] 
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With these assumptions, the basic equations of motion and of 
energy are reduced to 


V% 
V *u 


(3) 


For the primary surface the temperature is a function of the axial 
position only, i.e. 


where fp is a reference temperature at z 

For the secondary surface the temperature is a function of the 
As- 
suming that the temperature is uniform across the secondary 
surface but varies along the heat-flow path, the temperature must 
obey the following heat-balance equation: 


O*t,, ol 
wk,, ( = +k ( ) 
oy? Or 


where + sign applies toz = 


axial position z and its location from the primary surface. 


+a/2 and — sign, r = a/2. 


Equations (1) through (5) may be transformed by introducing 


the temperature-difference function 6, defined as 
O=t— (e+ Ciz 

become 

v9 

Viu = C2/p 


Hence equations (1) and (2 


Cyu a 


with the following boundary conditions 
0 


6, 0, for y = VU, b only 


0°64 o¢ 
wk,,, ) = +k , for zx 
oy? Or 


Solution. The method of solving the foregoing equations is 
similar to the procedure recommended in reference [1], i.e., by 
8) resulting in 


combining equations (7) and 


VO = C,C2/ap 10 


\ 


where (,, C2, a, and yw are independent of the cross-sectional co- 


ordinates, r and y. The corresponding boundary conditions ex- 


pressed in terms of a single function 6 are 
v°0 z= +a/2; 


A = 0, b 


(2%) ( OY] ) 
wk . = +k, > = 
eCy* or 


The solution to equation (10) can be assumed to be made 
two functions 6, and 6, satisfying the following relations 


up ol 


V'O, = C,C2/ap 
v7. = 0 
with 
and 6, =0 


on all boundaries 


=0 on y =0,b 
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00 06 

wk,, = +k, ; z (14) 
oy? “ \ ox 

The functions satisfying the foregoing relations can be obtained 


by the use of double Fourier series and the method of undeter- 
Thus, 


6, = 16C,C.b* 1s cos 
apr 


¢. mn(m*y? + n 
m,n=1,3,5 ! 


mined coefficients. 


mnrr/a) 
2)\2 


sin (nmy/b) (15) 


16C,C2b5 


amryua a nT 
nT cosh 
h 2¥ 


NX nwy 


cosh sin 
h 


and the velocity function u is given by 


16C,b? 
ur 


m,n 


Results. A commonly used heat-transfer parameter is 


Nusselt modulus, 


Nu = (hD/k (18) 


where D is the equivalent or the hydraulic diameter and h is the 
average peripheral heat-transfer coefficient defined as 


(19) 
lefined as the dif- 


fluid and 
section. 


The temperature difference in equation (19) is 
ference between the mixed-mean temperature of the 
suriace 


the primary 


Hence 


temperature at the same 


At=06,. = f Budxrdy { 
: ‘ Ja 


Combination of the preceding equations 15) to (20 


(20) 


together with 
the equivalent diameter D yields the equation for the average 


Nusselt number, 


where 


n tanh ( 
Ua 
“7 


m,n=1,3,5 
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Fig. 2 Correction factor for finite fin conductivity 


Thus equations (21) and (22) define the Nusselt number for any 
aspect ratio and for any combination of the fin geometry and 
thermal-conductivity ratio. Let 


K = (wk,,/bk,) 


j 


Hence equations (21) and (22) can be expressed as 


Nu = (Nuj-a) F.(K,¥) 


where 
FAK, y) = 1/11 + R) 


Thus, when the thermal conductivity of the fin material is in- 
finite, the peripheral temperature variation becomes zero, by 
virtue of equation (5), indicating the case of uniform wall tem- 
perature at the same cross section. The limiting Nusselt number 
in this case is given by equation (21) with R = 0; the resulting ex- 
pression is identical to that given in reference [1]. The function 
Fk, y) = 1/(R + 1) 


dicates the over-all effectiveness of the entire surface, secondary 


is therefore the correction factor which in- 


and primary 
The values of the function F(A, ¥ 
shown in Table 1 for aspect ratios of y = 1, 2,4, and 5 for the fin 


have been calculated and 
parameter between 0 and 10. 


Table 1 Correction factors for finite fin parameters, laminar flow 


K y¥2=10 y¥=20 y = 4.0 5.0 
0 0 968 207 0.138 
| 0.993 710 645 
2? 0.996 823 777 
j 0. 998 SON 873 

10 0.999 Q57 O44 


587 
So 
O38 
YO6 
US6 


These data are shown graphically in Fig. 2 


It is interesting to 
observe that, for ordinary fluids, such as air, in a square channel 
with b/w ratio of 20 and the wall conductivity in the order of 20, 
the fin parameter would assume a value of the order of 25. For 
liquid water the A-value would be of the order of 5. For these 
values, the correction factors are close to 1, representing a small 
reduction of the Nusselt number. For liquid metals, however, for 
which the thermal conductivity is comparable to that of the wall, 
the resulting fin parameter may be of the order of unity. Hence a 
reduction in the heat-transfer coefficient may be quite appreciable 4 


Influence of Secondary Surfaces on Limiting Nusselt Number 
With Slug Flow and Constant Axial Temperature Gradient 


Although slug flow does not exist in nature 
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and is only valid as a result of the nonviscous assumption, calcula- 
tions of the Nusselt number based on the slug velocity distribu- 
tion have been found to be fruitful in estimating the contribution 
of molecular conduction in turbulent heat transfer for liquid 
metals. In the case of a rectangular duct such as shown in Fig. 1, 
the interconnecting ribs may cause circumferential variations in 
temperature and the ordinary boundary condition of a uniform 
wall temperature cannot be maintained. Again in the present 
problem, the use of a fin effectiveness will not lead to the correct 
estimate as the result of this section will demonstrate. 

Under the same assumptions as for the preceding section and 
with an additional condition of uniform velocity profile, the 
problem is simply to solve, 


V0 = Cw 


for which the boundary conditions along the primary 


secondary surfaces are 


6 


O78 


dy 


The solution is obtained by including 4, 


Thus 


part of the general solution 


nwr 


nT ij 
cosh sin 
} 


} ) 


nT 
n*(m*y? + n? sinh * 
«y) 


wk, 
nw 


\ bk r 


Employing the same procedure for evaluating the Nusselt num 


ber, the following expression is obtained: 


Nu 
‘ 


T 


Cc 


: 9 > 
m*n?(m*y? + n?) 


I quations (25 26) can be reart inged to read 


Nu = (Nu,w.)F,(K, ¥) (27 


where Nu... are the slug-flow Nusselt numbers for the conditions 
of uniform wall temperature with respect to peripheral positions. 
Their values for various aspect ratios can be computed from equa- 
tion (25) by letting R = 0, corresponding to infinite thermal con- 
ductivity for the fin material. These values are also contained in 
reference [2 A correction factor F, for finite fin parameters is 
similarly defined as 


(28 
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Fig. 3 Correction factor for finite fin conductivity for slug flow 


a Parameter 


and the numerical results tabulated in Table 2. The values also 
are plotted in Fig. 3. Again, a typical example involving a liquid 
metal as the heat-transfer fluid shows a A-value of the order of 
unity, for which the correction factor to allow for fin inefficiency 


may be very significant 


Table 2 Correction factors for finite fin parameters, slug flow 


k= (52) 20s s=10 5-20 


672 0 421 0 
RO5 0 764 
938 0 852 
955 0 915 
985 0 962 


163+ 
0.531 
0.674 
0 JOS 
0 905 


Combined Forced and Natural Convection for Laminar Flow 
in Vertical Channels With Internal Heat Generation 


Analysis. The main assumptions in this analysis are that the 
flow velocity is fully established and that the heat input from the 
The heat- 


input distribution varies in the transverse directions in such a 


boundary to the fluid is constant in the flow direction 


manner that the wall temperature becomes uniform in the trans- 
verse plane. Because of the uniform heat input with respect to the 
flow direction, the temperatures of the walls and the fluids in fully 
established regions become linear along the flow direction. Fluid 
properties are assumed constant except when the body force is con- 
sidered. Heat dissipation due to the viscosity effect is assumed 
negligible. The fluid may contain heat sources which are assumed 
to have arbitrary but identical distributions in all transverse 
planes. The flow is in the vertical upward direction, along the 
positive 2-axis 


Under the foregoing restrictions, the eq iations of motion 


temperature distribution are 
1 op ; 
7 Fs 
Oz 


ot 
—-@q 30 


and 


Vu = 


k VW%t = puC, 
v p " = 


The co-ordinate axes shall be those shown in Fig. 4. In equation 
(29) the term F represents the body force in the positive z-direc- 
tion 

To express the body force per unit volume, the fluid is assumed 
to obey the following linear law of density variation with tem- 
perature; 
t — t)B) 


p = p[l — (31) 
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Fig.4 Flow channel configuration for combined convection 
where £ is the temperature coefficient for density. The validity of 
this relationship is of course dependent upon the fact that no phase 
change of the fluid occurs, and that the temperature variation 
is small enough to make the temperature coefficient 8 applicable 
within the scope of the investigation. 
ture-difference function 6 defined as 


9 = 


Introducing the tempera- 


(32) 


(4) 


where f) represents the wall temperature at the reference station 

of z = Oand C, being the temperature gradient in the axial direc- 

tion, the state equation (31) is rearranged to read 
p = po[l — B(O+ Ciz)] (33) 

Noting that 

— (82) 


Pwo = poll (34) 


and assuming (C28 is less than unity within the 
problem, equation (33) can be approximated by 


scope of the 


p = p,(1 — BA) (35) 
Hence the body force in the positive z-direction is 
Fs, = —pg = —py(l — BA)g 


and the equations of motion and temperature are now 


1 [ oP ie 
aloe Tt Pil 88) 


Vv = [Q/k 


Vu 


[C.puC,,/k,] 


referring to the co-ordinate system shown in Fig. 4. 
In order to express results in general terms, the following di- 
mensionless parameters are used: 


uC,C,D? 
x, = 2/D; o = 0 k a; "| 
/ ky 
Op 
dD; L = »*| ps | Um bd 
Oz / 


F = Q/[pC,Ciu,] 
R, = [p*gC BD*4 kp} 


where D is the equivalent diameter; u,,, the mean fluid velocity 


Equations (37) and (38) are therefore simplified to read 


0* 
Ox;? 


o*h 
Oy? 


o* O*%” 
+ Rad 


. 


Ox;? oy,’ 


The boundary conditions are simply 
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=0, (a@/D); y=0, (b/D) (42) 


Solution: Since there are two simultaneous differential equa- 
tions, the method of undetermined coefficients using double 
Fourier series is convenient. The pressure-drop parameter L is 
independent of x and y and is therefore expressible in terms of the 
double Fourier series; hence 


“2 y 
2 
* I, 3 


m,n=1,3,5. 


Max, nT) 


. sin 
(a/D) (b/D) 


a b 
for x, >0 and 
D 


>m> 0. 


The heat generation function F, if its distribution across the 
section is known, can be expressed in the following manner: 


) 


MIX, 
(a/D) 


NW i; 


F sin 
(b D) 


sin 


(44) 


mn 


m,n=1,2,3,4 


where the Fourier coefficients can be obtained easily by the usual 
In the case of uniform distribu- 
tion of heat source, the function F is 


rule or some tabulation method. 


16F - 1 


7? 


MAX, ni; 


sin (45) 
(b/D) 


- sin - F 
m,n=1,3,5..™"% (a/2) 
Finally, the functions v and @ are assumed to have the following 
forms: 


ny 


in (46) 
(b/D) 


) 
. , MAX, 
> sin — 8 
(a/D) 


-1,2,3 

_ mr, | ony _ 
> ®,,,, sin sin 47) 
T 23 (a/D) (b/D) 


which satisfy the boundary conditions prescribed. 
Substitution of equations (43), (45), (46), and (47) into equa- 
tions (40) and (41) and equating coefficients of sin mmz,/(a/D) sin 


= 


m,n= 


ntry:/(b/D) terms, the following relations are obtained: 


(48) 


_=0 for m,n = even integers 


n°? 
(b/D)?* 
nmr? |?) 
(b D)? ( 
mr? 


F nr? 
( 16 (a/D)? (b/D)? 
~ \ar2mn j | mr? 
( 


| nor? 2 | 
, a/D)? (b/D)? | § 


{ 
The solution defined by equations (48), (49), and (50) satisfies 
the governing equations as well as the boundary conditions. 


mr? 
FR, + L | . 
(a/D)? 


( 16 ) 
7 ?mn mr? 
dr, + | 
( 


a/D)* 


for m,n = odd integers. 


Results and Discussions 


The Pressure-Drop Parameter. For prescribed values of the tem- 
perature gradient, heat-source intensity, and pumping rate, the 
Thus one may regard the 
heat- generation function F and the Rayleigh number as the input 
functions with L as the resulting function. 


values of F and R,, the resulting L-value is defined by the follow- 


pressure gradient is a fixed quantity. 
For arbitrarily given 
ing relation: 
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(51) 


f, vdrdy = f, drdy 


Substitution of equations (43), (45), (46), (47), (49), and (50) into 
the integral (51) and performing the indicated operation, the fol- 
lowing algebraic relation is obtained: 


L = 


64 l 
FR 
wr ‘ 2d . j mr? n°? || 
5 min? yRa 4 4 


mnewl.; . 
(a/D)? (b/D)? |f 
mr? n*r* 
(a/D)? (b/D)? 


n*x? }2) 
+ ‘ 
(b/D)* { 


where the indexes assume values of positive odd integers 


= m*r* 
3.5. m%? Ry + 
' (a/D?*) 


Inspection of equation (52) shows that the pressure-drop 
parameter is a linearly decreasing function of the heat-generation 
parameter F. Heat generation (positive F) reduces the pressure- 
drop requirement and heat absorption (negative F) increases the 
necessary pressure gradient. For an aspect ratio of unity, i.e., 
a = b = D, the following expressions have been obtained from 


equation (52) by assigning various values to Ry. The linear 


equations containing F and L follow: 
for Ry, 0 
1707 F for Ry, w* 
14.77 F for Ry, = 108¢ 
357 F for Ry = 100r'* | 


The foregoing parametric representations are graphically shown 


in Fig 5 





e Drop Peremcte- { 


Fig. 5 Nusselt number and pressure-drop parameter at various Rayleigh 
numbers. Square tube only. 


The Nusselt Number. 
the usual manner is given by the following expression: 


The Nusselt number when evaluated in 


Nu = (F — 1)/(—4@,,,) (54) 


where @,,, is the dimensionless mixed mean temperature difference 
between the fluid and the wall at the same location, i.e. 


j 
?... = f, gud A / f, aA 


for which the following expression is obtained: 
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64 l 


rr 


m,n=1,3,5 


mr? n*x? t 
a 
a/D)? (b/D)? { 


J R, + mr? 
= (a/D)? 


m*n? 
n*r? 


{1 m*r* 4 PR } 
“| (@/D)? (b/D)? “( 


n>? 3)? 
(b/D)? f 


For an aspect ratio of 1 and a heat generation parameter of zero 
and for selected R,-values, the Nusselt-number values are tabu- 
lated in Table 3, where the Nu at Ry = 0 is identical to the case 
of pure laminar forced convection. The values in Table 3 are 
plotted in Fig. 5. 


(56) 


Table 3 Infivence of free-convection on forced convection Nusselt 


number without heat source, square tubes 


Ra Nu 
0 3.61 
x 3.69 
10x 4.27 
100x* 9.46 


When dealing with the combined effects 
of free and forced convection of the present problem, it is possible 
to have a zero Rayleigh number by either one of the two conditions; 
8 = OandC, = 0. The former corresponds to an incompressible 
fluid and consequently no free-convection effects can occur. The 
latter case denotes zero axial temperature gradient. Under this 
condition, free-convective current can exist if the heat generation 
Q is not zero. 
perature gradient, the problem leads to a trivial solution. 
Inspection of the variables defined by equation (39), however, 
indicates that they do not admit the case of C; = 0 as a possi- 
bility since, for C; = 0, F and @ become infinitely large. For the 
case of pure free convection with zero axial temperature gradient, 


Zero Rayleigh Number. 


Of course, with zero heat generation and zero tem- 


the governing equations may be directly obtained from equations 
(37) and (38) as 
(57) 


V 2u 
VO= 


— pg80/ pu 
—Q/k, (58 


for which the solutions are simply 


@ 


16 1 


r? 2. . mr? nr? 
ip BPS HS -> wae + 
a? b? 


_ marr 
sin 


p.gBQ 16 
kU, 7? 


1 _ mmr nw : 
sin sin (60 


2. a mr? n*n?\? a l 
oy SS See mn + 
2 b? 


a* 


For a given Q, equation (59) defines the temperature distribution 
and equation (60) defines the induced flow rates. It is of interest 
to note that the specific heat has no effect on the present prob- 
lem. 

Negative Rayleigh number. The Rayleigh number is defined as 


Ra = [p'gC,CiBD*/k ju] 


A negative value for the Rayleigh number would result if the 
temperature gradient should become negative along the flow 
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direction. Under this condition, the flow field is referred to as 
heating from below which leads to the Rayleigh’s instability as 
discussed in reference [5]. 

Inspection of equations (46) through (50) indicates that for all 
positive values of R,, the solutions are valid. Whereas for nega- 
tive R,, there are numerous R,-values for which the solutions as 
given by equations (46) through (50) break down 


when 
R, + mr? 
. (a/D)? 


This occurs 


nx? "]2 
= () 
t (b/D)? 


167 


—-R, = 
_ (1 + y)* 


(m*y? + n?)? (61) 


For a square tube, the critical values of the Rayleigh number are 
4m‘, 1001‘, 6767. . It is difficult within the scope of the 
present paper to ascertain the implications of these critical values 
since the equations employed herein do not apply to this type of 
nonlinear study. 

Insulated Walls. The case of zero heat transfer from the walls can 
be represented by F = 1. For nonuniform heat generation, the 


equivalent condition will be fires = 2 dA. The Nusselt num- 


ber as defined by equation (54) will be zero. In this manner, one 


has the condition of “insulated” walls. Inspection of the tem- 
perature expressions (47) and (49) indicates, however, that the 
temperature gradient on the boundaries is not identically zero for 
the case of F = 1, In fact, the temperature gradient changes 
sign along the boundaries. This indicates a reversal in the heat- 
flow direction, resulting in zero net heat transfer while the walls 
are maintained at a uniform temperature without peripheral 
variation 

From a consideration of heat conduction, existence of heat-flow 
reversals requires that the thermal conductivity of the materials 
of the walls be infinitely large so that heat may be conveyed along 
the secondary surface without change in its temperature. On the 
other hand, the case of truly insulated walls is one for which the 
fluid-temperature gradient in the normal direction is identically 
zero on the walls. It can be shown that the condition of identi- 
cally zero normal gradient is not compatible with the requirement 
of uniform wall temperature with respect to peripheral locations. 
If the walls are truly insulated which can be represented by a zero 
wall conductivity, then the wall temperature must vary along the 
boundaries in the peripheral directions. The solution to this 
problem can be obtained by some approximate procedures such 
as the collocation method. 

Finally, it may be stated that, if the wall conductivity is neither 
zero nor infinity, the resulting temperature distribution for the 
case of zero heat transfer must lie between the two extreme cases 
discussed previously. It is apparent that in a circular pipe with 
symmetry, an identically zero temperature gradient is the only con- 
dition for zero heat transfer, while the peripheral variation of the 
wall temperature is zero, by reason of circular symmetry 


Conclusions 


The following conclusions may be drawn: 


1 For laminar flow in channels with secondary surfaces, the 
over-all limiting Nusselt number or heat-transfer coefficient for 
fully established flow conditions is greatly influenced by the fin 
parameter (wk,,/bk,), as shown in Figs. 2 and 3 for laminar flow 
and slug flow, respectively. 


2 For small values of the fin parameter, the peripheral varia- 
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tion of the wall temperature is quite appreciable. For rectangular 
channels, the controlling factors are the aspect ratio and the fin 
parameter. 

3 For heat transfer in vertical channels with combined-con- 
vection effects and heat generation within the fluid, the general 
trend is comparable to that in a circular tube investigated by 
Hallman and reported in reference [3]. 
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ADDENDUM 


After the manuscript of this paper was reviewed, the author 
received comments relative to the subject matter. Aside from 
editorial comments with respect to typographical errors and some 
obscure parts of the paper which the author has corrected and 
clarified, the following points are of significance: 


1 That, in the case of three-dimensional flow, the assumption 
of fully established flow does not necessarily lead to equations 
(2), (29), and (30) of this paper. 

2 That the case of zero Rayleigh number is not necessarily 
identical to pure convection flow. 

3 That equations (12), (13), and (14) can be solved by re- 
laxation methods. 


With respect to the first point, the mathematical definition of 
fully established flow is simply 0u/0z = 0. A recent publication 
by Maslen [8] gave proof that the assumption of fully developed 
flow implies zero velocities transverse to the flow axis for incom- 
pressible viscous flow in laminar motion. 

Regarding the case of zero Rayleigh number, the author has 
added a paragraph explaining the possibility of pure free-con- 
vective current with zero axial temperature gradient and heat 
generation. As to the third comment on relaxation methods, the 
author is aware of the availability of approximation methods. 
Since the exact solutions were obtained, no effort was made to use 
the relaxation methods in the present application. However, it 
may be mentioned that the solution of equations by relaxation 
methods might be indeed convenient in solving the problem of 
insulated walls as mentioned in the paper. Finally, the author 
would like to thank the reviewers for their comments which ma- 
terially contributed to the completeness of the paper. 
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DISCUSSION 
P.C. Lu’ 


The writer is particularly interested in the third problem dealt 
with in this paper. A thesis‘ has recently been submitted to 
Kansas State College, in which the writer solved the same prob- 
lem by employing the finite Fourier double sine transform. The 
results check with those presented in this paper, although dif- 
ferent parameters are used. Basically, the finite Fourier trans- 
form and the classical Fourier method are of the same origin: but 
usually the method of Fourier transform is considered to be more 
straightforward in procedure 


In the 


problem when the channel is of ring or full-sector cross section by 


thesis mentioned, the writer also was able to solve the 


using finite Fourier single sine transform 
With reference both to this paper and a previous paper® by 
the author, the writer would like to point out that an analogy 
exists between some special cases of the problem under discussion 
and the deflection of plates lving on an elastic foundation 
When ©, 0, 


exists just 


an analogy to the ordinary thin-plate theory 
8 In pure forced convection Results for triangular 


pipes have been obtained by observing this analogy. 
*(raduate Research Assistant, Department of Mechanical Engi- 
Kansas State College, Manhattan, Kan.; now, Graduate 
Department of Mechanical Engineering, Case Institute of 
Technology, Cleveland, Ohio. 
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neering, 
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Lu, “A Theoretical Investigation of Combined 
Forced-Convection Heat-Generating Laminar 
Pipes With Prescribed Wall Temperatures," 
State College, Manhattan, Kan. 

*S. M. Marco and L. 8S. Han, “A Note on the Limiting Laminar 
Nusselt Number in Ducts With Constant Temperature Gradient by 
Analogy to ASME, vol. 77, 1955, pp. 


Free and 
Inside Vertical 
thesis 


Flow 
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Author’s Closure 


The author is pleased to learn that the problems for combined 
convection in channels of ring or full-sector cross sections have 
been treated in Mr. Lu’s thesis. The uses of the classical Fourier 
method or of the Fourier transforms are extensively discussed in 
mathematical literature and is not the center of discussion here. 
It is, however, not impertinent here to recall the basic structure 
of the equations governing the third problem dealt with in the 
present paper. If one combines equations (40) and (41), assum- 
ing uniform heat generation, the following equation in terms of @ 
is obtained: 


V‘o T Ro = - L 


with the boundary conditions 


o=0 


V%~ =F 


To these equations, an analogous problem in the thin plate theory 
is obvious; namely, a plate resting on an elastic foundation with 
elastic constant R,, with zero deflection on the boundaries and 
acted upon by a constant bending moment (in the cass 


of plates 
of polygonal shapes). 


Solutions in this case can be best handled 
by the use of the energy method. 

The author would like to thank Mr. Lu for his discussion and 
Professor Hoagland of the Massachusetts Institute of Tech- 
nology for pointing out a misprint which was the result of an 
oversight in the original analysis. Fig. 2 and Table 1 have 
accordingly revised, in the final form. 


been 


Transactions of the ASME 





Heat Transfer to Water in 


S. LEVY: 
R. A. FULLER? 


Thin Rectangular Channels 


Heat-transfer coefficients are presented for water flowing vertically in thin rectangular 


R. 0. NIEMI: 
periphery. 


channels (0.1 X 2.5 in.) 18 and 36 in. long and heated electrically around the entire 
The range of vartables covered is: 65 to 200 psia pressure, 90 to 200 F sub- 
cooling, and 4 to 50 fps water velocity. 
along the narrow and wide faces of the rectangular test section. 
reported with steam blanketing occurring first at the corner of the test section. 


Heat-transfer correlations are given for data 
Burnout data also are 
The 


proposed correlating equation, valid at the narrow face of the test section, gives values 
considerably lower than those obtained in a circular pipe. 


= NEUTRON FLUXES, typical of nuclear test reac- 
tors, require unusually large specific powers; i.e., large heat pro- 
This requirement is 
often met by utilizing flat-plate fuel elements of the type shown 
in Fig. 1. The fuel elements consist of several thin plates of 
uranium-aluminum alloy clad on both sides with aluminum. The 
plates, held in position by two side plates, are normally cooled by 
downward flow of water at rather low temperature and pressure. 
The heat-transfer problems associated with this kind of fuel ele- 
ments are: 


ductions per unit volume of reactor core. 


1 Use of thin rectangular channels of bigh aspect ratio, meas- 
uring about 0.1 
2 Continuous variation of bulk-coolant properties in the flow 


direction. 


x 2.5 in. 


3 High film-temperature drop from surface of fuel plate to 
bulk of the coolant and, therefore, large changes of water proper- 
ties close to the heated plate. 

4 Entrance effects where for a short distance beyond the inlet 
of the fuel element the velocity and temperature profiles are not 
fully develope d. 
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5 Possibility of “burnout” or steam blanketing of the fuel ele- 
ment during transient power increases or flow reductions. 


Yet, as pointed out in the literature survey that follows, data 
are often lacking when two or three of the foregoing conditions 
exist simultaneously. In particular, rectangular geometries have 
not been studied for circumstances encountered in a nuclear test 
reactor. For this reason, late in 1955 the Atomic Power Equip- 
ment Department of the authors’ company initiated an extensive 
heat-transfer program. A heat-transfer and burnout loop was 
especially erected to study flat-plate fuel elements of the type 
used in test reactors. This paper covers the experimental results 
obtained. Following a summary review of available data, the test 
equipment is described and the test data are presented and dis- 
cussed. 


MEAT mo bate FUEL PLATE 


SIDE PLATE 
li 


Fig. 1 Typical flat-plate fuel element 





heat-transfer area, sq ft 
flow area, ft? 
specific heat, Btu/lb deg F 
hydraulic diameter, ft 
deflection, in. 
mass flow rate per unit area, 
lb/hr sq ft 
heat-transfer coefficient, 
Btu/hr sq ft deg F 
current in test section, amp 
= thermal conductivity, 
Btu/hr ft deg F 
= Nusselt number, nondimen- 
sional, (hD/k) 
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Prandtl number, nondimen- 
sional, (yc/k) 

Reynolds number, nondimen- 
sional, (@D/u 

average test section perime- 
ter, ft 

heat generation in test sec- 
tion, Btu per hr 

burnout heat flux, Btu/hr sq 
ft 

resistivity of test 
ohm-ft 

bulk temperature, deg F 


section, 


test-section, thermocouple- 


temperature reading, deg F 


Subscripts 


outlet-water 
deg F 


saturation 


temperature 


temperature at 
outlet of test section, deg 
FP 

surface temperature of test 
section, deg F 

water velocity, fps 

total deflection, in 

absolute viscosity, lb/hr ft 


bulk 
film 


wall 
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Literature Survey 


The problem of heating and cooling fluids inside tubes has re- 
ceived considerable theoretical and experimental attention [1] .* 
Test data for turbulent flow have been correlated in several ref- 


erences. They include the equations of the following: 


1 The Dittus-Boelter [2] type where all physical properties 
are evaluated at the bulk fluid temperature 


hD (22) (my ‘ 
= 0.0243 
k M k 
h G\%* (cp\'? 
.D = 0.0265 (/ ) (*) 
k 7 k 
) 


2 The Colburn [3] form where properties are taken at the so- 


for heating, and 


for cooling 


called film temperature or arithmetic mean of surface and bulk 


temperature 
D DG\°* {cu\ /* 
AD. 0.023 ( > (+ 
ky By, kJ, 


3. The Sieder-Tate [4] kind which utilizes both bulk and sur- 
face properties 


AD DG 0.8 . : . ole 
: 0.027 ( (“ (*) (3) 
k uM k My 


Yet, as pointed out by McAdams [1], equations of the foregoing 
A coefficient of 0.023 has 
Similarly, a coefficient of 
0.023 is recommended [1] for Equation (3) and some data for 
air employ a value of 0.021 [6]. 


form are proving to be inadequate 
been suggested [5] for Equation (1) 


Recent tests [7] with gases would 
place the coefficient in the Colburn equation at 0.019. Even 
when only considering data for heating water, the correlations are 
not always in agreement. Test results for pipe flow at pressures 
from 500 to 2000 psi and with small property variation have 
been correlated by means of the Colburn equation with coefficients 
of 0.023 [8] and 0.019 [9]. Experiments at lower pressures and 
with large property variation include the data of Kreith and 
Summerfield [10] for circular tubes which agree with Equation 
(3); those of McAdams, et al. [11] for flow in an annulus which 
satisfy the Sieder-Tate equation with a coefficient of 0.026; and 
those of Stein, et al. [12] for an annular geometry which are best 
correlated by a Colburn equation with a coefficient of 0.020 [13]. 
Equivalent correlations for turbulent flow in rectangular channels 
are practically unavailable. Besides the data of Bailey and Cope 
14) which are somewhat uncertain due to a nonadiabatic mixing 
device at the outlet, there are the test results of Kays [15] for 
heating air in a rectangular passage with an aspect ratio of 5.85 
Kays’ values are for Reynolds numbers of less than 10,000 and 
fall 24 per cent below those predicted from Equation (2). Also, 
Washington and Marks [16] reported heat-transfer coefficients 
in rectangular channels substantially below the predictions of 


(v2 2 


Equations (1 , or (3 


Test data for Reynolds numbers less 
than 10,000 may be in a transitional region and the turbulent- 
flow, heat-transfer-correlating equations may not be applicable 

Similarly, entrance effects in rectangular channels have not 
been studied previously. Results, however, can be found for cir- 
cular tubes. Recent test results for water flow in pipes have been 
presented by Hartnett [17] and they agree with the analysis of 
Deissler [18]. 

Heat-transfer burnout data approaching the range of interest 


in test reactors have been obtained by McAdams, et al. [11] and 


> Numbers in brackets designate 


130 


References at end of paper. 


MAY 1959 


Gunther [19]. McAdams, et al. [11] have reported burnout-flux 
values for degassed distilled water flowing upwards in a vertical 
annulus containing a central tube heated electrically. Pressures 
ranged from 30 to 90 psi, water velocities from 1 to 12 fps, sub- 
cooling from 20 to 100 deg F. The diameter of the annulus was 
not varied in the tests and the data apply to a 0.25-in. central 
heater in a 0.77-in. glass tube. Gunther [19] carried out similar 
measurements in a horizontal glass passage '/, in. high, */\¢ in. 
wide, and 6 in. long containing a horizontal electrically heated 
metal strip at the axis. Water velocities varied from 5 to 40 fps, 
subcooling from 22 to 282 deg F, pressures from 15 to 160 psia. 
These tests, however, do not cover the case of downward flow of 
water in very thin rectangular channels with high heat production 
in the corners of the rectangular section (to simulate gamma and 
neutron heat generation in the side plates of Fig. | 


Description of Equipment 


Test Loop. The heat-transfer burnout loop is shown in Fig. 2. It 
is built entirely of stainless steel and designed for 1000 psi. A 
schematic representation of the experimental system is shown in 


Fig. 3. It consists of: 


i Primary loop, where water is circulated in a closed system 
by means of a centrifugal pump through a preheater, along the 
steam pressurizer connection, upward or downward through two 
parallel test sections, and on to a vertical or horizontal heat ex- 
changer. 

2 Demineralized water system, utilizing a demineralizer and 
associated cooling equipment to yield the desired water purit) 

3 Heat-exchange system, consisting of a horizontal heat ex- 
changer for low temperature and low pressure, and a vertical heat 
exchanger for high temperatures and net steam generation 

t Test sections, where heat is generated by electrical heating 
of the test section walls up to 9 ft in length 

5 Power supply, using six General Electric d-c are welders 
rated at 40 volts and supplying 1540 amp each to the test section. 

6 Instrumentation and controls, consisting of temperature, 
pressure, flow, current, and voltage-measuring equipment 


—_— 4@ 
\ 


Fig. 2 Heat-transfer and burnout loop 
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Fig. 3 Flow diagram (heat-transfer and burnout loop) 


Test Sections. A typical test section is shown in Fig. 4. It con- 
sists of a rectangular channel about 0.100 in. « 2.500 in. fabri- 
cated of stainless steel. A copper flange is attached to each end of 
the test section and it is connected to the current-carrying bus 
bars by means of flexible connections. Two side plates spanning 
the entire length and width of the test section are provided to re- 
sist internal pressure forces. The side plate as well as the rest of 
the loop are electrically insuiated from the test section and copper 
flanges. 
The heating-element size varied slightly from one test section 
to another. The exact dimensions of test sections 4 and 11 are 
shown in Table 1. Test section 5, also listed in Table 1, is a dis- 
torted test section used to find the effects of dimension variation 
in rectangular geometries. Its dimensions are given in a later 
part of the paper. Test sections used for burnout tests measured 
0.118 in. X 2.405 in. They were also 18 in. long and 0.028 in. 
thick. On the other hand, test sections 4, 5, and 11 were 36 in 
long and 0.056 in. thick. 
Twenty-eight thermocouples were installed at various posi- 
tions on the test section. The test sections also were instru- 
mented with two pressure taps at each end of the heating element. 
The thermocouples, insulated from the test section by 0.0015-in 
thick mica tape, were distributed lengthwise, and around the pe- 
riphery of the rectangular cross section. Their exact location is 
given in Fig. 5 which illustrates a typical set of thermocouple meas- Fig. 4 Heat-transfer test section 
urements. It should be noted that thermocouples labeled C, 
and C, were placed at the centers of the narrow sides of the test the 2.500-in. faces of the heating element. Thermocouples Mr 
section, while thermocouples Mr, Mz, §;, and S: were located on and My, are on the center line of these two faces, and thermo- 
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Table 1 


Run 


no, 


Test section Heat flux Q/A, 
no. Btu/hr sq ft 
4 265000 
550000 
267000 
572000 
268000 
574000 
269000 
579000 
173000 
162000 
, 010000 
990000 
983000 


fps 
70 
10 
70 
70 
30 
20 
5.40 
40 
70 
7.50 
40 
10 
38.380 


0.118 x 2.405 in. 


1] 
0.114 x 


815000 
630000 
530000 
999000 
547000 
440000 
740000 
1040000 
802000 
600000 
447000 
1018000 
27000 
60000 
103000 
154000 
446000 
449000 
331000 
240000 
822000 
244000 
113000 
174000 
65000 
59000 
58000 
107000 
122000 
159000 


10 
2.349 in, 
90 
20 
36.00 
00 
40 
90 
90 


50 
40 
54 
75 
76 
70 
40 
5.55 
62 
9O 
50 
95 
95 
90 
55 
32 
76 
68 
05 
90 


Cl i 1D 
CSN 


so) 


wwwwfenonmnwc 


5 ‘ 650000 
: 308000 
440000 
615000 
302000 
1000000 


For dimensions see 
Fig. 8 


* Upward waterflow 


couples 5, and S, are on either side of the center line My and 
0.735 in. away from it. Bulk temperatures were measured at the 
inlet and outlet of the test section and are assumed to have varied 
linearly from one end to the other. 

The test sections were equipped with a burnout detector. This 
burnout contro] system was capable of turning off the power sup- 
ply before peak fluxes could damage the test section. It consisted 
of balancing the resistance of the exit end against the resistance of 
the rest of the test section. As burnout occurs, the bridge be- 
comes unbalanced and shuts off the electrical power supply. 

Loop Operation. The experimental procedure consisted of the 
following: 

1 Filling the system with demineralized water. 

2 Starting the circulating pump and pressurizing the system to 
the desired pressure 

3 Applying the desired power to the test section after adjust- 
ing the flow and pressure setting. 

4 Regulating the inlet-water temperature by changing the 


secondary flow in the heat exchanger or utilizing the preheater. 
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40 


5.00 


Heat-transfer runs for water flow in thin rect i} 


Water 
velocity, 





Ratio of surface to 
bulk temp, 
T/T», deg F/deg F 
348-1. 282 
715-1 .468 
416-1 .344 
768-1 . 500 
400-1. 345 
895-1 .608 
483-1 .395 
782-1. 560 
468-1 .383 
238-1. 200 
09-1. 542 
08-1. 542 
04-1. 542 


Bulk Reynolds no., 
Nreculk) 
104000—124000 
103000--152000 
104000—125000 
104000-—-15 1000 
93000-—1 12000 
94000-—14 1000 
89000-—1 10000 
79000-— 126000 
41000-56000 
104000—117000 
110000--195000 
109000— 196000 
106000— 190000 


Bulk come, Ts, 

deg 

112-130 
112-152 
112-132 
112-154 
113-134 
112-158 
116-140 
113-168 
113-146 
113-124 
114-187 
114-186 
113-184 


we we ee 


118-212 
126-194 
113-172 
126-209 
97-141 
90-124 
129-184 
102-175 
107-161 
94-133 
86-114 
101-165 
79-92 
79-114 
82-144 
72-97 
91-143 
92-169 
85-138 
78-117 
118-190 
73-127 
68 
70 
66-4 


73000-— 146000 
81000- 135000 
72000- 118000 

110000—200000 
85000—128000 
78000--1 11000 

125000-191000 

100000—186000 

113000—192000 

103000—150000 
93800-—127000 

110000-—194000 
10600-12400 

9000-13300 
8400-17300 
29400-39700 
53200-87900 
36400-72500 
34800-59600 
32000—48800 
95200—168300 
22900—40900 
21600-29400 
22100-33500 
18800-22700 
9300-13700 
7900-12400 
7700-16300 
7300-17900 
7100-21400 


31—1.589 
930-1 .535 
975-1 .587 
175-1.635 
990-1 .632 
869-1. 550 
853-1 .599 
380-1.753 
360-1 .870 
190-1. 829 
000-1 .720 
395-1. 849 
278-1.172 
648-1 .316 
080-1 . 432 
723-1.391 
090-1. 645 
400-1. 658 
180-1 .623 
00-1. 545 
10—1.585 
26-1. 562 
735-1. 465 
03-1 .555 
455-1. 345 
612-1. 302 
701-1 
195-1 
315-1 
645-1 


99 
33 
? 


3 
4. 
42 
37k 


Dw NM RK eK me WN NNN NRK NRK NNN NNR ee bo bo 


105000-157000 
101000—125000 
96000-—13 1000 
105000—157000 
83000—107000 
107000—193000 


605 
485 
390 


985-1 
630-1 
492-1] 
630-1 .445 
385-1 .275 


116-191 060-1.345 


BND bet eet et et et 


5 Recording surface temperature, system pressure, pressure 
difference across the flow orifice, water inlet and outlet tempera- 
tures, current and voltage to the test section, and pressure drop 
along the test section. 

6 Reducing the data. The inside wall temperature 7, is 
computed first from the measured value Tp at the outer surface 
according to the expression [10] 

3.4137? r 
To -T. = a 

2p* i 


= current in test section, amp 

= resistivity of test section, ohm-ft 
k = thermal conductivity, Btu/hr-ft deg F 
p = average test-section perimeter, ft 


where 


In the present case there is no need to correct for the variation of 
electrical resistivity and thermal conductivity with temperature 
since their ratio for 3048S, used in the test section, remains prac- 
tically constant and varies only 0.33 per cent, well within the 
range of the accuracy of the data. 
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Table | (cont'd) Next, the mean temperature of each segment of the test section 


is determined and the corresponding heat output obtained from 
the measured current and an extrapolation of the measured re- 
sistance at room temperature. The fluid temperature at any 
29 115 | position along the test section is derived by integrating the heat 
34 | 116 generated to the desired point. A knowledge of the heat flux, 

] 117 | bulk and.surface temperatures in turn determines a heat-transfer 
lie t coefficient. 


tia | = i 
120 h = (Q/A)/(7, — T;) 


121 From the Sieder-Tate equation 


Heat balance, per [ RUN | symeBot | 


cent 


Run 8(0.3) 
9(0.6) 
10(0.1) 
11(0.5) 
12(0. 43) 
13(12.5)* 
14(14.4)* 


15(6.2) 
16(22.2)* 
17(4.1) 
18(1.3) 
19(1.9) 
20(1.4) 





Error Analysis. The heat-transfer coefficient is calculated from 


122 | : h = 0.023 (k/D) (DV/j)®* (cu/k)'/* (s/w) 
123 
124 | 
125 
126 | (i) Since the measured and calculated heat-transfer coefficients 
i27 | a differ, the factors involved in the experiments and in the correlat- 
> ing equation must be examined. In the experimental heat- 
126 a transfer coefficient, the factors are: 
129 


30 | . 
i31 r 


os h = 0.023 k (p/p)* (cu/k)'/* (u/pw)-'* (V°-*/D9-*) 
* Flow based 
on electrical Q. 


Heat flux, Q/A 
Surface temperature, 7’, 


eeoour+ie Piolie de BPO xX OP 














Bulk temperature, 7, 
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Fig. 5 Typical set of thermocouple measurements for test section No. 4, run No. 11 
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The heat flux Q/A is determined from the heat input Q and 
the heat-transfer area A. The heat-transfer area can be meas- 
ured and calculated to a high degree of accuracy. An error of 
0.020 in. in a 5-in. perimeter corresponds to an error of 0.4 per 
cent. The heat input was determined from electrical measure- 
ments of the voltage and current with instruments of '/, per cent 
of full-scale deflection error. Parallax reading errors were 
minimized by the usual reflecting mirror in the instrument. At 
a half-seale reading the instrument errors would be '/, per cent 
and a reading error of 1 per cent give an error of 1'/, per cent. 
Non- 
uniformity of thickness of the test section would yield variations 
in /*2 heating. Commercially drawn tubing is uniform to within 


1 mil, but the forming operation used to make the test section 


This error, however, would vary from reading to reading. 


probably produced a nonuniform thickness at the corners. Or- 
dinary visual inspection did not yield any detectable thickness 
changes 

The temperature measurements are expected to be the critical 
ones. The factors involved in the temperature measurements will 
be considered 


1 Thermocouples. The thermocouples were made by the 
General Electric Engineering Laboratory using a helium are- 
welding procedure. These thermocouples were calibrated to an 
accuracy of | deg to 400 F 

2 Error due to Conduction Loss. There is a 1'/-mil thick- 
ness of mica electrical insulating tape between the outside surface 
of the test section and the thermocouples. Heat conduction due 
to external losses would give a temperature drop across the mica 
tape. To minimize the drop, a 2-in-thick block of magnesia in- 
sulation was installed over the thermocouples. It should be em- 
phasized that the temperature error introduced by conduction 
losses reduces the temperature and yields increased values of the 
heat-transfer coefficient. Radiation losses enter only in increas- 
ing the insulation exterior film heat-transfer coefficient and, since 
the insulation is so thick, radiation effects are small 

3 Thermocouple System. An ice-bath cold junction was used 
test section No. 4 A Leeds and 


Northrup precision laboratory-type potentiometer made specifi- 


for the series of tests with 


cally for thermocouple installations was used. The instrument it- 


self is accurate to 0.1 per cent. Later, an automatic Brown re- 
corder was substituted and temperatures were recorded graphi- 
cally with an automatic plotter. The data duplicated the ice-bath 
precision-potentiometer arrangement and further test data were 
The guaranteed to '/, 
per cent of full-scale deflection Up to 600 F the error would be 
For a 50-F temperature differential this would cor- 
respond to a 3 per cent error 


recorded graphically recorders were 
| deg F 
Higher temperature differentials 
would have a correspondingly smaller percentage error. 

t Wall Temperature-Drop Correction. Since the 


wall temperature was measured, the temperature drop across the 


outside 


test section wall must be known. The correction was given by 


3.413 [7 
2p*h 
The J? error is about twice the error in current J and therefore is 
about | per cent maximum 
O.8 per cent The ratio (r/k 


The maximum error in p? would be 
is determined from other experi- 
menters’ data 


The error in thermal conductivity is the con- 


trolling error (n estimate of the accuracy of (r/k) is 1 per cent 
The maximum error would then be about 2.8 per cent. A wall- 
temperature drop of 100 F would yield a 3-deg-F error. 

Summory. The errors in determining the heat-transfer coefficient 


are: 


= error in reading volts and amperes = 1.5 + 1.5 


(a) (Q A 
= 3 per cent 
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(b) Surface temperature, 7), 


(1) thermocouple, 1 deg 

(2) conduction loss, '/; deg 

(3) thermocouple recorder, 1'/2 deg 

(4) temperature drop through wall, 3 deg for 100-F dif- 
ferential. 


Error for 100-deg-F wall conduction drop = 6 deg. 

The surface temperature rise is of the order of 200 F or 
greater when the wall correction is 100 deg. The error is 3 per cent. 

The total cumulated error could be 6 per cent. 

In the calculated heat-transfer coefficient, the factors involved 


are: 


(a) k(p/p)*, (cu/k)'/*, (w/pw)"4, all dependent on tempera- 
ture. The values are determined from other experimenters’ 
Accuracy estimated at greater than 1 per cent 

(b) V, velocity. A calibrated orifice was used to measure 
flow. The initial accuracy was about '/, of 1 per cent but with 
operation may decrease in accuracy to about 1 percent. A read- 
ing error of 0.05 in. at a 30-in. reading corresponds to an error of 
'/, of 1 per cent. 


values. 


The flow-measurement accuracy is between 1 
and 2 per cent. 


Now, the velocity is determined in the test section from 
W = (pAV )oritice = (PA V test section 


Errors in determining the flow area would affect the velocity. 
The heat-transfer coefficient will have 0.8 the error in the velocity 

(c) D, hydraulic diameter. The hydraulic diameter could be 
determined to about 2 mils accuracy in the unpressurized state. 
The 
test section was backed with a '/;.-in-thick hard-asbestos gasket, 
which was backed with a steel supporting section about 
thick Compression 
of the backing gasket and deflection of the support plate would 
allow the channel spacing to change. 


The channel size was not measured in the pressurized state 


in. 
The whole assembly was bolted together 


The gasket compression can be estimated. For compressed 
asbestos, a 100-psi pressure will give a 3 per cent compression 
Deflection y is: y = (0.03) (0.0625) = 0.001875 
Gasket) y = 0.002 in. 
Deflection of backup plate: Assume plate to be simply sup- 


ported at bolting lines. The plate will act as a uniformly loaded, 


: ( 5 Wi 
— maximum) = * 
384/ El 


d= 


simply supported beam. 
/ deflection 


bd 


where = | in. (assumed), 


Then, for 100-psi pressure 


( 5 ) (100 K 3)(3)8 Pe 
= = ().0027 in 
384 1X ('/,)3 


(30 * 108) 
12 


The total of the plate and gasket deflections is 
Y, = Y¥(total) = 2[0.002 + 0.0027] = 0.0094 i: 


If the middle third is assumed to have the maximum deflection 
and a linear taper of deflection to zero at the ends, then 


w Y, oe , 
wes * ase ee ee 


2Y, 
Apr = we + w( 3 


Neglecting the slight increase in perimeter due to stretching of the 
test section, 
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4Ap 
Dy = e = 
Pp 


[- 


>) 
i= 


4(we + w(2Y,7/3)) 
2(w + ce) 


+ w(2¥7/3) |, 
wt+e * 


[2.5 & 0.108 + 2.5 (2 x 0.0094/3)] } 
2.5 + 0.108 


Dy = 0.219 in. 


The undeformed test section hydraulic diameter was 0.207 in. 
The test section hydraulic diameter could have changed due to 
deflection of the backing plate and compression of the gasket: 


0.012 


— or +5.8 per cent 
0.207 


The velocity change is inversely proportional to the flow-area 
change. 


l l 
— 100 - 
+100 (4¥ " (4 we + w(2Y, ») 
a ia ! 


we 
= —6.7 per cent 


Then, the error in the calculated heat-transfer coefficient due to 
test-section deformation is 


—0.8(6.7) + —(0.2) (5.8) = 6.5 per cent at 100 psi 


The deformation of the test section is the one factor which would 
reduce the heat-transfer coefficient. 


Test Results 


Tables 1 and 2. 
Table 1 lists the tests performed to obtain the heat-transfer co- 


All the heat-transfer runs are summarized in 


efficients while Table 2 shows the runs which determined burnout 
heat flux. 
Over 


40 runs were carried out to measure heat-transfer co- 


efficients in thin rectangular sections. The range of variables was 


Water velocity, fps 
Heat-transfer rate, Btu/hr sq ft 


2.90 — 45.5 fps 
58,000 — 1,040,000 
Bulk water temperature, deg F 66 273 


Ratio of surface to bulk temp, 


deg F/deg F 1.17 2.65 


With the exception of four runs marked with asterisks in Table 
1, all runs were performed with downward water flow. Also, 
boiling was not permitted in any of the tests in Table 1. 

A tota! of 16 burnout runs were performed with downward 
water flow. They cover the following set of conditions: 


6 — 53 
60 — 140 
90 — 200 


Water velocity, fps 
Outlet pressure, psia. 
Subcooling, deg F 


Only 16 runs were completed because rupture of the test section 
happened every two to three runs despite the burnout detector. 
Data of interest involved such high water subcooling and large 
water velocities that excessive thermal shocking of the test section 
resulted after the detector turned off the power supply. 

In all such cases of rupture of the test section it was noted that 
the failure was localized in the corner. The burnout data as re- 
ported here, therefore, apply to steam blanketing at the corner of 
a thin rectangular section with large beat production at the corner. 
This phenomenon was verified further during a voluntary burnout 
test carried out without the detector to check its accuracy 


Discussion of Results—Heat-Transfer Coefficients 


A plot of all the heat-transfer coefficients on the 2.5-in. faces is 
given in Fig. 6, with the exception of those for test section No. 5. 
The data, including those points in the entrance region, are 
plotted in terms of Sieder-Tate parameters and compared to the 
original Sieder-Tate Equation (3). 
that: 


Examination of Fig. 6 reveals 


1 The heat-transfer coefficients are 30 to 45 per cent below 
those predicted from Equation (3), the deviation increasing with 
Reynolds number. A lower coefficient of 0.023 in Equation (3), 
as suggested in reference [1], would still yield values 15 to 30 per 
cent too high. One concludes that rectangular channels for the 


conditions tested exhibit poorer heat-transfer characteristics 


This reduction in heat- 
transfer coefficient is small enough not to lessen the ability of 


than those accepted for circular pipes. 


flat-plate elements to operate at very high heat-transfer rates. 

2 A correlation with a constant exponent on the Reynolds 
numbers is not adequate. Even though simultaneous considera- 
tion of all runs indicates that the Reynolds number exponent may 
be below 0.8, close examination of the runs in Fig. 6 shows that 
the exponent for each run is larger than 0.8 and varies from run to 
run, increasing slowly with reduced heat rates \ major reason 
for the variation in the exponent could be the continuous change 
in fluid properties as the water and surface temperatures vary 
along the heated length 


In particular, the viscosity continually 


Table 2 Burnout data for subcooled water flow in thin rectangular channels 


Outlet 
Run Heat generation pressure, 
no. Q, Btu/hr 
1B 368604 
2B 325000 
3B 432000 
4B Run no. 4 
5B 965000 
6B 1058000 


Heat flux Q 
Btu /hr sq 
600332 74.0 
530000 70.0 
735000 67.0 


temp 7'sat, 
deg F 
306.7 
302.9 
300.0 


psia 


1575000 
1725000 


75.0 341.2 
75.0 307 .6 


7B 879000 
8B 659000 
9B 527000 
10B 398000 
11B 346000 
12B Run no. 12 
13B 501000 
14B 676000 
15B 388000 
16B 525000 


1430000 
1075000 
859000 
648000 
564000 


75.0 
70.0 
70.0 
65.0 
65.0 


307.6 
302.9 
302.9 
298.0 
298.0 


817000 
1100000 
633000 
857000 


120.0 341.2 
120.0 341.2 
115.0 338.1 
138.0 351.9 
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Saturation Outlet bulk Subcooling 
; temp T'our, T'sat Tout, velocity 


Water Q/A 

400,000 + 4800 Q/A 
V, fps (Fon = 7, Fa ~- Tous 
218.0 88 29 0.7270 6768 
204.5 OS 05 0.6075 5386 
173.3 126.7 2.25 0.7290 5801 


deg | deg F 
7 


‘ 
‘ 


no burnout at 50.5 fps, 1,708,000 Btu/sq ft hr 


167.7 173. 2.2 278 9078 
137.6 170 53.0 118 10147 


164 } 204 

164.8 < 9025 
156 21.4 7464 
159 D. 2 5570 


1I31.é ¢ 5470 


8714 
6523 
5492 
1070 
4280 


no burnout at 45 fps, 1,725,000 Btu/sq ft hr 


139.6 201 

148 192. : 
202. § 135.: 
196 155.§ 


5960 
8940 
6031 
7528 


1050 
5730 
1678 
5497 
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decreases in the direction of flow, i.e., the corresponding water- 
velocity distribution changes from inlet to outlet. These 


ondary thermally produced’”’ hydraulic entrance effects could ac- 


“‘sec- 


count not only for the Reynolds number exponent variation with 
heat rate, but also for the spread in the experimental measure- 
ments 


3 Entrance effects are noticeable in each run of Fig. 6. The 
heat-transfer coefficients at the inlet of the test section correspond 
to the lower Reynolds numbers values of each run. They are 
seen first to decrease, then reach a minimum, and finally increase. 
Che thermal-entrance length may be defined as the distance from 
start of heating to the point where the Nusselt number increases 
is the Reynolds number, raised to a fixed power. This thermal- 
entry length is seen to be about one fourth of the test section and 
to vary b | 


it slightly with water velocity. Each of the foregoing 


effects will now be considered in more detail. 


Effect of Rectangular Cross Section. 
efficients 


As already noted, lower co- 
How- 
ever, proper understanding of these reduced heat-transfer condi- 
tions must await measurement of velocity and temperature dis- 
To date, velocity meas- 
urements have been carried out only in triangular ducts in the low 
teynolds number range; less than 25,000. 


vere obtained in the rectangular test sections. 


tribution in rectangular cross section. 


The most recent tests 
of Eckert and Irvine [20] reveal that in a limited range the flow 
may be turbulent near the base of the triangular duct, but laminar 
near the apex of the cross section. A similar pattern of laminar- 
turbulent flow or of low turbulence—high turbulence flow coex- 
istence may be anticipated at the edges of a rectangular cross 


section. The higher surface temperatures measured at the edges 
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coefficients for thin 


rect 





of the present heating elements tend to verify this assumption. 
A plot of heat-transfer coefficients at the center of the short 
side of the rectangular channels is given in Fig. 7. Two values 
are shown for each run corresponding to positions 12 and 24 in 
away from the inlet of the heated section. 
are noted for test section No. 4. 


Higher coefficients 
These may be due to slight dif- 
Test section No. 4 was made up of two L 
strips and had very sharp corners, while test section No. 11 was 
fabricated from a circular tube and had rounded edges. In both 
cases, however, the heat-transfer coefficients are considerably 
lower than those from Equation (3). 


ferences in geometry. 


The values are 60 to 70 
per cent below the Sieder-Tate prediction and about 15 to 30 per 
cent below the data given in Fig. 6. The exact extent of this re- 
gion of reduced heat transfer was not determined. However, Fig 
5 indicates that it must have been confined to the edges of the sec- 
tion since thermocouples §,, S., and Mr, were reading the same 

Another way of tracing the lower coefficients of Fig. 6 to the 
rectangular geometry was obtained with test No. 5 
Test section No. 5 was a distorted test section deformed by ap- 
plying internal pressure to it without installing the backing 
plates. It was recognized early in the program that fuel-element 
channels cannot be expected to remain flat in the reactor. Some 
bowing will happen as pressure differential and heating are ap- 
plied to them. A series of runs therefore was performed with test 
section No. 5 to determine the seriousness of small variations in 
geometry. 


section 


The exact changes in geometry were not determined 
from run to run. Once the test section was originally deformed, 
the backing plates were placed in position. However, the backing 
plates were flat and did not match the test-section curvature so 
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that further geometry distortion may have happened with 
changes in system pressure and temperature. 

Relative dimensions of the test section were obtained after the 
last run and are given in Fig. 8. They were taken 2, 12, 20, and 
32 in. from the inlet of the test section. An average cross section 
used to calculate the hydraulic diameter is also shown. It is im- 
portant to note that the cross section varies in the flow direction 
and that its average hydraulic diameter is 50 per cent larger than 
that of test sections Nos. 4 and 11. Typical heat-transfer data 
for test section No. 5 are given in Fig. 9. 
the average hydraulic diameter of Fig. 8. It is seen that the co- 
efficients are higher than those of test sections Nos. 4 and 11. 
The increased roundness in the geometry, therefore, brings the 
Further 
correlation of the data in Fig. 9 cannot be attempted owing to 
velocity changes and variation in cross section from run to run. 
The hydraulic diameter is not known for the tests plotted in Fig. 
9 and the spread between the data indicates that it must have 
varied with time. 


They are based upon 


heat-removal properties close to those in circular pipes. 


Data from test section No. 4 were 
used to determine the best type of correlation. 


Type of Correlating Equation. 
Plots correspond- 
ing to Equations (1), (2), and (3) were attempted. Fig. 10 uti- 
lizes the Sieder-Tate parameters, Fig. 11 uses film properties, 
while Fig. 12 is based upon bulk conditions. Heat balances are 
given in Table 1 to show the relative accuracy of each run. In 
three runs shown with check marks, the heat balances did not 





pH 
2 


agree within 10 per cent and the flow rates used in the Reynolds 
numbers were based on the electrical input rather than the orifice 
reading. 

Figs. 10 to 12 show that Sieder-Tate and Dittus-Boelter types 
of equation bring the runs closest together, and are to be preferred 
to a Colburn kind of equation. Relations identical to Equations 
(1) and (2), but with reduced coefficients probably could be used 
to correlate satisfactorily all the runs of test section No. 4. How- 
ever, upon closer examination it is noted that each run satisfies a 
different equation, and in particular the exponent of the Reynolds 
number varies considerably. The exponents are all greater than 
0.8 reaching the highest values when properties are evaluated at 
film conditions. Intermediate exponents are obtained for the 
Sieder-Tate parameters, and lowest values result when all proper- 
ties are taken at bulk conditions. For all three types of correla- 
tion, the exponent on the Reynolds number decreases with heat 
flux. For instance, from Fig. 10, runs 18, 19, and 20 at about 
1,000,000 Btu/hr sq ft, have an exponent of 1.1, while runs 8, 10, 
12, and 14 at 270,000 Btu/hr sq ft have an exponent of 2.2. As 
previously noted, the water-velocity distribution is expected to 
vary over the entire test section due to the constant change in 
water temperatures. It is anticipated that these circumstances, 
not previously treated analytically, could account for the varia- 
tion and increased values of the exponents. 

Entrance Effects. The thermal-entrance effects are illustrated in 
Fig. 13, which shows runs 11, 14, and 19 of test section No. 4 
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The entry length for these three runs is about 10 to 12 in. The 
corresponding entry length-to-diameter ratio is 55 to 65 and is 
considerably larger than the value of 15 reported by Hartnett for 
pipe flow [17] at this Reynolds number. It is important to note 
that Hartnett’s data apply to a fully established velocity distribu- 
tion, while in the present test section only 1 in. of calming hy- 
draulic section was provided. 
flow, however, 
the hydraulic 


Deissler’s analysis for turbulent 
indicates that even with simultaneous initiation of 
and thermal boundary layers, the ratio of thermal- 
entry length to diameter should not exceed 15. It is expected 
that Deissler’s results for pipe flow cannot be applied to rectangu- 
lar cross sections. For instance, in the case of triangular ducts, 
Eckert and Irvine [20] found from pressure-gradient measure- 
ments that the flow did not develop completely at a distance of 70 
diameters from the inlet. 
Even though the obtained entry lengths are considerably larger 
than those of reference [17], the effects of Reynolds number upon 
entry length appear preserved. Hartnett [17] obtained only a re- 
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Fig. 12 Correlation of heat-transfer data for rectangular test section 
No. 4 in terms of bulk properties 


duction of 30 per cent in entry length when decreasing his Rey- 
nolds number from 100,000 to 10,000. The present data shown 
in Figs. 13 and 6 indicate that the entry length varies but slightly 
with Reynolds number, decreasing slowly with Reynolds number. 
Values of 40 to 60 diameters are recommended in design for 
Reynolds number between 10,000 and 100,000. 

Flow Direction. Four test runs were carried out with upward 
As shown in Fig. 6, reversal of flow direction did 
not affect the heat-transfer coefficients. 


flow of water. 
This can be expected 
since at the high Reynolds numbers of the tests the natural-con- 
vection effects are negligible. 


Discussion of Results—Burnout Heat Transfer 


Insufficient data on peak heat fluxes were obtained to warrant 
an attempt at dimensionless correlations. Instead, data were 
plotted using groupings developed in references [11] and [19]. 
Fig. 14 shows a plot of the experimental results as compared to 
the equation proposed by McAdams, et al. [11]. Fig. 15 gives a 
similar comparison between the present data and Gunther’s [19]. 
Both figures reveal that the obtained burnout fluxes are ap- 
proximately one third those previously reported. However, as 
it was already pointed out, the data given in Table 2 give the 
heat-transfer flux needed to produce burnout at the corner of the 
test section. Since the heat-transfer coefficient at the shorter face 
is considerably lower than on the larger face of the section, lower 
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burnout fluxes may be expected. In other words, without heat 
generation in the corners, much higher heat fluxes would be at- 
tained before steam blanketing occurs. 

Examination of Fig. 14 shows that the use of an equation of the 
form proposed by McAdams gives an over-all correlation of the 
data within plus or minus 25 per cent. 
tion is 


The recommended equa- 


(2) = 0.29[400,000 plus 4800(T.a1 — 7,))V'/* (5) 
A BO 

Application of Equation (5) to a typical test reactor will now 
be considered. At the hottest point in the reactor the maximum 
heat flux is assumed to be 1,150,000 Btu/hr sq ft, corresponding 
to a heat generation in the fuel of 1.38 & 10° Btu/hr cu ft. If it 
is assumed that 7 per cent of this heat is generated in the water 
and side plates immediately adjacent to the fuel plates, it is found 
that the heat generation in the side plates is 2.44 « 10’ Btu/hr cu 
ft. If heat from the side plate is transferred in its entirety to the 
water flowing along the fuel plates, the maximum heat flux at the 
corner of the rectangular flow channel is then found to be 560,000 
Btu/hr sq ft. Fig. 14 gives a corresponding burnout flux of 
1,400,000 Btu/hr sq ft based upon water velocity of 35 fps and 
subcooling of 200 F. This means that the reactor power must 
increase to at least 240 per cent before any possibility of steam 
blanketing needs to be considered. Similarly, it is found that flow 
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Fig. 15 Burnout fluxes at various pressures and degrees of subcooling 
in rectangular test section 


reduction to 22 per cent of the original value decreases subcooling 
enough to produce burnout at the normal corner flux of 560,000 
Btu/hr sq ft 


Conclusions 


1 Reduced heat-transfer coefficients were obtained for turbu- 
lent water flow in thin rectangular channels measuring about 0.1 


<x 2.5 in. The heat-transfer rates are 30 to 45 per cent below 
those predicted from the Sieder-Tate equation and 15 to 30 per 
cent below most of the reported values for circular pipe flow. 

2 The use of bulk properties or Sieder-Tate parameters gives 
best correlation of the data from several runs. However, equa- 
tions with constant exponent on the Reynolds number are in- 
adequate. Test runs, when examined one at a time, have ex- 
ponents above 0.8, the value of the exponent increasing as the 
heat rate is reduced at constant water velocity. 

3 Heat-transfer coefficients decreased at the short side of the 
rectangular section. Values 15 to 30 per cent below those on the 
2.50-in. face were measured. 

4 Thermal entry lengths-to-hydraulic diameter ratios vary 
from 40 to 60 for Reynolds number between 10,000 and 100,000. 

5 Burnout heat fluxes were measured with heat generation 
Burn- 


out occurs at the corners and the values there are correlated by a 
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around the entire periphery of the rectangular channels. 


1959 


McAdams type equation. Values one third of those obtained for 
circular pipes were measured. 
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DISCUSSION 
S. J. Green* and B. W. LeTourneau® 


The authors are to be commended for undertaking a difficult 
experimental program in a relatively unexplored area. However, 
in view of their somewhat drastic conclusions concerning the 
relative inefficiency of rectangular channels as compared with 
circular tubes from the heat-transfer standpoint, a rather critical 
examination of the data would appear to be in order. 

The authors’ conclusion that the heat-transfer coefficients for 
15 to 30 
per cent) than those for circular tubes is in disagreement with all 
applicable data known to the writers. In particular, if only 
those data of Washington and Marks [16] with Vp. > 10,000 are 
considered, their values satisfy the Dittus-Boelter type of equa- 


rectangular-channel geometry are substantially lower 


tion [authors’ equation (1) ] using an average coefficient of 0.0239 
Also, it might be noted that additional rectangular-channel heat- 


transfer data, presented at the same meeting as the paper under 


discussion [21 |* agreed well with the Colburn equation [authors’ 
equation (2 

Also, 
lished flow of air 
through square, rectangular, and triangular ducts at Reynolds 
1000 to 330,000 and at temperatures from 80 to 


For rectangular channels, their data were correlated by 


») 


Livingood [22 


Weiland, 


and 


Lowdermilk, and have pub- 


heat-transfer triction-factor data for 
numbers from 
1320 F. 
an equation which was the same equation as that used for corre- 
lating their round-tube data 

It might have been desirable if the authors had conducted a 
few heat-transfer runs with round tubes to show agreement with 
published data 
mental validity of their data on rectangular channels 


This would have helped to confirm the experi- 


The authors estimate that their accumulated error in inside 


surface temperature 7',, could be as much as 6 per cent, but do 


not indicate the error in heat-transfer which would 
result from such an error in 7',. It is worth while to estimate 
this error for run 11, which is described in Fig. 5 and Table 1 of 
the paper Near the middle of the channel we can take, approxi- 
mately, 7) = 360 F, 7, = 130 F, and 7',/7, 1.7 
conditions, 7’, = 221 F + 6 per cent or 221 4 and 
T., — Tg = 91 + 13.2 deg F or 91 F + 14 Thus 
in this case a 6 per cent error in wall temperature would result in a 
14.5 per cent error in heat-transfer coefficient 


> 


flux (+3 


coefhicient 


For these 
13.2 deg F, 
5 per cent 


If the errors in 
per cent) and bulk temperature (not 
cluded, then the possible error 


stated) are in- 
in experimental heat-transfer 
coefficient will be near 20 per cent It is noted that the 6 per cent 


error in 7’, would be 8 per cent if the 1'/2 per cent error in cur- 
ised rather than '/2 per 


cent in the authors’ estimate of the wall temperature-drop error 


wu 


rent, given elsewhere in the paper, were 


Further, the possible error in the thermal conductivity of the 
than the 
11 this would cor- 


test section wall is believed to be +5 per cent rather 


| per cent stated in the paper [23 For ru 
respond to an additional possible error of 0.05(360 — 221) = 
7 deg F. This could result in an additional 8 per cent error in the 
heat-transfer coefficient 

In the case of the calculated heat-transfer 
at LOO psi 


compression, which is of such a nature as to de 


the 
due to the gasket 
rease the calculated 
coefficients, thus bringing the experiment il and calculated values 
into closet 


coefficients, 


largest error is that (6.5 per cent 


The combination of these errors is of 
if the consistent rather than 


random, the conclusion of heat-transfer coefficients 


agreement 


such a magnitude that error were 


15 to 30 pel 
cent lower than those for circular tubes could be entirely due to 
4 Bettis Plant (operated by the Westinghous« 
for the U. 8S. Atomic Energy 
\T-11-1-Gen-14), Pittsburgh, Pa. Mem. ASMI 
Bettis Plant, Pittsburgh, Pa. Assoc. Men 
®* Numbers in brackets designate Reference 


Electric Corporation 
AEC 


Commission under Contract 


ASMI 
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experimental 
geometry. 


error rather than to the rectangular channel 

It is noted that although the authors measured both the volt- 
age and current, only the current and published values of resistiv- 
ity were used to calculate the heat generation. Would the au- 
thors like to comment on why this technique was used and what 
was the agreement between the three possible methods (E/, 
I?R, and E?/R). 

The writers would be interested to know the results of any 


Here the 


inlet and exit water thermocouples gives a measure either of 


isothermal measurements any difference between 


errors in calibration or of estimated heat loss. Similarly, a com- 
parison between the water and outside-wall thermocouples de- 
Finally, 


urements of the exit-water thermocouples under net 


termines the accuracy of the wall thermocouples. mens- 
steam- 
generation conditions also would confirm the accuracy of these 
thermocouples. Were such measurements made? 

The authors state that the experimentally determined ex- 
ponent on Reynolds number, for a given run, is higher than 0.8 
and reaches values as high as 2.2 for some runs. This is very 
difficult to believe. If the velocity profile were to approach the 
extreme value associated with complete hydraulic roughness, 
then the Reynolds analogy could predict exponents as high as 
1.0, and such values might not be surprising at very high Reyn- 
olds numbers Values as high as 2.2, however, suggest possible 
errors due to experimental techniques, such an electrical thermo- 
couple pickup. The authors’ explanation that the continuously- 
varying bulk temperature and hence velocity profile can account 
for such high exponents is difficult to understand, especially 
since the highest exponents are associated with the small tempera- 
ture gradients, and the lowest exponents with the larger tem- 
perature gradients. 

The paper states that the deviation of the data below equation 


3) increases with increasing Reynolds number. It 


would be 
worth while to compare this trend with trends predicted by the 
available semitheoretical methods those 
Deissler [24] and by Martinelli 


that the Martinelli analysis predicts heat-transfer coefficients 


such as given by 


(25|. It might be noted here 
for rectangular channels which are essentially the same as those 
for round tubes, particularly for Np, > 0.7. 
The authors’ 


is not clear 


method of illustrating entrance effects, Fig. 13, 
It is usual to plot Aioe/ho versus L/D, where 
hy = local heat-transfer coefficient 

= fully established heat-transfer coefficient 
L= 


D hydraulic equivalent diametet 


distance from inlet 


Since, in this experiment, the fluid properties vary significantly 
with length, it might be well to correct ho to local fluid properties 
One thermal-entrance-effect 
information known to the writers is the data of Weiland and 


Sams [26] 


source of rectangular-channel, 
While in their experiment fully established flow was 
not reached, the data for both average and local heat-transfer 
coefficients in short rectangular channels were slightly higher 
than the values predicted from established round-tube correlations 
for the same L/D ratio 

With regard to the burnout data, it would be interesting to 
calculate the circumferential heat-flux distribution in the ree- 
tangular channels, assuming uniform volumetric heat generation 
and using conventional analog or digital techniques for solving 
two-dimensional heat-conduction problems. Nucleate boiling 
wall superheat correlations may be used, or the surface tempera- 
ture may be assumed constant. (To calculate the circumferential 
heat 


mates of the local heat-transfer coefficients relative to the average 


flux distribution under forced-convection conditions, esti- 


could be taken from the work of Eckert and Low [27], who indicate 
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that the local heat-transfer coefficient approaches the average 


value at a distance of 0.05 to 0.10 of the equivalent diameter from 
1 90-deg corner.) _It would not surprise the writers if conventional 
burnout heat-flux correlations would fit the present data satis- 
factorily if the maximum local heat flux were used 

It would seem overly conservative in most cases to apply these 
data directly to rectangular coolant channels in nuclear reactors. 
In conventional reactor fuel plates, the fueled region is of slightly 
less width than the coolant channels. While there is a significant 
umount of gamma heating in the extreme ends of the fuel plates 
ind in the side plates, this geometry and heating rate are nor- 
mally such that the heat flux through the narrow edges, and 
the portions of the wide sides nearest the corners, is much less 
than that opposite the fueled region 

For this reason, the rectangular channel test sections used for 

electrically heated burnout tests [28, 29] at the Bettis Plant 
have been designed with a reduced metal cross section in the 
corners, such that the heat flux on the narrow edges and on the 
part of the wide faces nearest the corner was about 20 per cent of 
that in the main part of the channel. With this type of channel 
and at pressures from 2000 to 600 psia) the burnout heat flux 
was found to be in fact higher than that for a round tube at the 
same conditions, and the burnouts did not occur preferentially 
it the corners 

Three minor errors which might be pointed out are that the 
ordinate scales of Figs. 7, 9, 12, and 13 are incorrect, the left- 
hand side of equation (4) should read 7’, T,, and the value of 
yy in the calculation of the adjusted D, should be 0.0094 rather 
than 0.094 
The 
ratio of heat-generation to heat-flux values indicates a heat- 
transfer area of 0.613 sq ft, while the test-section dimensions given 
indicate a heat-transfer area of only 2(0.118 + 2.045) (18)/(144 

0.541 sq ft for average heat flux, or an equivalent heat-transfer 
(test section volume thickness) of (0.028) [2(0.118 

1(0.028)] (18)/(144) (0.028) = 0.555 sq ft for local 
values of heat flux far from the corners, assuming uniform volu- 
metric heat generation. For runs 15B and 16B, the inlet tem- 
peratures necessary for consistency of the data are below freez- 
ing. For example, in run 16B the 196 F outlet bulk temperature 
at 138 psia corresponds to an exit enthalpy of 164 Btu per lb. 
The m-ss-flow rate may be calculated as W = pAV = [(7.8) 
(3600) (0.118) (2.045) (60.3)]/144 = 2840 lb per hr, so that the 
channel enthalpy rise is Q/W = 525,000/2840 = 185 Btu per 
lb. Thus the inlet enthalpy is calculated to be 164 185 = 

21 Btu per Ib! 

The writers conclude by pointing out that it is very difficult to 
obtain accurate heat-transfer coefficients for water at high Reyn- 
It is 
necessary to use the most refined experimental techniques, in- 
cluding the use of individually calibrated thermocouples for both 
surface and water-temperature measurements, extreme precau- 


Table 2 also contains some apparent numerical errors. 


irea 


2.045 


olds numbers, where the coefficients are very large. 


tion to eliminate heat losses, thermocouple electrical pickup, 
and thermocouple lead-wire conduction effects, minimization of 
wall-temperature-drop correction by the use of thin, high-thermal 
conductivity walls, and precautions to insure accurate knowledge 
of the channel dimensions during test. In the paper, the large 
cumulative effect of the many experimental errors on the accurate 
determination of foreed convection heat-transfer coefficients may 
be seen clearly 
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Surface Temperatures,”’ 


Authors’ Closure 


The authors wish to thank Messrs. 3. J 
Tourneau for their extensive comments. 


. Green and B. W. Le 
Because many of the 
questions they raise deal with the “lower’’ heat-transfer coeffi- 
cients reported in the paper, it is worth while to restate here that 
the present work deals with the simultaneous effects of (1) a rectan- 
gular cross section of high aspect ratio, (2) large variation of bulk 
properties in the flow direction, (3) large film temperature drops, 
and (4) undeveloped flow and thermal conditions at the inlet of 
the test section. Experimental and analytical studies are lacking 
in the case of a water-cooled channel for the specific test condi- 
tions reported herein. It is, therefore, important to avoid using 
preconceived relations or explaining the reduction in heat-transfer 
coefficients in terms of any one of the four test variables outlined 

Green and Le Tourneau mention the work of references [16, 21, 
and 22]. These tests do not exactly correspond to the conditions 
at hand. For instance, in references [21 and 22] the data are con- 
cerned with flow of air in rectangular channels of aspect ratio of 
only 5to1. Further, in reference [22] the flow is fully developed 
at the entrance of the test section and the surface conditions 
rough 

The discussers’ attempt to explain the low heat-transfer coeffi- 
cients in terms of experimental errors is not justified for several 


reasons: 


(1) Green and Le Tourneau assumed a 6 per cent error in the 
wall temperature. The 6 per cent mentioned in the paper is the 
sum of a 3 per cent error in surface temperature plus a 3 per cent 
error in heat flux measurements, as stated in the previous para- 
graph. Instead of the 13.2 deg F error calculated, one has a maxi- 
mum error of 5.5 deg F in film temperature drop. 

(2) Random errors such as reading errors should be neglected 
since they can lead only to data scatter. For instance, the 
1 '/, per cent error used for the current includes an assumed 1 per 
cent reading error 

(3) The errors stated for the thermal conductivity are maxi- 
mum errors according to reference [23] and Dickinson and Welch 
concluded that the accuracy of their equivalent Equation (4) was 

+3 per cent in spite of the quoted accuracy of the thermal con- 
ductivity. 

(4) The errors were maximized by the authors whenever they 
contributed to reducing the heat-transfer coefficient. The gasket 
compression calculation is typical of such instances. Further, 
errors contributing to an increase in the heat-transfer coefficient 
such as conduction along thermocouple wires and through the test 
section insulation should not be added when considering the maxi- 
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mum errors because they lead to reduced coefficients; in reality 
they should be subtracted. 

(5) The same experimental methods have been applied to 
circular pipes in another test loop and satisfactory comparison was 
obtained with the Colburn equation. 


The discussers note that the heat generation was calculated 
from values of current and local resistance. This technique was 
selected to account for local variation in heat flux along the test 
section. In every run the total heat generation calculated by this 
method was checked against the value obtained from voltage and 
current readings. The agreement was always excellent. It is 
also worth noting that the temperature instrumentation was 
checked out at ambient conditions before startup and that agree- 
ment was obtained within expected thermocouple calibration 
error within one degree. Similar checkouts were not performed at 
higher temperatures and the discussers’ suggestion is a worth- 
while one to consider in future experimental work. 

The variation in exponent of the Reynolds number with heat 
flux or ratio of surface to bulk temperature could be prematurely 
attributed to variation in the experimental errors with heat flux. 
Indeed, at the lower heat flux the film temperature drop is small 
and more susceptible to error but even at low heat flux the errors 
are not large enough to explain the high exponent on the Reynolds 
number. It is worth noting that in reference [30] the exponent on 
the Reynolds number varies and is greater than 0.8 ai the exit end 
of the test section. In Fig. 5(a) of that reference the exponent 
starts to change at Reynolds number of 50,000 and the effect of 
heat flux upon exponent is not clear, some runs even showing an 
increase in exponent at the lower Reynolds range with heat flux. 
Close examination of the same figure also reveals good agreement 
with the data shown in Fig. 12. For instance, at a bulk Reynolds 
number of 150,000 for a 2.1 value of 7,,/7T, (about an equivalent 
1.3 ratio in °R/°R) the ordinate in Fig. 12 is about 240 in agree- 
ment with the value given in reference [30] for a 1.4 temperature 
ratio (°R/°R). This agreement partially serves to establish the 
validity of the present tests and data. It also points to property 
variation as the potential reason for the reduction in the heat 
transfer coefficients. Unfortunately, the effects of the property 
variation could not be as firmly established as in reference [30]. 
The use of film properties did not bring all the data together as in 
reference [30] and the other variables used in the tests cannot be 
excluded as potential contributors to the reduction in heat-trans- 
fer coefficient. Theoretical solutions such as those given in refer- 
ences [24 and 25] are of no assistance in clarifying the picture. 
The solutions are based upon simplified models and upon pipe 
velocity distribution so that they can be expected to yield results in 
igreement with those for round tubes i 
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The discussers note that the authors report greater deviation of 
the data from Equation (3) at increasing Reynolds number. The 
authors’statement refers to Fig. 6 where the effects of all the test 
variables are compounded. On the other hand, examination of 
one run at a time in Fig. 12 shows smaller deviation of the data 
from Equation (3) at higher Reynolds numbers in agreement with 
the trend noticeable in reference [30]. Also, the entrance effects 
were not illustrated by means of usual plots because the authors 
do not know how to exactly calculate the fully established heat- 
transfer coefficient or to accurately account for fluid variation. 

With regard to the burnout data it is very hard to see how the 
local heat flux can be three times as large as the average heat flux. 
Conduction effects along a 28-mil wall could never produce such 
extreme peaking. It is also important to note that the proposed 
data apply to test reactors as noted in the first sentence of the 
introduction and not to power or propulsion reactors referred to 
by the discussers. 

The authors cannot disagree with the discussers’ comment that 
large heat-transfer coefficients are difficult to measure. They 
feel, however, that they have achieved data repeatability (see 
Fig. 12) and that the scatter is not necessarily tied to experimental 
errors and techniques but rather to the large number of variables 
considered. They would recommend separation of the variables 
and the use of another method besides electrical heating to study 
high heat-transfer coefficients. 

The authors thank the discussers for pointing out the errors in 
Figs. 7, 9, 12, and 13, Equation (4), and the value of Yr. The 
last two typographical errors were corrected in this last proof 
reading. 

However, Table 2 stands correct as is. Due to a typographical 
error the burnout test section dimensions appeared as 0.118 in. X 
2.045 in. The dimensions are 0.118 X 2.405 in. as now shown 
in the text (or as one could have deduced from Table 1, Fig. 8, and 
even Table 2). 


in Table 2. 


With these dimensions, the heat fluxes are correct 
In runs 15B and 16B, the use of the dimension 
2.045 in. and the lack of an allowance for an error in the heat 
balance at low velocities are responsible for the below-freezing 
The authors feel that 
they have considered enough variables without the added compli- 
cation of icy water. 


temperatures computed by the discussers. 
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Combined Geometric and Network 
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Director, Heat and Mass Flow 
Analyzer Laboratory, Columbia 
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Mem. ASME prises a 


Analog Computer for Transient Heat Flow 


A method to analyze three-dimensional transient heat flow is described, which com- 
continuous 


resistance medium but discrete capacitances. After de 


scribing the principles of the technique, the several components and their criteria are 
discussed, and the accuracy is illustrated by comparing the results of several simple prob 
lems with analytically determined values. 


1 Background and Purpose 


Sas term “‘analog computer’’ as used in this paper 


describes a device which permits the solution of a problem in on 
field of physics by interpretation of measurements carried out in 
other field 
solved by 


an ‘analogous’ For example, a problem in heat con- 


duction is appropriate interpretation of electrical 


measurements on an analogous system The basis lor such com- 
puters is the identity of mathematical expressions describing the 
problem and its “analog’’; if the same equation describes several 
fields (as the Fourier equation, which describes heat conduction 
diffusion, and electric currents in noninductive cables), a multiple 
ind the choice in which field to make measure- 


strongly 


inalogy exists, 


ments is influenced by the convenience with which 


measurements can be made 
The 


makes the electrical analog computer 


relative ease and accuracy of electrical measurements 


a versatile and powerful 
The 
While the paper deals 


with analogs for heat flow, the remarks have (in view of multiple 


tool.' One can distinguish between two types of analogs: 


geometric and the network computer 


analogies) broader applicability 


In geometric analogs the body in which heat flow occurs (hence- 
forth called “the object’’) is represented by a “model’’ which is 
geometrically similar to the object and related to the latter by a 
geometric ? suggested the electrolytic 
tank, metallized Beuken [3] de- 
scribed a method to represent finite boundary resistance values in 
Langmuir’s tank, while introduced a method to in 


( lude suC h resistances whe nh using metallized paper 


scale factor. Langmuir [1 


Bruckmeyer [2] used 


paper ; 
Kayan [4 


In network analogs the continuous system is broken down into 


discrete sections, with properties lumped at center or surface of 


Such analogs have found wide and deversified ap- 
But resistance networks with 
capacitances connected to the nodes and capable of handling 
Beuken [5 


which was later extended 


each section 


plications here only the use of 


transient problems are of interest first introduced 


this ton hnique 


ind used by the author 


16) 
Che geometric analog is limited to the solution of steady-state 
problems 


whereas the main limitation of the limped network 


It should be remembered here that the term analog is frequently 
used to describe differential analyzers, particularly electric differential 
analyzers. In this paper the term analog computer will be used only 
1s defined above and will always refer to electric analogs. 

? Numbers in brackets designate References at end of paper 

Contributed by the Heat Transfer Division and presented at the 
Annual Meeting, New York, N. Y., November 30—-December 5, 1958, 
of Tae AMERICAN Society OF MECHANICAL ENGINEERS, 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, June 13 
1958 Paper No. 58 4.69, 
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ies in its complexity if three-dimensional (or even complex two 
dimensional) shapes are to be handled. 

It is tempting to combine the two techniques, and this has been 
suggested by the author [7]. A brief report of an attempt to carry 
out such combined study has been given by A. V. Clark [8 
This paper presents results of work by this laboratory on such 

combined method,”’ carried out over many years and still con- 
tinuing 


2 Principle and Problems 


The Principle. Since no medium is known whicl 


would show appropriate electric capacity against a reference, a 


resistance 


method was sought in which the resistance medium is continuous 
but capacitances are introduced in lumped form. Also, if there is 
t boundary conductance between the surface of the body and the 
ambient, such conductance can best be represented by discrete 
external resistors. This method, though on the surface very 
simple, presents a number of difficulties both conceptual and ex- 
perimental. 

The two conceptual difficulties are both concerned with the 
.ransition from the continuous to the finite system: From ca 
pacitance to body resistance and from the latter to boundary re- 
sistance. The problem of transition from lumped capacitance to 
continuous resistance is explained, using a simple case (Fig. | 
for the solution of which no computer of any kind is required 
The problem is that of heating a slab, initially at zero tempera- 
ture, by imposing on one surface a step increase of temperature 
the other surface faces the reference ambient across a boundary 
conductance; the sides of the slab are perfectly insulated, so that 
heat flow is one-directional (perpendicular to the surface). In 
the mixed 


systems (continuous resistance-lumped capacitance 


handled by 


the problem can be the arrangement shown sche- 
matically in Fig. 1. 

\ tank T is filled with a resistance medium R in which connec- 
tors C are immersed, each connector leading to a capacitance C, 
(Cy, to Cy). The power supply is connected to electrodes E+ 
and Eb Initially the switch S is open and the system is not 
charged At time zero, S is closed and the condensers start to 
charge. But while the differential equation requires that at time 


T R 


Re 





system 
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0* capacitances C.-C, receive only an infinitesimally small 
current, in the system, Fig. 1, all capacitances receive finite 
values of current as soon as § is closed. 

The second difficulty related to the finite boundary resistance 
occurs both in steady and transient state, and is explained with 
regard to the ‘‘cold surface’’ of the model, Fig. 1. 

The methods [3, 4] are all similar in that the boundary re- 
sistance which by definition is uniform for the surface is broken 
down into a finite number of discrete resistances. The voltage 
distribution in the liquid represents the temperature distribution 
in the slab. In view of the one dimensional heat flow the voltage 
at any distance z from the (arbitrarily chosen) origin 0 should be 
the same for ail positions y. 

But the electrodes e;, e, . . . e, connecting the individual 
boundary resistance R; (film) to the continuous resistance R, 
which represents the body, result in a distortion of the voltage 
distribution particularly near e; . . . e,. In corners and in case 
of other two-dimensional problems this may cause serious errors. 

The Problems. Solution of a problem by means of the combined 
method, i.e., with distributed resistance and lumped capacitance, 
requires answers to a number of questions as follows: 

a. Resistance Material. The three main requirements for a 
suitable resistance material are: That it must easily accept a 
shape, that it must be usable with d-c, and that the resistivity 
must be in an acceptable range. 

The last requirement rules out all commonly used metals. 
Mixtures of mercury and glass beads proved unmanageable. 
Finally, a solution of tetra-butyl-ammonium picrate (hereafter 
abbreviated thap) in Dowtherm was tried, based on suggestions 
of C. A. Kraus and R. M. Fuoss [9], and found acceptable. 

The electric resistivity was measured in a glass cell with pal- 
ladium electrodes 114/32 in. apart; the electrodes were 1'*/s9 in. 
wide and immersed 1'/, in. deep using 200 v d-c. 

The resistivity changes also characteristically with tempera- 
ture and therefore the solution must be used at constant tem- 
perature. Fig. 2 shows the resistivity versus temperature for 
several values of concentration. If plotted on semilog paper (con- 
centration on logarithmic abscissa) the resistivity decreases at 
a temperature of 22 C from 140 megohm-in. (at 27 mg/1000 cc) 
linearly to 20 megohm-in. at 270 mg/100 cc; then the curve 
flattens out to reach a value of 3.9 megohm-in. at 618.4 mg/100 
ee. Atsuch high concentration solution of the liquid takes a long 
time and, although saturation may be reached only somewhat 
higher, it is not practical to use concentrations in excess of about 
600-620 mg/100 cc. Obviously, being a liquid the mixture 
satisfies the requirement of easy forming, and the resistivity range 
is appropriate. The third requirement, usability with direct 
current, requires explanation. 

The solution is slightly electrolytic; under influence of d-c a 
slow dissociation takes place. This makes itself felt in form of a 
deviation from Ohm’s law: A plot of current measurements made 
in the afore-mentioned cell at different voltages does not go 
through zero; instead, a finite voltage is needed to achieve any 
current flow. This limiting voltage is sometimes called the ‘“‘back- 
emf’’ and depends to a considerable extent on the electrode ma- 
terial. 

A number of electrode materials have been investigated, and 
the following values have been found for the back-emf, when used 
in thap: 


nickel 

platinum 

palladium 

platinized palladium 
palladinized palladium 
platinized platinum 


L471 v 
Lislv 
1.027 v 
0.546 v 
0.509 v 
0.444 v 


The difference between the performance of the last two metals 
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Fig. 2 Resistivity-temperature relationship for a solution of thap in 
Dowtherm 


is not big enough to warrant the much higher cost of the platinized 
platinum; palladinized palladium is suitable. 

b. The Electrodes. The choice of the electrode material can- 
not be separated from that of the resistance material, Thin 
sheets and strips can be used for the surface electrodes, which can 
be bent so as to fit any desired surface contour. The “electrodes’’ 
connecting to the capacitances present a special problem, and are 
therefore given a special name (capacitance connectors) and are 
dealt with later. Suffice it here to say that for the sake of economy 
in the triai runs nickel electrodes and steel capacitance connectors 
were used; in spite of this conscious choice of less than the best 
materials results were good (see Section 4). 

c. The Tank. Any material that does not react with the re- 
sistance liquid and which has sufficiently high electric resistivity 
is acceptable. Because of the ease of working to any shape, wood 
was used and made impermeable on the inside by an appropriate 
coating, Bakelite BL3128 of the Bakelite Corp., and Weldwood 
Waterproof Resorcinal glue made by the U. 
were found acceptable; 


8. Plywood Corp., 
since the former requires baking at 
elevated temperatures the experiments described in Section 4 
were carried out in tanks coated on the inside with several layers 


of Weldwood. 


Capacitance Connectors. In steady-state (electrolytic) tanks, 
wires can be used as probes to determine the potential distribu- 
tion because no current flows through the probe. But in tran- 
sient experiments current flows through the capacitance connec- 
tors, and thus the resistance of the connector becomes significant; 
the resistance is that between the liquid and the connector. 
If a wire is used as connector, the resistance from the latter 
to the bath approaches infinity as the diameter approaches zero. 
To understand this, consider a hollow cylinder of liquid re- 
sistance material; the outside diameter can be at any reasonable 
distance from the axis of the cylinder, while the inner surface is 
given by the diameter of the wire. The resistance of this liquid 
cylinder is proportional to the logarithm of the diameter ratio 
which approaches infinity as the wire diameter approaches zero 
Conversely an increase of the diameter of the capacitance con- 
nectors shorts out increasingly large areas of the resistance 
medium. The diameter of a cylinder (in case of two-dimensional 
work) or of a sphere (for three-dimensional work) has been 
selected to date on an empirical basis only; the experiments are 
reported under Section 4. 

Boundary Electrodes and Boundary Resistance. Accepting the 
concept of discrete resistors to be connected by individual elec- 
trodes, the question arises as to how many electrodes should be 
used; obviously, theoretically, the greater the number, the better 
the results. But as the number increases, the required value of 
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Fig. 3 Electrodes and separation 





the individual discrete resistance increases, thus putting a prac- 
tical limit on the number of electrodes. The question thus be- 
comes, how large an error is introduced by a finite number of 
electrodes? Fig. 3 is an enlarged view of a straight surface, a few 
electrodes, and the surrounding resistance medium. Contrary to 
theoretical requirements each electrode imposes over its length 
the condition of an isothermal face. Thus even for a given num- 
ber of electrodes, the size of /., of the electrode, and thus the un- 
If the 


distance between two electrodes, 14, becomes small, the resistance 


covered space /4 between them remains an open question. 


of the liquid connecting the two electrodes also becomes very 
Elec- 


trodes can be placed nearer if between the electrodes an isolating 


small; thus through the liquid a near-short is created. 


fin is placed which protrudes into the liquid. Then the size, li, 
of the insulator is open 

In connection with this electrode-boundary resistance problem, 
the optimum solution can be sought from the viewpoint of a cor- 
rect value for the rate of heat flow (represented by the current) or 
the temperature distribution either at the surface itself or in the 
body near the surface (represented by voltages). These problems 
are mentioned here only for the sake of completeness. Work on 
these aspects is in progress and not part of this paper. 


3 Test Body 


In order to examine the usefulness of the several techniques 
selected, it was decided to solve several problems of two and three- 
dimensional nature, for which analytical solutions were available 


or were developed These problems are listed hereafter. 


a) A 
mensional problem is initially at 


square bar of sufficient length to be treated as two-di- 


zero temperature. At time 
zero the four surfaces are subjected to a step increase in tem- 
perature. What is the temperature rise in the body? 

b) The same bar is heated from two adjacent sides, while the 
two other sides were kept at zero temperature 


f \ 


ing the six surfaces to a step increase of temperature. 


cube, initially at zero temperature, is heated by subject- 


d) The same cube is heated by providing for a step increase 
on tour laces, plac ed sy mmetrically, while the other two faces are 
kept at zero temperature 

The two-dimensional problems were analyzed in a wooden tank 
with an inside opening of 8 X 8 in., with 100 capacitance connec- 
tors spaced at 0.8-in. center distance. Steel cylinders of '/s, '/«, 
ind */,-in. diameter were used as condenser connectors 

The three-dimensional tests were carried out in a cubic tank of 
S-1n 


inside length. Steel spheres of */,-in. diameter were pro- 


vided with */\-in. holes along one diameter; the spheres were 
strung on glass rods and held at proper distance from each other 
by glass tubes slipped over the rods. The rods in turn were held in 
place by a board provided with regularly drilled holes at the bot- 
Teflon in- 


sulated wire was soldered to each sphere, and all wires of one 


tom of the tank and by a spacer plate at the top. 


vertical row of spheres were brought to a plug-in connector. A 
total of 125 spheres provided 5 nodes in each direction. Because 


of symmetry the cubic tank was used to represent '/, of the 
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actual cube. Thus, for problem c, three faces joining in one corner 
were equipped with electrodes and subjected to the step increase 
while the insulated. In 
problem d, investigating again one eighth of the cube, two “heat- 


in voltage, other three faces were 
ing faces’’ are now adjacent to the cold face; thus at some edges 
An edge separating a heated from a cold 


face belongs to both and should at the same time be connected to 


discontinuities occur. 
full voltage and to ground. To overcome this difficulty in the 
tank, such electrodes were cut back at the edges, so that the 
do not touch. 


“zero electrode’ and the ‘“‘hot electrodes’’ 


4 Results 


Two-Dimensional Tests. In order to separate the errors in- 
curred by lumping and by the bypassing of the capacitance con- 
nectors by the liquid resistance medium computing experiments 
were compared with the analytical solutions, and in part with the 
solution obtained on a lumped resistance-capacitance network. 
Fig. 4 indicates the position of the points in question; Figs. 5 and 
6 show the temperature (or voltage) rise as function of Nyro for 
problem a. Results with */sin. capacitance connectors check 
3 and 


Also for point E the 


very closely with analytical results, while curves based on ! 
'/,-in. connectors are much less accurate 
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Fig. 4 Position of measured points (square bar) 
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Fig. 5 Temperature versus time; points B and E (problem a) 


Transactions of the ASME 

















T 








+ 
i 


a 
cee Es 





z= 
+ 


$—$—$4 44+ ___4._+__4}- ++ 















































+—_+—_-+ 


























colculoted 


3/6" electrodes 


“4 electrodes 


° 1/6" electrodes 

















Neo , FOURIER NUMBER 1.0 


time; points A and C (problem a) 
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Fig. 7 Temperature versus time; points E and G (problem b) 


check between the lumped network and the analytical values are 
very good, particularly for 10 nodes. Fig 
problem }, points E and G. 


7 shows the results for 
Errors are appreciably smaller than 
in problem a; but results appeared encouraging enough as to 


warrant proceeding to three-dimensional problems 
Three-Dimensional Tests 

Heated From All Sides The potential 

rise of points on the space diagonal was computed analytically 

and compared with values found from the model tests. The 

points on the space diagonal are numbered 5-5-5 (next to the hot 

corner) to 1-1-1 (next to the center of the cube) 


a. Cube Proble me ¢ 


Results are 
shown in terms of fractions of the steady-state temperature of 
each point. Results at point 5-5-5 are poor, improving rapidly 
as one progresses toward the inside (see Fig. 8 
values of a 5 node R-C network are included. The error in case 
of the tank method far exceeds that of the lumped network. In 
Fig. 9 two further pairs of curves compare the analytical with tank 
values, the latter taken by means of a probe at a point A very 
near the ‘‘hot corner,’’ and at a point B very near the center of 
the cube; The 
A large. 
But it can be seen that the major influence is in the position rather 


For comparison 


neither point had a capacitance connector. 
errors at point B are very small, in case of point 


than in the presence or absence of a capacitance connector. 
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Fig. 8 Temperature versus time; node 5-5-5 (problem 
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Fig. 9 Temperature versus time; points A and B (problem d) 


Obviously, the proximity of the electrodes is detrimental to 
If the potential at a given point, e.g., at node 5-5-5, 
should be determined, one should see to it that this point is not 


at the first node. 


accuracy. 


But an increase of nodes is expensive; e.g., if 
point 5-5-5 should be the third instead of the first node from the 
corner, the total number of nodes (provided uniform lumping is 
used) would be 154 = 1335! One way to provide for uneven 
lumping is to provide two tanks made to uneven scale: a smaller 
tank representing the hot corner, including the node 5-5-5, and 
the other representing the balance of the cube. 

The relationship of the two tanks is shown in Fig. 10; the insert 
tank had a side of five inches and was equipped with nine electrode 
connectors. In the large tank a solid cube of 1.6-in. side length 
was put in the hot corner to displace the liquid resistor. The 
corresponding faces a, b, e, f; a, b, ¢, d; and a, d, h, e, on both 
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Fig. 10 Main tank and insert tank 
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Fig. 11 Transfer electrodes on insert tank 
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Fig. 12 Comparison of results for node 5-5-5 (problem c) 


tanks were covered with metal squares as shown in Fig 11 for the 
insert tank Corresponding plate s were connected with wires 

In the insert tank nodes along the diagonal were characterized 
by double numbers 11-11-11, 22-22-22, and 33-33-33, the latter 
being in the new hot corner, and the first near the point a, in Fig. 
10 Point 22-22-22 corresponds exactly to point 5-5-5 in the 8-in. 
tanks 


in. tank (repeated from Fig. 8) with those found in the insert tank 


Fig. 12 compares the results for this point found in the 8- 
and with analytically determined values. The insert tank gives 
satisfactory results 

Fig. 13 shows that results for the two other nodes are less satis- 
factory. The large error is easily understandable for point 33-33- 
33. With regard to the innermost point 11-11-11, it is surprising 
that the direction of the error is reversed as compared with node 
1-4-4 in the 8-in. tank. In the latter (Fig. 8) the potential value 
in the model was at any time No below the analytical value, while 
for node 22-22-22 the analytical value was below the model value. 
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Fig. 13 Temperature time curves for points 11-11-11 and 33-33-33 
(problem c) 
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Fig. 14 Temperature versus time; points 1-1-1 (problem d) 


\ possible explanation may be seen in the use of the plates, as 
shown in Fig. 10, which provide a relatively short path from the 
power electrode on faces ee, fz, ge, he, and be, ¢s, f2,' ge to the capaci- 
tance connector. 

b. Cube Heated From Four Sides—Problem d. Figs. 14 and 15 
show the results for the several nodes, again in comparison both 
Point 
5-5-5, now receiving heat only from two instead of from three 
sides, heats up slower in comparison with conditions of problem c 


with analytical results and with lumped. network results. 


(compare Fig. 8). In light of the more favorable results of this 


problem, no tests were run with the small insert tank. 


Conclusion and Further Work 


The problem mentioned with regard to the finite boundary re- 
sistance and the connectors is obviously also relevant for the 
and the insert tank. A 
study along the lines indicated in Section 2 is under way. Further 
problems with cubes, possibly in connection with node and elec- 
trode size should be studied, including arrangement and size of 
the insert tank. 


connections between the main tank 


Also the replacement of cylinders and spheres 
by square bars and cubes (for symmetry in square sections) 
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Fig. 15 Temperature versus time; points 2-2-2; 3-3-3; 4-4-4; and 5-5-5 


might be considered. The concept of separate tanks connected 
across sets of corresponding electrodes can be extended to make 
possible handling of three-dimensional problems in bodies con- 


sisting of materials with different conductivities 
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DISCUSSION 
R. P. Benedict® 


In this paper, the author introduces an interesting technique 
for approximating solutions to transient heat-transfer problems. 
The author has indicated not only the various advantages in- 
herent in the subject geometric-network analog, but he points out 
at least three basic difficulties. These are: 

1 Back-emf phenomena which opposes impressed emf. 

2 Necessary transition from lumped-capacitors to distributed 
resistances. 

3 Unavoidable currents through capacitance connectors which 
complicate connector designs. 


Since the writer has worked with the Liebmann type of re- 
sistance network analog‘ and with the electrolytic tank analog,5 
it might be appropriate to make a few remarks along these lines. 

While electrolytic tank analogs offer the distinct advantage 
of distributed resistances and facility in representing complex 
geometric boundaries, their use is nevertheless attended by many 
complications. Foremost is the chemical reaction which always 
takes place at the electrode-electrolyte interface as an outcome 
of the transition from an electronic to an ionic mode of electrical 
conduction. As mentioned by the author, this polarization effect 
means that locally the electrolyte is not an ohmic conductor, a 
fact quite important when it is remembered that all the electric 
field equations on which the various analogies depend are based 
on Ohm’s law. It is evident that the polarization effect must be 
minimized. This may be accomplished by proper choice of elec- 
trode, electrolyte, and power supply. The following suggestions 
are offered in this connection: 

a Inexpensive electrodes of practically zero transition re- 
sistance may be obtained by coating electrodes of any material 
with a colloidal dispersion of graphite such as Aquadag. 

b Tap water has been found to be quite satisfactory as an 
electrolyte. 

ce To avoid the back-emf phenomena associated with d-c 
While 


the use of a-c is accompanied by its own peculiar polarization 


tanks, a-c power supplies usually have been employed. 


phenomena (i.e., the appearance of an unwanted capacitive re- 
actance) this effect may be reduced by going to the higher fre- 
quencies (500 to 1000 eps). 

It is realized that the R-C networks, as used in this paper 
with the geometric analog, require a d-c power supply. However, 
if a Liebmann resistance network were used in place of the R-C 
network in the combined geometric-network analog, certain ad- 
vantages would accrue; for in bypassing polarization effects, we 
avoid also the transient measuring problems associated with 
capacitors charging and discharging. This latter advantage is 
realized because in the Liebmann method, time is handled as a 
stepwise variable so that at any instant we have a steady-state 
problem 

It is true that by 
electrolytic tank we do not avoid the use of lumped circuit ele- 


combining a Liebmann network with an 
ments, nor do we resolve the problem of unavoidable currents 
through the connectors from these lumped circuit elements to the 


electrolyte. In fact, in view of these difficulties and considering 


3 Development Engineer, Westinghouse Electric Corp., Lester, Pa. 

‘R. P. Benedict, ““Liebmann Network Approximations to One- 
Dimensional Transient Heat Conduction Problems,’’ ASME Paper 
No. 58-——A-221. 

*R. P. Benedict and C. A. Meyer, ‘‘Electrolytic Tank Analog for 
Studying Fluid Flow Fields Within Turbomachinery,” ASME Paper 
No. 57—A-120. 
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VARIOUS SOLUTIONS TO TWO DIMENSIONAL TRANSIENT 
HEAT TRANSFER IN A HOLLOW CIRCULAR CYLINDER 





a + —- 
ANALYTICAL SOLUTIONS 
ANALOG SOLUTIONS 





DIMENSIONLESS TIME (+ 
Fig. 16 


also such items as uniform resistivity, ph, leakage, and evaporation 
of the electrolyte, as well as model building, leveling, and size of 
connectors, it is the feeling of the writer that at present the in- 
creased difficulties encountered probably more than offset the 
gains to be realized by the use of the proposed geometric-net work 
analog 

As a matter of interest, a problem similar to probiem a of 
the paper was solved by a two-dimensional Liebmann network 
Programming consists simply of evaluating resistance ratios at 
each of the nodes of the problem. Results are presented in Fig. 16 
and these are seen to compare favorably with the numerical solu- 
tions (programmed for high-speed digital computers) of Dr. Y. 8 
Tang of the Westinghouse Electric Corporation 


H. P. Fullerton’ 


The development of techniques for applying the electrolytic 


*Switchgear and Control 
Philadelphia, Pa 


Division, General Electric Company, 


150 


MAY 1959 


tank to three-dimensional transient problems is a valuable con- 
tribution in a difficult field. The author and his associates are to 
be commended for their ingenuity in meeting the problems which 
they faced and for the guidance this paper extends to others who 
may attempt work in this area. 

One wonders, however, whether in using the insert tank in 
problem c, care was taken to choose the resistivity for the insert 
tank so as to retain the proper scale relation. No mention is 
made of this in the paper and it appears that this may account for 
the anomaly at point 11-11-11 mentioned in the text and shown 
in Fig. 13. If electrolyte of the same resistivity as that in the 
main tank were used in the insert tank, the effect would be the 
same as though a fluid of lower resistivity had been placed in that 
portion of the main tank. 
tential field would result 


A corresponding distortion of the po- 


Author's Closure 


Mr. Benedict suggests the use of ‘“‘steady-state’’ 
ated 


tanks, oper- 


with a-c, using the Liebmann techniques. The writer's 
enumeration of difficulties encountered in such procedure seem 
sufficient to make the R-C method of the paper the more at- 
tractive. This leaves as alternative for the handling of three- 
dimensional problems a straight Liebmann technique, using a 
taking the 


place of capacitances, or the combined technique suggested in 


three-dimensional resistor grid with “drain resistors”’ 
this paper. If the three-dimensional shapes are simple, e.g., 
cubes, the Liebmann method may become attractive, even if the 
number of resistors required would be large But for irregular 
shapes, or for regular shapes with irregular boundary resistances, 
the number of resistors required for the Liebmann technique 
becomes prohibitive and the method proposed in the paper seems 
to be clearly in the advantage; the good results for point E 
Fig. 7) indicate that it will not be necessary to place too many 
capacitors, and the use of insert tanks (Fig. 10 may also be used 
at points of irregularities of shape 

The reason for choosing simple geometries and two-dimensional 
cases as examples in the paper is that only for such cases are 
Therefore the ex- 
ample in Mr. Benedict’s Fig. 16 is not helpful for a comparison 


analytic solutions available for comparison 


even of two-dimensional cases between the combined techniques 
and the Liebmann method. The latter is in distinct disadvantage 
for complex two-dimensional shapes 

Mr. Fullerton is correct in saying that the different size of the 
insert requires a different time constant of the system in the insert 
However, this time constant can be achieved by changing resis- 
tivity or capacity (capacitance per unit volume). Experi- 
mentally it was easier to use the same resistance liquid as in the 
main tank, changing the capacity The capacitance of each 
node in the insert tank was reduced to '/, of the value of that in 
the large tank, corresponding to the square of the length ratio 
In this manner the Fourier numbers (dimensionless times) in 
main and insert tanks remained identical 

As mentioned in the paper, the anomaly at point 11-11-11 is 
probably due to the proximity of the two electrodes with dif- 


ferent potentials. 
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Heat Transfer From an Annular 
Fin of Constant Thickness 


The two-dimensional solution for annular fins, given in graphical form with the ef- 
ficiency plotted versus the design-parameter group, (b — a)[(2h)/(kw)]'/*, supplies 
needed design information for fins of small height-to-thickness ratio and large height- 
to-inner radius ratio. 


The one-dimensional solutions previously given for annular 


fins are accurate for height-to-width ratios of the order of 10 or more, while the two- 
dimensional results for rectangular fins are useful as approximations to annular fins 


when curvature 1s not large. 


With height-to-width ratios less than 10 and for annular 


fins with large curvature, design of fins can be computed with the results presented in 


this paper. 


A PROBLEM ot practical engineering interest is that of 
the diffusion of heat through an annular fin of constant thickness. 
Interest in fins is prevalent in the many fields of engineering, such 
as air conditioning, heat-exchanger apparatus, electric motors, 
and electronics, where, using extended surfaces, engineers have 
been able to build more compact and efficient heat-transfer ap- 
paratus. 

The effect of adding fins to a surface immersed in a fluid de- 
pends upon a number of factors of which the most important is the 
increase of exposed surface area. For heat transfer from an in- 
terior isothermal surface to an external fluid, there are two fac- 
tors which determine the rate of heat exchange: 


1 Resistance to heat flow through the convective film of the 
ambient fluid. 

2 Conductive resistance to heat flow within the metal fin it- 
self 


Numerous solutions to the problem of diffusion of heat in an 
annular fin of constant thickness have been presented, but all 


cases refer to one-dimensional heat flow [1—5].2 An assumption 


1 This paper is based on an MS thesis submitte 


by H. H. Keller to 


| 
the Mechanical Engineering Department of the University of Pitts- 
burgh. 


? Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division and presented at the 
Annual Meeting, New York, N. Y., November 30—-December 5, 1958, 
of THE AMERICAN Society OF MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME He adq larters, August 5. 
1958. Paper No. 58—-A-116 
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inherent in these solutions is that the fin thickness is sufficiently 
small compared to its height that temperature gradients normal 
to the length direction of the fin may be neglected, For fins of 
large height and small thickness this assumption produces negligi- 
ble error. For fins with a small height-to-thickness ratio, the 
error introduced by the one-dimensional assumption becomes im- 
portant. 

Avrami and Little [6] have applied a two-dimensional analysis 
to a flat rectangular fin. Their results should approximate closely 
an annular fin of small height-to-thickness ratio and large inner 
diameter; effectively, a fin where the curve is unimportant. The 
work reported here supplies the solution for annular fins of small 
height-to-thickness ratio where curvature is important; Le., 


where the fin height is comparable to the inner diameter 


Treatment of the Problem 


The diffusion of heat in an annular fin of constant thickness re- 
duces to the problem of obtaining the steady-state temperature 
distribution in a finite, hollow, circular cylinder representing the 
annular fin. Several restrictions and assumptions are inherent 
in solving this fin problem: 


1 The flow of heat and temperature distribution throughout 
the fin are independent of time. 


, 


2 The fin material is isotropic. 
3 There is no heat source in the fin. 
1 The flow of heat is both radial and axial; i.e., two-dimen- 


sional. 


These establish the Laplace differential equation as describing 





surface of end of fin = 
surface of side of fin 

inner radius of fin 

outer radius of fin 


heat-transfer coefficient an te 


modified Bessel function of first (h 


kind of order n and of argu- 
ment r 

Bessel 
second kind of order n and of 


modified function of 


argument z 
thermal conductivity of fin ma- 


heat-flow rate 


dimensionless heat-flow rate 


ratio of radius of any point in 
height of fin, r = 
a) 


- a 
function of (a. wu, 2B, ); 
hb a 


equation (12 


temperature 


fin thickness 


roots of characteristic equation 


radius of any point in fin 4) 


ratio of fin height to thickness 

axial position of any point in fin 

dimensionless axial position at 
any point in fin, z = 2’/w 

ratio of heat-transfer coefficient 
times fin height to thermal 
conductivity of fin material, 
pw = h(b — a)/k 


terial 
an integer 
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dimensionless temperature dis- 
tribution in fin 


temperature at base of fin 


fin efficiency 
function defined in equation 
(8) 





the heat flow within the fin, which in dimensionless form is given 
by 
1 OV OV 
+ ~+ § 


= 0 (1 
or? r Or oz? ) 


The restrictions on the boundary conditions for the fin are: 

1 The heat-flux rate from the fin surface at any point is 
directly proportional to the temperature difference between the 
surface at that point and the ambient fluid; i.e., Newton’s law 
of ‘cooling applies with a constant convection coefficient. 

2 The ambient temperature of the fluid is constant at an ar- 
bitrary temperature of zero. 

3 The base of the fin is isothermal at temperature V,’ 

4 The fin is of constant thickness. 

These restrictions lead to the following dimensionless-boundary- 
condition equations: 

At the base of the fin 


At the end of the 


At the sides of the fin 
oV(r, 0 


or 


oVir, 1 
or 


and 


The solution to equations (1) through (5) yields the dimen- 


sionless temperature distribution in the finite hollow cylinder 
representing an annular fin [7] 


where the la, } are the positive roots of 
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The heat rate at the base of the fin can be found by integrating 
the Fourier law for conducted heat-flux rate across the base of the 


fin 
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Fig. 1 Annular fin of constant thickness 
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The fin efficiency, used as the basis for determining a fin design 
is defined as the ratio of actual heat-flow rate to the heat-flow rate 
that would exist if the entire fin were at the temperature of the 


tka 


.) 


base of the fin. This definition in equation form is 


from which 
Q’ 4rak VBS 
The fin efficiency then becomes 
28S 2Q 
6) jb+a b) 
+" ; 
' a aB\ 


jo+a 


Y la , aB{ 
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Numerical Reduction of Solution 


The fin efficiency g was evaluated with the geometrical and 
operating variables given, i.e., b, a, h, and k given; this required 
the determination of g per equation (13). ‘The dimensionless 


heat-flow rate Q was found in equation (10) to be 


a > 


n=1 


t= (3) ef -m(5)y 2” 


\a,,? + p? + 2uB} yp (; : 2 n) 


where 


ot a, J) a, K.( a,b 
b—a/ \b-—a b-a 
c a,b \} r a,a oe ( a,b 
k(t + Ks (;2") . 7 aes 
r im 1 ( a,b t 
b a b-—a { 


b-—a 


(15) 


and 


a,b 
—a 
a,b 
p= Is ( (16) 
b a b a/\ 


The Bessel functions contained in Q are evaluated by interpola- 
tion of tabulated values when the value of the argument is less 
than 20. For values of the argument greater than 20, an asymp- 
totic approximation of the Bessel functions [8] For 


} Was used. 
z> 20 


K,(z) ¢ ( = 


The roots of the characteristic equation, equation (8), were 
computed roughly in graphical form and the values were then re- 
refined by Newtonian interpolation. A graphical plot of the right- 
hand and left-hand sides of equation (7) 
uw = 0.2and B = 10. 

A simplification of the solution may be made by noting that the 
expression 


a,a j a,b E a,b \I . 
—}<a.K - — pKy - _— K, 
a,J; (; = -) ) 1 (; is .) Le (; e -) a, ( 


is shown in Fig. 2 for 


aa j . § a,b . ( a,b \| ‘ ( 
0 = 7 K = Ko 4 - 
. (; _ .) ‘ (5 ~ .) - b a ) f b a 
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Fig. 2 Graphic solution of (a’/8) = (2ua’)/|(a’)? — u*] for u = 0.2, 
8=10 


may be approximated closely by 


a,a 
a, Ky 
b-—a 


a,a 
Ko ( — ) 
b-—a 


and may be approximated by 


lb-a 


? 
2 a 


a, (a, > 20) 


The results are plotted in Figs. 3, 4, 5, 6, 7 for five values of 


bh/a. The selection of a/(b — a) determines b/(b 


- a) and b/a. 
The values chosen for a/(b — a) are 


a 
= 0.2, 0.5, 1.00, 5.00, and 10.00 
b—a 


These values cover a range from tall fins on small-diameter tubes 
to short fins on large-diameter tubes. 

The values chosen for 8 were picked to cover a practical range 
from a short wide fin to a long and narrow fin. 


for B are 


The values chosen 


B = 1/3, 1.00, 3.00, 10.00, and 50.00. 


The values chosen for 4 were such as to keep the magnitude of 


J a,,a a,b \| 
) \ Onl (; ~ -) + blo (; — .)t 


( a,a ) is il ( a,b \| 
b—-a 7 on i - OF) 
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Fig. 3 Efficiency of annular fin of constant thickness 
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Fig. 4 Efficiency of | 
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Fig. 5 Efficiency of annular fin of constant thickness 


(2u8)'/* (where [2u8)]'/* = b — a[2h/kw]'/*) within the limits of 


0.05 to 10. For (2u48)'/? less than 0.05 the efficiencies are almost 
unity and for (2u8)'/* greater than 10 the efficiencies are 
almost zero 


Discussion 


The results of this paper have been compared to the results ob- 
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tained by Avrami and Little [6] for the two-dimensional rectan- 
gular bar and to the results obtained by Gardner [5] for a one- 
dimensional annular fin. 

Comparison with the results obtained by Avrami and Little in- 
dicate that theirs could be used in designing an annular fin of small 
height-to-width ratio and small b/a ratio; as the b/a ratio in- 
creases, the results of Avrami and Little deviate from those ob- 
tained here. The values are compared in Table 1. 
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Fig. 6 Efficiency of lar fin of constant thickness 
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Fig. 7 Efficiency of annular fin of constant thickness 


Table 1 Comparison of two-dimensional rectangular and radial cases his results are satisfactory for large height-to-width ratios 


(b — a) ¢( two-dimensional) Shown in Table 2 is a comparison of Gardner’s one-dimensional 


2h \'/2 -, b-—a values with the two-dimensional values of this paper. As can be 
(= ) a.” ” . seen, Gardner’s results apply with decreasing accuracy as the 
00 0.50 0.925 O15 732 length-to-width ratio is reduced below 10. 

00 2.00 0.480 176 151 

00 5.00 0.196 182 038 


pa 


10 0 50 0 915 O06 6R5 Table 2 Comparison of one and two-dimensional radial cases 


10 2 00 0 168 445 150 — ¢ (two-dimensional 
10 5.00 0.190 178 038 b b-—a 
“ 

20 0.50 0.912 902 665 ‘ (one- 
20 2 00 0.465 138 140 *  dimen- = 50 
20 5.00 0.190 178 038 cw sional ) 

0.920 915 906 
00 0.50 0 895 880 600 0.470 168 145 
00 2.00 0.390 0.370 126 0.195 190 178 
00 5.00 0.150 0.138 032 

0.920 912 902 

5.00 0.50 0.875 0.855 560 - : 0.465 165 138 

00 2.00 0.340 0.318 110 0.190 182 172 


+00 5.00 0.118 0.110 027 
0.900 895 SSO) 


*bh/a 1.00 infers rectangular results of Avrami and Little;  ‘ ; 0.390 390 370 
compare efficiencies given for rectangular resultsat b/a = 1.00with =‘ 0.150 150 138 
efficiencies given for radial results at other values of b/a, at con- 
stant b a/w and constant (b — a)(2h/kw)'/:. 3 0.875 875 855 

K 7 0.340 0.340 318 

( ‘omparison with the results obtained by Gardner indicate that 3. . 0.120 0.118 0.110 
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A Variational Method for Fully Developed 
Laminar Heat Transfer in Ducts 


A variational method ts presented for determining fully developed velocity and tempera- 


ture distributions for laminar flow in noncircular ducts. The heat addition to the fluid 
ts taken to be untform in the axial direction, but a variety of thermal boundary con 

ditions are considered around the periphery of the duct cross section. Several illustrative 
examples are given, and comparisons are made which show good agreement with available 
exact solutions. These examples include ducts of rectangular and circular-sector cross 


sections. 


Introduction 


Tes UTILITY of variational methods as a tool for 
the solution of complex problems in the theory of elasticity and 
other fields of applied science is well established. The power of 
the variational technique has yet to be brought to bear on prob- 
lems in heat transfer, except for recent work on unsteady con- 
duction, e.g., Biot [1].! 

In this report we are concerned with forced convection heat 
transfer, and our aim is to formulate and apply a variational pro- 
cedure for solving the problem of fully developed heat transfer in 
noncircular ducts. To demonstrate the method, velocity and heat- 
transfer solutions will be found for laminar flow in passages of 


rectangular and of circular sector cross section. In all cases, the 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division and presented at the 
Annual Meeting, New York, N. Y., November 30—-December 5, 1958, 
of THe American Socrety oF MECHANICAL ENGINEERS 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, August 4, 
1958 Paper No. 58—A-92. 


Nomenclature 


heat input per unit length will be taken as constant along the 
passage, but a variety of thermal conditions around the boundary 
of the cross section will be considered. 

Attention is focused here on noncircular ducts because of their 
increasing application in the design of high performance heat ex- 
change equipment. In particular, such noncircular configurations 
may be found in the cores of nuclear reactors, where fuel elements 
may be arranged in various patterns. 

A brief description of the general variational procedure will be 
given in the next section, and this will be followed by the specific 
problems. More detailed information on the mathematical 
fundamentals may be found in a number of texts, for example, 
Hildebrand [2]. 


The Variational Procedure 


General Considerations. The basis of the variational method to 
be applied here may be stated as follows: Suppose that the physi- 
eal principles involved in a problem are expressed mathematically 
as a differential equation, e.g., a conservation law. Then, the 
solution of the problem is given by the function which satisfies 





cross-sectional area of duct 
half-length of the short side of temperature; 
rectangular duct, or half-length 


of the side of square duct 


perature; 


half-length of the long side of ocity 
rectangular duct 

1 dp variable 

a constant, . 

friction factor 

specific heat at constant pressure 

hydraulic diameter, 44/wetted 
perimeter 


radial co-ordinate; ro, duct radius 


w! 


fluid velocity; 
representation 


gw dz X dimensionless co-ordinate, z/a 

co-ordinate along the short side of 
rectangular duct, or along the 7-7, 
side of square duct Q/8k 

dimensionless co-ordinate; y/b for 


dimensionless angular co-ordinate, 
7,,, bulk fluid tem- 6/6 

wali temperature angular co-ordinate in circular-sec- 
a, mean fluid vel- tor duct; 65, half of opening 
angle, radians 

of a dependent absolute viscosity 
fluid density 


dimensionless temperature, 
®,, bulk fluid tem- 


perature 


integrand of the variational inte- 
gral 

fully developed heat-transfer 
efficient 

variational integral 

thermal conductivity 

Nusselt number, hD,/k 

static pressure 

heat-transfer rate to fluid per unit 
length of duct 

dimensionless radial co-ordinate, 
r/o 
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rectangular duct or y/a for 
square duct 

co-ordinate along the long side of 
rectangular duct, or along the 
side of square duct 

axial co-ordinate along duct length 

dimensionless temperature in cir- 


w 
cular-sector duct, : 


Q/k6o 
bulk fluid temperature 
aspect ratio of rectangular duct, 


b ‘a 


dimensionless temperature, 
. 
Q/A(y + Dk 


perature; $,, average wall tem- 


: d,, bulk fluid tem- 


perature 


partial differentiation 
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both the differential equation and the given boundary conditions 
Now, according to the calculus of variations, there corresponds 
to this differential equation an integral with a very particular 
Namely, 


mum when it is evaluated using that function which satisfies the 


characteristic that the value of the integral is a mini- 


differential equation and boundary conditions. As discussed later 
these boundary conditions must be of certain permissible types. 
When other functions satisfying the boundary conditions but not 
the differential equation are introduced, the value of the varia- 
tional integral exceeds the minimum. These properties suggest a 
procedure for obtaining approximate solutions of differential 
equations by utilizing the variational integral. The simplest 
scheme would be to select a number of trial functions which satisfy 
the boundary conditions, and then introduce them successively 
into the variational integral. That function leading to the least 
value of the integral would then be regarded as most closely ap- 
proximating the correct solution. More sophisticated schemes, 


e.g., the Ritz method, provide systematic ways of finding trial 
functions with good minimum properties. 

We now proceed to express mathematically the discussion of 
the preceding paragraphs. The initial item to be considered is 
the correspondence between a differential equation and its varia- 
tional integral. For convenience, rectangular co-ordinates will be 
ised 

Consider a function V which may depend upon z and y; for 


example, Vo may represent the velocity or the temperature. 
Suppose that the physical problem has been mathematically 
formulated, so that we know the differential equation which V 
must satisfy. Then, according to the calculus of variations, we 
must first find (by inspection, experience, or other means) a func- 


tion Fiz, y, | | } . which when inserted into the relation? 


) ( or) o ( *- ) oF 
dz \OV, oy \OV*, oV 
When F is known we can 


ields the differential equation for V 
write the variational integral /, 


] = Sf, F(z, y, V, V,, V,)dz dy 


The variational method states that the solution of the differential 
equation for V also gives a minimum value for J. 
(s a simple illustration, suppose that V is governed by 


\’ , 0 


re 


From experience, we propose that the F function be 


h 


Inserting this expression into equation (1) and differentiating 
The 


variational integral corresponding to equation (3a) is then written 


vields equation (3a so our choice for F is a proper one 


in accordance with equation (2) as 


I . i + on 
. V,)?]dx dy 


Leaving this illustration and returning to the general situation, 
it is important to note that the correspondence between varia- 
tional integrals and differential equations is contingent upon par- 
ticular boundary conditions being satisfied. Specifically, V must 
be either constant on the boundary or else the normal derivative 
of } We will later show that these 
restrictions provide no obstacle for the types of problems con- 
sidered in this report 


must be zero on the boundary. 


? This is often called the Euler Equation 
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Once a variational integral has been found, it may then form 
the basis of an approximation method for finding V. The proce- 
dure is to choose trial functions for V which make the integral as 
small as possible (the exact solution gives the absolute mini- 
mum). A well-known technique for minimizing a variational 
integral is called the Ritz method. According to this approach, 
V is written as a linear combination of functions, each of which 
satisfies the boundary conditions, i.e., 


oVo(z, y) + Ai Vilz, y) 


AsV: Zz, y 


With this choice of V, the function F may be evaluated and the 


integration in equation (2) carried out. The result is an expres- 


sion containing the Ap, . To achieve a minimum value 
of the integral J, this expression is differentiated separately with 
This 


provides (r + 1) simultaneous equations which can be solved for 


respect to each A, and each result is set equal to zero 


the unknown constants Ao, ... A,. In practice, this procedure 
may be applied successively using an increasing number of terms 
in the series (4) until convergence to a desired accuracy is 


achieved. The convergence is hastened by choosing the form of 
the V, in accordance with any intuition or knowledge that one 
may have about the physical form of the solution. 

Application of these ideas to forced convection flow and heat 
transfer will be made in later sections. 

The Conservation Equations. Heat transfer involving fluid mo- 
tion is governed by three conservation laws: Mass, momentum, 
and energy. For the situation of fully developed heat transfer in 
passages, the velocity is taken to be unchanging along the length 
Under these conditions, the conservation of mass equation is 
identically satisfied, and conservation of momentum reduces to a 


single equation whic h, in rectangular co-ordinates, Is 
1 dp 
oy? a dz 
This equation expresses the fact that the pressure forces must ex- 


The 


viscous flow boundary conditions require that the velocity vanish 


actly balance the viscous forces in a nonaccelerating flow 


on all solid surfaces. 
Now we turn to the energy equation for fully developed flow, 
which, when written in rectangular co-ordinates, is 


oT ‘T’ O77 
peu 1 — - 
<< oy* 
Fluid property variations and viscous dissipation have been 
uniform heat input per unit 
length, the condition of fully developed heat transfer is that the 
oT’ /dz, is 


neglected. For the situation of 


axial temperature gradient, constant. An over-all 


energy balance gives 


oT Q 


oz pAic, 


where A is the cross-sectional area of the passage and Q is the heat 
input per unit length. 
gradient and 0*7'/dz? 


Using this expression for the temperature 
0, equation (6) becomes 


Q v7” OPT 
-afk — 
A oz? oy? 
The thermal conditions around the boundary of a given cross sec- 
tion may take on a variety of forms, e.g., uniform temperature, 
uniform heat flux, partially insulated, ete. Specific choices will be 
made in the examples given later. 
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For cylindrical co-ordinates, the conservation laws analogous 


to equations (5) and (7) are 


O*u 1 Ou 2) dp 
r or 2 dz 


( or7 oT ) " 
=k (7a) 
or? og? 


Variational Principles for the Velocity and Temperature Distributions. 
We are now in a position to write the variational integrals cor- 
responding to the conservation laws of the previous paragraphs. 
We deal first 
According to 
the rules for the variational procedure, the initial task is to find a 


Consider a rectangular duct as shown in Fig. | 
with the differential equation (5) for the velocity 
function F which, when inserted into equation (1), yields equation 
is the 


5). It is easily verified that the following choice for F 


proper one 


The variational integral corresponding to equation (5) is then 


eta +b 
i= (u,? + yw 2? + QWCujdz dy 
—a b , 


The physical requirement that the velocity be zero on the 
boundaries satisfies the mathematical restrictions on the bound- 
ary conditions needed to assure correspondence between the dif- 
ferential equation and the variational integral 

Next, we go to the energy equation for the rectangular duct. 
By inspection, the function F which leads to equation (7) is found 
to he 


and the corresponding variational integral is 


~+a +b one tod (@ +. 
] = io t+ 2, + 24 - -~T \drdy 9 
b kA a 

e a -o 

It has already been noted that in order for the variational integral 
9) to correspond to the differential equation (7), the permissible 
boundary conditions are that the temperature is constant or that 
the normal derivative is zero. Since the zero derivative case is 


trivial (all walls insulated), equation (9), as it stands, is only 


, 
usable for ducts with uniform peripheral wall temperature. To 
treat a situation where the peripheral heat flux is specified, i.e., 
will be made in 


prescribed normal derivative, a modification 


the formulation. This will be demonstrated later when illustra- 
tive examples are given. 

Finally, we give the form of the variational integrals in cylindri- 
cal co-ordinates. The specific geometry to be considered is a 


) 


passage whose cross sect 1On, 48 show nin Fig «18a cireular sector. 


Fig. 1 Co-ordinate system for rectangular duct 
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Fig. 2 Co-ordinate system for circular-sector duct 


The variational integral for the velocity, corresponding to equa- 
tion (5a), is 


eT t 
I | f u,? + 
0 — 


Further, the variational integral corresponding to the energy 


Q u 
T+ 2 
Ss (*) li 


The comments with respect to boundary conditions which have 


+ 2 «) rdrd@ (8a) 


equation (7a) is 


r| rdrd@ (9a) 


been made in connection with equations (8) and (9) also apply to 


equations (8a) and (9a 


Applications of the Variational Procedure 


Velocity Distribution in a Square Duct. It is seen from the energy 


equation (7) or its variational integral (9) that a knowledge of the 
velocity distribution is required before the temperature problem 
can be solved. Hence we shall first direct our attention to the 
velocity distribution. 

In solving the velocity problem, we are able to benefit from 
knowledge accumulated in the solution of the torsion problem in 
elasticity theory. There is a direct correspondence between equa- 
tion (5) and the differential equation for the stress function in the 
Further, equation (8 
‘energy integral’’ for the stress function. 


torsion problem. 


is analogous to the 
Thus, from the results 
given by Timoshenko (reference [3], p. 283).for torsion of a 
square bar, we can write a first approximation to the laminar 
velocity distribution in a square duct as follows 


1 


(x? — a*)(y? — a®) [0.5844 a® + 0.1185(2? + y?)] 
0.2808 a® 


(10) 


To test the accuracy of this expression, a comparison with the 
exact velocity solution [4] is made in Fig. 3. Profiles are plotted 
along a co-ordinate axis and along a diagonal. The agreement is 
generally very good, except for deviations near the corners. 
With the thought of improving the approximation, the varia- 
tional method was reapplied using an expression of slightly greater 
complexity than equation (10). The steps in the calculation will 
be shown in some detail, so that future calculations which proceed 
in the same manner may be described only in outline form. The 
profile chosen for the improved approximation is 
*y(y2? — a*)[Ag 4 (11) 


A,(x? + y?) + Asr*y?] 


The factors (x? — a?) and (y? — a*) assure that each of the three 
terms satisfies the condition of zero velocity at the boundaries 
ta,y = +a. 


r= ! The constants Ao, A, and A; remain to be de- 
termined. 
Substitution of equation (11) into the variational integral (8) 


leads to the result 
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EXACT SOLUTION(4) 


VARIATIONAL METHOD 
1" APPROXIMATION 


2 APPROXIMATION 


5 
y/o 


(a) Profile along axis of symmetry ( =0,0< y < 1) 
a 


EXACT SOLUTION (4) 


VARIATIONAL METHOD 
i" APPROXIMATION 
2™¢ APPROXIMATION 


4 5 
a/o of y/o 


(b) Profile’along diagonal (x = y) 


Fig. 3 Laminar velocity distribution in a square duct 


0.71111 0.29799 A,*a* + 0.010642 A,*a"” 


A,?a* + 


0.48762 AgA,a‘ 0.040635 AyA,a* + 0.077400 A,A.a* 


. Oc / ‘ 
eCl\, atest ie 


To find the values of Ao, A,, 
we differentiate equation (12 


and A, which will make / a minimum, 
with respect to each A and set the 


resulting equations equal to zero. This operation yields 


»/ 


A 


> 1.4222 Aga* + 0.48762 


1,‘ 


0.040635 Aza*t 


0 — 0.48762 Aya*? + 0.59598 A,at 


16 
+ 0.077400 A.a* 4 
45 


0 — 0.040635 Aya? 0.077400 A,a*t 


+ 0.021284 A,a* + 


Simultaneous solution of this group of linear algebraic equations 


provides the following values 


0.082030 C 0.24610 ¢ 


‘ a‘ 


0.58984 CC’ 
’ 


a’ a 


rhe mean velocity is computed from its definition, giving 


*a 
A,a® + 
5° 


0.28111 Ca‘ 


After substituting the A values into equation (11) and dividing 


by the mean velocity, we find 


a*) (0.58984 a* 


+ 0.082030 a*(x? + y*) + 0.24610 x4y?] (lle 


Velocity profiles evaluated from equation (lla) are compared 


with the exact solution in Fig. 3. The agreement is very good, 


160 / may 1959 


and it is seen that the exact solution is approached more and more 
closely as more terms are added to the approximating poly- 
nomial for u. 

Now that accurate velocity distributions have been computed 
for the square duct, we may proceed to solve the heat-transfer 
problem. 

The Square Duct With Uniform Peripheral Temperature. The first heat- 
transfer situation to be considered is a square duct whose walls 
ave a uniform temperature around each cross section, while the 
heat transfer per unit length in the axial direction is constant. 
Under these conditions it is clear that the wall temperature will 
increase from section to section along the length of the passage. 

For this situation, we may use the variational integral (9 
without essential modification. For convenience, we introduce 
dimensionless variables 


'-T 
Q/8k 


w X 


, 


ind rewrite equation (9) as 


s+] +1 u 
in f ?,? 4 ?,? 1 4 
J -1 —1 = 


where the condition at the duct boundaries is that ® = 0 


») dX d) 9a 


When 
the temperature around the boundary is uniform, the temperature 
distribution across the section would be expected to have a shape 
not too different from the velocity. So, as a first approximation 
for ®, we select a simple polynomial having a form similar to equa- 
tion (10), i.e 


’ 


@ = (X2 — 1)(¥?2 — 1) [eo + «(X? + Y? 

Then, this expression for ® is introduced into the variational in- 

After carrying 

there 
un- 


tegral (9a) along with u/d from equation (10). 
out the integration and minimizing (by 


results two linear algebraic equations for computing the 


differentiation), 
known ¢ and c;. When these constants have been determined, 
the expression for ® becomes 

md = (X? 13a) 


j 


1)(¥? — 1)[—0.9276 + 0.1368(X? + Y? 


Once the temperature distribution is known, we may proceed 


to find the Nusselt number. Using the definitions® 


AD, n Q 


Nu h= — — — (14) 
Perimeter X (7, — 7, 


k’ 


+ The perimeter appearing the definition of h is the heated perimeter. 
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and noting that for a square duct: D, = 


2a, perimeter = 8a; 
the Nusselt number expression becomes 


a Me 
Ah -7T) & 


yu 


The dimensionless bulk temperature ®, may be computed from 
the usual definition 


1 1 
>, = fi f ®(u/a)dX dY = —0.5555 
0 0 


Then, the final result for the square duct with uniform peripheral 
wall temperature is 


Nu = 3.60 (15) 


An exact calculation by Marco and Han [5] gives a value in pre- 
cise agreement (to three significant figures) with equation (15). 
Hence there is no need to use more complex approximations for 
® or for u/%. 

The Square Duct With Uniform Peripheral Heat Flux. The next situa- 
tion to be considered is the square duct with a uniform heat flux 
around the walls of each cross section, as well as uniform heat 
input per unit length in the axial direction. The condition of uni- 
form heat flux is equivalent to the requirement that the normal 
derivative of the temperature be constant at all the duct walls. 
This type of boundary condition is not consistent with the varia- 
tional integral as given in equation (9). However, there is no 
difficulty in deriving an alternate variational integral which will 
apply for the situation under consideration. 

To accomplish this, we write the temperature 7 as the following 
sum 


Q 


T = 7,+ — 
16ka? 


(16) 


Now, the condition of uniform heat flux requires that 
(oT On)waii = Y Ska 


where n is direction of the (outward) normal to the duct walls. 
It then follows by differentiation of equation (16) that 07',/dn is 
identically zero on all the boundaries. The differential equation 
for 7 is found by substitution in equation (7), giving 


o°7; r aT, Q ( 1) a7) 
Oz? oy? ~ fat \ a , 


The variational integral corresponding to this differential equa- 
tion is found by inspection to be 
- 1) 6] dX dY 


+1 +1 ’ 
I= [ ce (diy)? + ( 
J-1 J-!1 = 
(18) 


where @; = 7,/(Q/8k) and X and Y are dimensionless co-ordi- 
nates. That the differential equation (17) actually corresponds 
to the integral (18) is assured because 7; possesses proper bound- 
ary conditions. 

To use the variational integral (18) we must select a trial func- 
tion which has zero normal derivatives at the boundaries. 
example, we might choose 


For 


dg: = (X® — 1)"%¥2 — 1)"[dp + d,(X? + ¥%)] (19) 


This, however, would be an unfortunate choice, since not only do 
the normal derivatives vanish at X = +1, but also 
¢; itself vanishes. Hence equation (19) does not provide any 
freedom in the variation of the temperature around the bounding 
walls and, in fact, this variation is already specified by the last 
term on the right side of equation (16). 


+1,Y = 
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As a more suitable first approximation, we select the following 
expression for ¢;: 


di = eo[(X? — 1)? + (¥? — 1)?] + e,(X? — 1) ¥2 — 1)? (20) 


Since the simple velocity expression of equation (10) led to favora- 
ble Nusselt number results in the last example, we shall use it 
again, at least for a first approximation. Inserting equations 
(20) and (10) into the variational integral (18) leads to 
I =— (. 8 + = + = af + 2.7743 eo + 1.6741 
= — —@¢ —= ¢@ 2. i e 
at eo eee a si 
Differentiating with respect to eg and e, to determine the values 
which minimize J gives the two equations 
al 32 nae 
a 0 — 12e, + — + 2.7743 = 0 
% 5 
32 512 


Oo— €o 


1.6741 = 0 
5 105 * F 16 


which are satisfied by e = —0.16032, e, = —0.13290. Using 
these constants in conjunction with equations (16) and (20) pro- 
vides the temperature distribution as follows: 

T a* + ¥* 


= > + 


% = O/8k 2 


ra 1)?] 
3 a= \2 


— 0.13290(X? — 1)*%(Y¥? — 1)? (21) 
The Nusselt number and heat-transfer coefficient are again de- 
fmed as in equation (14), except that 7’, must now be regarded 
as the average value of the wall temperature which varies around 
the periphery. Introducing the dimensionless temperature @ and 
using a bar to denote an average, the Nusselt number becomes 


9 


Nu = 


Pw — M 
The average wall temperature and the bulk temperature which 


appear in equation (22) are found by integrating the temperature 
distribution (21) as follows: 


1 1 l 
$= | f (>) x-i dY + [ (ra | = 0.5814 (23) 
” /J0 /J0 
l 1 u 
>, = f f o( ~ } dX dY = —0.06803 (24) 
0 Jo . 


With these values, the Nusselt number result is 


Nu = 3.09 (25) 


To check the accuracy of this result, we proceed to reapply the 
variational method with more complex expressions for the ve- 
locity and temperature distributions. For the velocity, equation 
(11a) will be used, while the temperature distribution is written 
as follows: 

xX? +. ¥? | 
—s + Pi 

x* + ¥* . , 

= ——— + fa [(X* — 1)? + (¥# — 1)?] 


+ fx(X? — 1)2(¥? — 1)? + fal(X* — 1)* + (¥? — 1)4] 


Evaluating the variational integral (18) and following the same 
minimizing procedure as before, it is found that 
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fy = —0.17532, f; = —0.14003, fs = —0.020761 


The values of d,, and @, needed in computing the Nusselt number 
may be evaluated from the defining equations (23) and (24), 
giving 


d, = 0.5826 o, = 0.0654 


Finally, substitution into equation (22) leads to 


Nu = 


3.09 26 


to that ob- 
tained by use of the simpler velocity and temperature profiles [see 


This result is identical (to three significant figures 


equation (25 Hence there is no need to continue the problem 
by using still more complex profiles. 

The square duct with uniform peripheral heat flux was studied 
by Cheng [6] using separation of variables, and this led to rather 
involved results in the form of infinite series. He computed a 
Nusselt number of 3.82, but this value appears inconsistent with 
his results for rectangular ducts of higher aspect ratios and hence 
may be in error. A computation by Irvine [7] using a finite dif- 
ference procedure yielded a Nusselt number of 2.9, a deviation of 
6 per cent from the current results 

Flat Rectangular Duct With Uniform Peripheral Heat Flux. Next, we 
turn to flow in a flat rectangular duct and, in particular, to a duct 
b/a = 10. The situation of uniform 


peripheral heat flux will be considered first. 


with aspect ratio ¥ = 
An approximation for the velocity distribution in a flat rectan- 
from the solution of the cor- 


[3], p. 285), i.e 


gular duct can be taken directly 
responding torsion problem (ref. 


” 


1S » 
Y 5 -Y) 
3 hae ee 
y? 


a 2 l 


— 0.632/y 
where X 


It is noted that this form is not symmetric in Y, and as written 


r/aand Y = y/b 
only applies when Y is positive. A comparison with the tabula- 
tion of the exact solution given in reference [4] shows that equa- 


tion (27 


is a good representation of the velocity, and its use poses 
no special problems if care is exercised with regard to signs. 
To formulate the problem for a rectangular duct with uniform 
peripheral heat flux, it is only necessary to retrace the steps al- 
ready taken for the square duct 
tion (16), 


In a manner analogous to equa- 
we write the temperature 7’ as 


= (4° + 9F* 


ya 


(28) 


For Y = l It is 
easily verified that 07',/On must vanish on all the boundaries of 


square duct), equations (28) and (16) coincide. 


the rectangle in order that the condition of uniform peripheral 


oT | ¢ 
heat flux be satisfied (ie, = . 
On | wall 4(a + b)k 


the variational integral corresponding to equation (18 


+] e+! 1 
l= if | (ou 1 ar 
0 J 0 Y 


Further, 


becomes 


i) 6, | dX dY (29) 


We are now in a position to use this 


j 


( 
where @ = 7’; : 
$(1 + 


integral to find the temperature distribution. 


y)k 


As a first approximation for @,, we select a function of the same 


form as has already been used for the square, i.e., 


gd; = go(X? 1)? 1)? (30) 
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This expression is introduced in integral (29) along with the 
velocity (27). After minimizing the integral, the constants in the 
resulting simultaneous equations were determined for an aspect 
ratio y = 10. This yielded the following temperature distribu- 


tion 


X? + 10Y? 


~ Q/A + dR 


- 0.11703(X? 


+ 


1)? — 3.0660(¥? — 1 


0.038328(X? 1)2(¥? — 1)? (31) 


The Nusselt number is again defined by equation (14), where 


7, is interpreted as the average wall temperature. For a rectan- 


“ 


gular duct heated on all boundaries 


fa ( . 
¥Y + 


In terms of @, the Nusselt number becomes 


D, = ), heated perimeter = 4(a + 
l 


4 Y 
Ge ~GQit 


The average wall temperature and the bulk fluid temperature 


may be computed from 


s I 
?, = 
a-+ b J0 


a 
— e f Waals 
0 
l l 
7 f (d)xe:d¥ + ( ) ya: dX 33 


) « 
lel, 
f f ax a¥ 34) 
0 o © 


The numerical values of these quantities corresponding to equa- 
1 I 1 


tion (31) for @ and equation (27) for u/i are 


d, = 0.9472, %, = —0.2537 
The corresponding Nusselt number is 3.03. 

In addition to the Nusselt number, which represents an averag- 
ing of the heat-transfer performance around the periphery, it is of 
interest to know the wall temperature variation around the 
boundary. This has been evaluated from equation (31) and is 
plotted in Fig. 4. The hot spot in the duct corner is strikingly 
displayed there. 

Higher order approximations for the temperature profile were 
then employed to investigate the accuracy of the predictions of 
Nusselt number and wall temperature variation. Since the tem- 
perature variation in the y-direction is much larger than that in the 
z-direction, the original expression (30) was modified by adding 
For the 
second approximation, a term g:(Y? — 1)? was added to equation 
(30). After evaluation of the constants go, . 
of variational method, the expression for the temperature dis- 


polynomials in y to provide more degrees of freedom. 
.- gs, by application 


tribution became 


T X?+ 10¥? 
~ Q/41 + yk 2 
— 0.11703(X? — 1)? — 6.9554(Y? — 1)? 


— 0.038328(X? — 1)2(¥? — 1)? — 3.8894(Y? — 1)8 


o 


or 
(oo 


(34) leads to 
This yields a 


Computation of d, and @, from equations (33) and 
values of 0.67779 and —0.58018, respectively. 
Nusselt number of 2.89, which is 4.5 per cent lower than that 
found for the first approximation. 

This change in the Nusselt number suggests that still another 
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Fig. 4 Wall temperature distribution for fully developed heat transfer 
in a rectangular duct of aspect ratio 10 


term be added to the approximation for the temperature distribu- 
tion. After affixing the term g,(¥* — 1)‘ and following through 
the variational procedure to determine the constants, the re- 
sulting temperature distribution is 


- X? + 10Y? 


— 0.11703(X? 


— 0.038328(X? — 1) ¥? | 


12.05(Y? — 1 1.9697(12 — 1)4 (36 


This gives an average wall temperature, @, = 0.6143, and a bulk 
temperature @, = —0.6379; with a corresponding value of 2.90 
for the Nusselt number. The very close agreement of this result 
with that of the previous approximation indicates that further 
It is 
worth while noting that this Nusselt number agrees closely with 


refinement of the temperature distribution is not necessary. 


the value 2.88 given by Cheng’s [6] exact solution 

The wall temperature variations corresponding to the succes- 
sive variational calculations are shown in Fig. 4. The figure il- 
lustrates the manner in which the sequence of solutions con- 
verges to the final result. 

As has already been pointed out in reference [6], the Nusselt 
number for rectangular ducts with uniform peripheral heat flux 
remains at a value of about 3 as the aspect ratio increases from 1 
to 10. This is in contrast to the situation of uniform peripheral 
wall temperature (with constant heat input per unit length), 
where the Nusselt number increases from 3.60 to 6.79 as the as- 
pect ratio changes from 1 to 10, and approaches 140/17 (8.23) as 
the aspect ratio approaches infinity. The Nusselt number for 
flow between infinite parallel plates with uniform axial heat flux 
is also 140/17. So, the parallel plate channel represents the 
limiting situation for the rectangular duct with uniform periph- 
eral wall temperature. An important practical consequence is 
that the parallel plate model cannot be used to predict hot spot 
corner) temperatures when the heat flux is uniform around the 
periphery. 

The Fiat Duct With Partially Insulated Walls. The severe corner hot 
spots which are indicated in Fig. 4 are in some part due to the 
stipulation of uniform heat flux on the short end walls. 


applications, e.g., is likely 


In some 


reactor core assemblies,‘ it that the 


* An often-used reactor core assembly is composed of parallel fuel 
plates, producing essentially uniform heat flux, supported at their 
edges by unheated structural materials 
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short end walls would essentially be insulated. We now proceed 
to investigate, using variational techniques, the situation where 
there is uniform heat flux along the long sides of the rectangle, 
with the short sides being insulated. The aspect ratio y will again 
be taken as 10. The analysis follows closely that which has been 
already presented for the flat duct with uniform heat flux all 
around the periphery. 
To reduce the problem to one which is suitable for the varia- 
tional method, we write the temperature as 
T,+— x (37) 
8yk 
where X represents x/a as before. In order to satisfy the condi- 
tions that 
oT Q oT 
= = at z= +a, —— 


+ =( at 
ox 4bk oy 


y= +b 


it is necessary that 07',;/On be zero on all the boundaries. 
convenience, we may rephrase equation (37) as 


For 


vy + 1 X? 
o = o: + y 2 


/ ¢ 
where @ = rls 1 j and gq = P| rT id rs 
(l + Wr r YA 


easily verified that the variational integral for @, is again given 


It is 


by equation (29). Inasmuch as ¢; must satisfy the same equa- 
tions and boundary conditions as in the preceding illustrative 
problem, we may make direct use of the previous results. Then, 
the result for the temperature distribution corresponding to y = 10 
may immediately be written as 


11 X? 
1Q 2 


— 0.11703 (X? 


— 10.291(Y? —1)? — 0.038328(X? — 1): 


— 12.05(Y2 — 1)* — 4.9697(Y? — 1)* (39) 


The average wall temperature and fluid bulk temperature can 
be computed from equations (33) and (34), giving: d, 
—1.3038 and @ = —2.1095. The Nusselt number is then 
evaluated from the basic definition (14). After noting that the 
heated perimeter is 4b, the Nusselt number expression becomes 

. 4 
Nu = = 


4.96 
du = Pr 


This value for the Nusselt number lies between those for uniform 
peripheral heat flux and uniform peripheral wall temperature. 

The wall temperature variation as computed from equation (39) 
is plotted in Fig. 4. On that figure there also appear curves 
giving the corresponding results for uniform heating all around 
the periphery. A comparison may be made at the same heat in- 
put per unit length (i.e., same Q). It is seen that insulating the 
end walls leads to a significant decrease in the hot spot tempera- 
ture at the duct corners. 

Velocity Distribution in a Circular Sector. The problems heretofore 
considered have been concerned with rectangular geometry. 
Now, we turn to a very different configuration: The case of heat 
transfer in a duct whose cross section is a circular sector. The 
co-ordinate system is shown in Fig. 2. 

The velocity distribution will be considered first. Using dimen- 
sionless variables 9 = 6/@ and R = r/ro, we write a first ap- 
proximation for the velocity distribution as 


u = (8? — 1)R2(1 — R*)(m + mk) 


+ (04 — 1)R*(1 — R*)(ma + mR) (40) 
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Fig. 5 Velocity distribution along symmetry line (@ = 0) for laminar 
flow in ao circular-sector duct, opening angle 33.4 deg 


where the n m, are constants to be determined by applica- 


tion of the variational method. Equation (40) is introduced into 
the variational integral (Sa), and the minimizing process is carried 
out by differentiating with respect to each of the m’s and setting 
the resulting equations equal to zero. In this way, four equations 


are provided for computing the m, . m,. The numerical cal- 


culations were carried out for a duct of 33.4-deg opening angle so 
that a direct comparison could be made with the exact solution of 
relerence ite 


The mean velocity may be computed by integrating as follows 


re Oe / To 6 
fi f ur dr d@| fi fi r dr d@ 
0 0 / 0 0 


Then, the normalized velocity distribution resulting from the 
variational proc edure is 


(41) 


4 l 


——- {(@ — 1)R? 
7 0.007575 


1 — R*)(—0.05724 — 0.004605 R) 


(4 1)R*(1 — R*)(0.01675 — 0.033873 R)} (40a) 


The velocity distribution along the symmetry line (@ = 0) has 


been calculated from equation (40a) and is plotted in Fig. 5. 
Also included is a curve representing the exact solution [8]. The 
agreement is seen to be very good 

The relationship between the friction factor and Reynolds 
number may also be evaluated once the velocity distribution is 


known. Based on equation (40a), we find that 


c,*Re = 53.79 


The exact solution gives a value of 53.75 for the constant. 

Circular Sector With Uniform Peripheral Temperature. Having deter- 
mined the velocity distribution, we now turn to the heat-transfer 
problem. We shall consider the situation where the walls have 
a uniform temperature around each cross section, while the heat 
transfer per 


Under these 


unit length in the axial direction is a constant. 
conditions, the variational integral of equation (9a 
applies without modification 

For convenience, let 7' in this equation represent the difference 
between the temperature at any point in the fluid and the wall 
temperature 


consequently, 7’ = 0 on the boundaries. It is 


expected that the reduced temperature distribution would have a 
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shape generally similar to that of the velocity. So, we select 
a function of the same form as equation (40). This expression, 
along with the velocity (40a), is introduced into the variational 
integral (9a). The end result of the computational process is 
T-T, 3 = 


6? — 1)R*(1 — R*)(0.037613 + 0.120612) 
Q/kOo 


+ (@* — 1)R%1 — R*)(—0.010928 — 0.010200 R 42 


The bulk fluid temperature Is obtained by integrating as follows 


my 
hal ” uB r dr dO 
0 0 
- a 
f f vi ur dr dé 
0 0 


Finally, the Nusselt number may be evaluated from the basic 
definition 


8, = —0.01629 


14) after observing that the equivalent diameter and 
heated perimeter for a circular sector are 2re/(O + 1) and 2ro(4o 
+ l 


, respectively. So, 


’, ( 9 y 1 
Nu oe -{ —— 
4+ 1 Oy 


This is about 3 per cent higher than the value given by the exact 
solution in reference [8]. 


Concluding Remarks 


The specific problems considered here are meant to illustrate 
the method and by no means exhaust the situations to which the 
variational technique may be applied. In addition to numerous 
other laminar flow configurations, it may be possible to make 
application to the turbulent situation, provided that suitable 
choices ar’ made for the trial functions. 

The authors are currently extending the variational method to 
the thermal entrance region of ducts. 
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DISCUSSION 
Andrew A. Fejer’ 


The variational method developed in this interesting paper 
represents a valuable addition to techniques being used currently 
in solving heat-transfer problems. The authors demonstrate that 
this method can be used effectively to calculate heat-transfer 
rates and temperature distributions in laminar flows occurring in 
ducts of various shapes and for various types of thermal boundary 
conditions. And, as indicated by the authors, the method will, no 
doubt, be found useful also in connection with other types of heat- 
transfer problems. However, there are some questions suggesting 
themselves regarding the results obtained in some of the illustra- 
tive problems. In the case of the rectangular duct with uniform 
peripheral wall temperature and uniform axial heat flux, i.e., for 
oT,,/Ox = OT ,,/dy = 0 and dT7,/dz = c, the authors find that 
the Nusselt number increases from 3.6 at aspect ratio 1 to 6.79 at 
aspect ratio 10, approaching 8.23 as the aspect ratio goes to ~. 
This result is in good agreement with the results of the numerical 
solution of Clark and Kays.® 

On the other hand, for uniform peripheral and axial heat flux, 
i.e., O7',/O0n = const and d7',/dr = const, the authors find a 
Nusselt number of about 3 over a range of aspect ratios from 1 to 
10. Itis difficult to reconcile this finding with the fact that for in- 
finite aspect ratio (i.e., the two-dimensional case) both temperature 
and heat flux are constant along the periphery of the duct and 
Nu = 8.23 in this case. In this connection it may be of interest 
to examine the form of the Nusselt number as defined by the 
authors for the case of uniform peripheral heat flow: 


where ¢, is the average wall temperature along the duct periph- 
ery. In view of the fact that Nusselt number may be considered 
as a ratio of the temperature gradient at the wall to the tempera- 
ture difference across the fluid in the duct, it would, perhaps, be 
more meaningful physically to define a local Nusselt number along 
the periphery and to use its average to describe heat-transfer 
processes when peripheral variations of temperature occur. 


William C. Reynolds’ 


Comments based on research work in progress on Foundations of 
Irreversible Thermodynamics under sponsorship of Air Force Office of 


Scientific Research Contract AF49(638)-201. Reproduction in whole 
or in part for any purpose of the U. S. Government is authorized. 


The Ritz-Galerkin method is indeed a powerful one, and it is 
gratifying to see it used effectively in convection problems. 

To the uninitiated, variational methods may seem rather like 
magic, particularly if the physical significance is not apparent. 
The variational formulations of elasticity phenomena can be de- 
veloped from notions of “minimum strain energy’’; one might 
expect that the variational equations for laminar flow might also 
be derived from some sort of a physical extremum principle, and 
this is in part true. There is a particular type of variational 
problem called the “isoperimetric”’ problem in which an extremum 
is sought for one integral, subject to the condition that another 
integral is a constant [1]. Classically one seeks the curve of fixed 
length which encloses the biggest area, and the isoperimetric solu- 
tion is the circle. The hydrodynamic portion of the laminar flow 
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problem [authors’ equation (8)] can be formulated in this manner, 
with an extremum being sought for the integral 


D= f, M(u,? + u,*?)dA = min 


subject to the restraint imposed by continuity that 


f, pudA = Au,p = w = const 


By formulating the problem in this manner a certain physical 
significance can be seen; the rate of dissipation of mechanical 
energy D [2] is to be a minimum, subject to the condition that 
the total flow rate is constant. 


(43) 


(44) 


That is to say, one examines all 
possible velocity profiles which have the same total flow rate, and 
zero velocity at the wall; the actual velocity profile is the one for 
which the rate of dissipation is the least! This physical formula- 
tion is quite appealing, and not too unreasonable; it can be looked 
upon as a “path-of-least-resistance”’ statement. Several well- 
known laminar flow problems have been handled by W. L. Rogers 
[3] using minimum dissipation as a basic postulate instead of the 
momentum principle. These include essentially the whole class of 
fully established laminar flows. 

This physical picture may provide the insight which will 
eventually let turbulent flow problems be treated in the same man- 
ner. It may be easily shown from friction factor measurements 
that the dissipation in a given tube for turbulent flow is considera- 
bly greater than the dissipation in laminar flow at the same flow 
rate. Malkus [4] has proposed a theory for fully developed tur- 
bulent flow between infinite parallel plates in which he extremizes 
the dissipation rate. He suggests that the flow is such that the 
dissipation rate isa maxrimum. It may well be that it is in reality 
again a minimum, this time selected from a different family of ad- 
missible flows. As an analogy, consider the classical “soap film’’ 
problem in which a soap film is stretched between two co-axial 
opposed wire hoops; elasticity theory says that the soap film is 
shaped such that the surface area isa minimum. The variational 
formulation yields a calenary as the surface of revolution for the 
film connecting the two hoops. However, as the hoops are moved 
farther and farther apart, the catenary film becomes unstable; 
the film eventually breaks, and the soap forms a thin circular mem- 
brane inside each hoop. The resulting discontinuous surface 
consisting of two disks is another surface of least area, but it is 
decidedly different in character. Perhaps the same sort of thing 
happens in transition from laminar to turbulent flow. Initially 
the laminar flow is the one of least dissipation; as the flow rate is 
raised, the laminar flow becomes unstable, and another flow of 
minimum dissipation appears. In a sense the laminar dissipation 
is the “minimum-minimum”’ and the turbulent dissipation is the 
“maximum-minimum.” It can be hoped that continued effort 
will lead to a variational formulation of turbulent flow phe- 
nomena which will complement the laminar formulation pre- 
sented by the authors. 

The heat-transfer part of the problem may also be formulated 
in an isoperimetric manner. If the flow is isothermal, the varia- 
tional formulation is equivalent to seeking the temperature profile 
which provides minimum production of entropy [5, 6] subject to 
the condition that the heat input is fixed. But this is not very 
satisfying, since heat transfer requires temperature variations. 
Perhaps continued effort in the area of irreversible thermo- 
dynamics will shed more light on the physical significance of the 
variational formulation for the temperature profile. 

Finally, it should be pointed out that the authors fail to verify 
that their integrals are minimums and not maximums. This is 
quite important, for the relative accuracy of solutions can be 
readily determined by comparison of the values of the cor- 
responding integrals, the best solution being the one for which the 
integrals are least (or biggest, as the case may be). The 
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condition that the integrals be a minimum is given by the second 


arvation (|7| and in the authors’ nomenclature is 


Fy,v, > 0 
Fyyvy > U 


By inspection of the integrals of the authors’ equations (8) and 


(9), it is seen that both extremals are indeed minimums 
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Myron Tribus’ 


rhe method used by the authors is of greater generality than 
they claim. Two problems which can be solved by their methods 
but let 


It is only 


ire as follows 


Consider the square duct shown in Fig. 1; 
there be a nonuniform temperature on the periphery. 
necessary to invent a function f(z, y) such that f(z, y) is analytic 
inside the duct and is equal to the prescribed wall temperature on 
the duct wall Then one considers that the t mperature is com- 
posed of two parts according to 7’ = 7' f substituted in Equa- 


tion (4 


to give 


(° l'; <2 é F ( o*/ o*f 
/ + = 4 
Ox? oy? Ox? oy* 


for which the boundary condition is 7); = 0 on the walls. Instead 


of Equation (9), one obtains 


e444 Se : r ‘ = ( j Q ) 
7 ia ' " \ a Ak 
wy. wie , 
2 . Tp drdy (45 
Ox? or? f 


which may, of course, be solved by the methods given in the 


paper. Therefore the methods given can be used in the case 
where the wall temperature is nonuniform 
The authors’ method may also be used for nonhomogeneous 

boundary conditions. Again consider the duct of Fig. 1 but let 
there be external thermal conductances such that around each 
boundary 

oT 

hT 
on 


For example, leth = h along zr = +a,h = hy 
h = h 
To handle this kind of boundary condition substitute 7 = 


along z = —a, 


along y = +a,h = hy along y = —a. 


Or, y Ty r,u)in Equation (7) to get 


gV27T; + 2V9°VT: + T1029 — Gz, vy) = 0 (46) 


* Associate Professor, Department of Engineering, University of 
California, Los Angeles, Calif. Mem. ASME. 
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where 


Dr. Sparrow has pointed out to me that, if we choose for our func 
tion ga function such that dg/@n = gh on the contour, the bound- 
ary condition on 7, becomes 07)/0n = 0. This condition on g 1s 
the same as 0 In g/On = h, hence it is usually not difficult to find 
the needed function g. It can be verified by the use of Equation 


1) that letting 


avg .. 
] ZT , 


T.,?) + GT; - 
leads to Equation (46). For the case considered a suitable func- 
tion for g 1s 

l 


Ing = . *hs (y — a)*he + (x 


ta 


48) 


Using Equations (47) and (48), the authors’ methods may there- 
fore be applied to nonhomogeneous boundary conditions. By 
appropriate choice of the function g, other more complicated 
handled. It is probable that 
methods for doing so are already well described in the literature of 
other fields. 
common and the authors have rendered an important service in 


boundary conditions may be 


The use of these methods in heat transfer is not 


bringing this powerful technique to the attention of the readers 
of this journal. Applications to steady-state conduction are of 
course obvious. The writer looks forward to learning how the 
methods of Biot can be applied to convection-problems and hopes 


that the authors will report on this work soon. 


Authors’ Closure 


The authors would like to thank the discussers for their very 
worth-while comments and additions to the content of the paper. 
We agree with Professor Fejer that physical reasoning is de- 
ceiving in regard to the limiting case of infinite aspect ratio for a 
Evidently, the heat addi- 
tion at the short sides of a flat rectangular duct influences the 


duet with uniform peripheral heat flux. 


fluid temperature distribution in such a way that the tempera- 
tures along the long sides of the duct become quite nonuniform 
Thus, in the limit, the case of infinite parallel plates is not ap- 
proached, since for this configuration the wall temperature would 
be uniform even for uniform wall heat flux. The Nusselt number 
defined in the paper is one of practical interest in that it provides 
a means for computing the average heat-transfer performance 
Actually for the uniform peripheral wall heat flux case, the details 
of the wall temperature are of greatest significance, and it is for 
this reason that Fig. 4 has been included. 

Professor Reynolds’ comments on the general physical in- 
terpretation of variational methods are very interesting, especially 
in relation to turbulent flows. We are looking forward to seeing 
more of this work in the literature. 

We would like to thank Professor Tribus for adding appreciably 
to the generality of the variational treatment by considering other 
types of wall boundary conditions. There are several qualifying 
remarks which appear appropriate to Professor Tribus’ condition 


oT 
on 


—AT =0 (49) 


First, it should be noted that A is a conductance; that is, a heat- 
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transfer coefficient divided by a thermal conductivity. Also, the 
negative sign in the equation implies that n is the inward normal 
rather than the outward normal used in the paper. Of even 
greater importance are the restrictions which are imposed on the 
external temperature by the requirement of uniform axial heat 
transfer. To explore this, we rewrite condition (49) displaying the 
external temperature 7’, and using k and / for thermal conduc- 
tivity and heat-transfer coefficient, and n for the outward normal, 
mt 

oT | os aes 

k + h ( i — 7 

on 


(49a) 


e 
w 


To prepare uniform axial heat transfer, the integral of the local 
heat flux around the duct periphery must be a constant independ- 


-f i We 


So, we have 


ent of axial position 


oT 
k ds = constant = 
on 


Journal of Heat Transfer 


Now, the wall temperature is of the form 
oT 
Oz 


z+m(a2,y) 


where 07'/0z is a constant for fully developed heat transfer. Thus 
the integral on the right side of equation (50) yields the condition 


oT ” 
h z+ T, | ds = constant 
Oz 


In order for this to be satisfied, we have the restriction that 7, 
must vary linearly with z with the same slope as 7. This may 
significantly limit the practical applicability of Professor Tribus’ 
development relating to the condition (49). 

In reference to his question about future work, we hope, as 
previously mentioned, to show how variational methods can be 
applied to the thermal entrance region of ducts 
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The Effect of Surface Roughness on the 
Convection Heat-Transfer Coefficient for 


R. T. LANCET 


Los Angeles, Calif. 
Assoc. Mem ASME 


Fully Developed Turbulent Flow in Ducts 
With Uniform Heat Flux 


Experimental data are presented for the heat-transfer coefficient and friction factor in a 


smooth: and a rough duct with a hydraulic diameter of approximately 0.035 in. 


The 


flow was fully developed and turbulent, and the heat addition was uniform over the 


length of the tube. 


The rough tube indicated appreciable increases in heat-transfer co- 


efficient and friction factor. The smooth-tube friction factors corresponded to rough-tube 
values, indicating the difficulty involved in obtaining smooth surfaces for very small 


ducts. 


I. SOME industrial applications it becomes necessary 


to make heating or cooling ducts so small that the roughness re- 
sulting from normal manufacturing techniques becomes large 
with respect to the tube radius. In such cases the heat-transfer 
coefficient and pressure drop may increase appreciably over that 
for a smooth pipe. In still other cases, it might be desirable to 
provide a rough surface to increase the heat-transfer rate where 
space is more important to power loss as a result of increased 
pressure drop. 

The heat-transfer data for rough pipes is meager and of neces- 
sity limited to specific types of roughness. To date the primary 
type of roughness studied was of a continuous-ring type or thread 
type [5and6].?_ These tests indicated an increase in heat-transfer 
rate but also gave an even greater increase in friction factor. 


The object of this experiment was to obtain some data on the 


' This work was performed for General Electric Company and as 
part of a Master of Science thesis for the University of Cincinnati. 

* Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division and presented at the 
Annual Meeting, New York, N. Y., November 30—December 5, 1958, 
of Tae American Society or MecHanicaL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, August 5, 
1958 Paper No. 58 A-122 


Nomenclature 


effect of roughness in a very small duct and particularly to evalu- 
ate the effect of finite protrusions as compared to the continuous- 
ring type of roughness. 


Description of Equipment 


A cutaway view of the test section is shown in Fig. 1. The 
plenum chamber is preceded by a dense screen to provide a rela- 
tively uniform inlet velocity. It is insulated on the inside by 
teflon (which has a thermal conductivity of about 1/60 of the 
surrounding steel) to keep the air-temperature rise from the 
thermocouples to the test piece to a minimum, Fig. 1. Further- 
more, the outside of the plenum chamber is not insulated to 
provide a heat sink for the heat conducted to the plenum chamber. 
Armstrong A-1 cement was used to insulate the air while passing 
out of the plenum chamber to the test piece, Fig. 1. Three 
chromel-alumel thermocouples were installed in the plenum 
chamber as well as two static-pressure taps. One thermocouple 
was tack-welded to the outside for check-calculation purposes. 

The smooth test piece consists of a 4.5-in-long tube with a 
0.025-in. K 0.125-in. duct with 0.050 in. thick walls, rounded 
at the outside corners with approximately a 0.05 -in. radius. The 
surface roughness in the direction of air flow varied from 3 to 6 
microin. (based on an arithmetical average) and 21 to 24 microin. 


in the direction normal to the air flow. The former corresponds 





fluid free-flow area 
tube-wall cross-sectional area 
tube-wall log-mean area 
specific heat at constant pressure 
effective hydraulic diameter 
Fanning friction factor 
gravitational constant 
heat-transfer coefficient based on 
local wall-surface temperature 
and mixed fluid 
temperature 
k fluid thermal conductivity 
k, = tube-wall thermal conductivity 
k, = mean protrusion height 
L = tube length 
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Mach number 

Nusselt number (based on local 
mixed fluid temperature and 
local wall temperature), Nu,, 
= h,,Dy/k 

heated perimeter, for rough sur- 
face use Maximum perimeter 

Prandtl number, Pr = uC p/k 

static pressure 

total pressure 

corrected radial heat flow 

radial heat flow per unit area 

maximum tube radius, for rough 
tube taken as 1/2 duct maxi- 
mum height 


Reynolds number, Re = paD,/pu 
tube-wall thickness 
fluid temperature 


wall-surface temperature at maxi- 
mum radius 

outside-tube-surface temperature 

mean axial velocity 

cooling air flow 

co-ordinate parallel to pipe axis 

Darcy friction factor 

absolute viscosity 

kinematic viscosity 

fluid density 
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Fig. 1 


to a radius to mean protrusion height ratio of 
(the radius 
tube 


ipproximately SOO 


considered smooth ised here was half 


The 
and brazed together along with a 0.030-in-thick disk on each end 


which will be 


the maximum duct height). was made in two pieces 
Gaskets were 
and the 


heavy back-up plate 


to fasten to the plenum and mixing chambers 
placed between the and the 


piece was screwed in place using a relativel) 


test piece chambers test 
The back-up plate had four concentric rings to hold the nichrome 
strips which provided radiation and convection insulation since 
the spacing between the rings was sufficient to adequately inhibit 
free convection. 

The tube 


which were 


was instrumented with twelve thermocouples, two of 
The 


wrapped with two insulated 0.013-in-diam constantan wires, 


lost during assembly, Fig. 1 tube was then 


paralleled for the desired resistance. This was then wrapped with 


asbestos chord to provide insulation and make the unit more 
rugged during assembly. The tube 
tric sheet (0.003 


final layer was instrumented with 10 ther 


was wrapped with 4 concen- 
thick). The 
mocouples and wrapped 
had 8 coils of 
terminal boards 


chord 


lavers of nichrome ind 0.005 in 


with 9 coils of nichrome wire (the second tube 


nichrome ribbon) each connected to over- 


head. This was again wrapped with asbestos and some 
cement applied. 

The rough tube was 4.5 in. long and the maximum width and 
and 0.045 in., 
approximately 0.010 in. XK 0.010 in. 
spaced approximately 0.010 in. 
of 0.025 in. X 0.125 in., Fig. 2 


the completion of the tests 


height were 0.145 in. The protrusions 


<x 0.010 in. 


respectively. 
were and were 
apart, leaving an open duct area 
Fig. 3 shows the protrusions at 
Based on the maximum radius, the radius to protrusion height 
ratio R was 2.5. 


primarily for ease of manufacture 


The relatively large protrusions were selected 


and to show whether the rela- 
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Fig. 2 Rough duct cross section taken normal to axis 
































tively high heat-transfer coefficients indicated as possible by a 
previous analysis (not discussed here) could actually be achieved. 
The wall thickness was again 0.050 in. (exclusive of the protru- 
rounded to about a 0.025-in. radius. 
The disks at the tube ends were increased in thickness to 0.060 in 
because 


sions) and the edges were 
of manufacturing and resulting sealing difficulty ex- 
perienced with the first tube. 

The tube was assembled in a fashion similar to the smooth tube 
with the exception of having 16 operating thermocouples, an im- 
proved guard heater including nichrome ribbon instead of wire, 
This was 
done in an attempt to improve the low Reynolds-number results 
(A comparison of Figs. 8 and 9 indicates that the end effects were 
reduced. ) 


and appreciable insulation around the guard heater. 


The main function of the mixing chamber was to provide an 
the 
the duct impinged on a flat mixing strip. 
thick L-605 sheet 


accurate value of exit-gas temperature. The air leaving 
A nozzle of 0.003-in- 


metal directed the air at three CA thermo- 
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Fig. 3 Protrusions at completion of tests 
couples The air then flowed over the thermocouple shank twice 
while following it path to minimize the heat losses of the tube exit 
The 


vas further insulated on the outside with about 1 in. of 


ur before reaching the thermocouples, Fig ] MUXng 
hamber 
insulation. For the second tube the insulation was improved and 
increased in thickness to about 2.5in. In addition, thermocouples 
were mstalled on the flange and outer wall of the mixing chamber 
The downstream static pressure was me sured and a static tap 
placed in the upstream wall of the mixing chamber as shown in 
Fig. | 
sottled air 
duced b 
The electrical system for the second tube differed from the first 


tests oaly in that 8 guard heater 


dew point, 60 F) was used. The pressure was re- 


two regulators to about 100 psia at the plenum chamber 


coils were used rather than 9 and 
that a constant-voltage transformer was included in the guard- 
heater system The 8 guard-heate r coils were controlled by 8 
separate variacs. The sensitivity of these secondary variacs was 


increased by feeding them by a Oth primary variac which was 


connected to the 110-volt a-ce supply The tube eiectrical svstem 
was similar to the guard-heater supply but had a single variac for 
control Ammeters and voltmeters were installed and the ex- 
perimental tube power was computed from the current and re- 


sistance of the constantan wire 


Reduction and Analysis of Data 


The actual dimensions of both the smooth 
Only the 


ised and they are given in Tables | and 2 


and rough tubes 
differed from the design value measured values were 

The evaluation of the local heat-transfer coefficient was made 
by considering a differential portion of the duct as shown in Fig. 4 


rhe details are given in [1 The resulting equation is 


Table 1 Smooth-tube dimensions 


Design Measured 
0.025 0.0229 
0.125 0.1297 
0.00312 0.00297 
0.300 0.305 
0.0417 0.0389 


Item 
Height, in. 
Width, in. 
Area, sq in. 
P»,in. 

Du, in. 


Table 2 Rough-tube dimensions 


Measured 
0.045 
0.025 
0.1474 
0.1274 
0.00318 
0.0418 
0.385 
0.305 


Item 
Maximum height, in. 
Minimum height, in. 
Maximum width, in. 
Minimum width, in. 
Minimum area, sq in. 
Minimum Dag, in. 
Maximum P,, in. 
Minimum P,, in. 


Design 
0.045 
0.025 
0.145 
0.125 
0.00312 
0.0417 
0.380 
0.300 


The second term on the right accounts for the net heat conducted 
along the tube rather than to the cooling air in the differential 
length dz. 

The heat generated per unit length was calculated from the 
measured resistance of the coil divided by the length of the wind- 
ing. The resistance change due to temperature was calculated and 
neglected when found to be negligible (the temperature coefficient 
for constantan at 25 C is 2 * 10-*ohm/ohm deg C 

The heat generated was then corrected by the over-all heat- 
balance deficiency. That is, the heat picked up by the cooling air 
plus that conducted out the ends of the tube was subtracted from 
the total heat generated and the difference was assumed to be lost 
The heat 
picked up by the cooling air was calculated from the inlet and exit- 


iniformly along the tube length to the surroundings 


ir temperature measurements, the measured weight flow, and 
the average specific heat. The corrections varied from 1 to 3 per 
cent for runs equal to and above Re = 8000, from 7 to8 per cent 
for Re 5000 and from 12 to 13 per cent for Re 3000. The 
corrections increased to 33 per cent for Re = 1000, and, as a result, 
The values for Re = 


ill values for Re < 3000 were discarded. 


3000 are included in Fig 5 as indicating the transition region 
and are not used further 

The heated perimeter was taken as the inside tube perimeter 
normal to the air-flow direction 


The 


equation 


wall temperature was calculated from the following 


Qt 
K w+ { ln 


The local cooling air temperature was calculated in a stepwise 
manner as shown in the following for a portion of the duct lying 
between two thermocouples a distance z apart as in Fig. 4 


% 
LLL. M44 
$s 
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WLLL 
© Th. 


ww I 














Fig. 4 Duct cross section taken parallel to axis 
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The resulting equation is 


ses 
vas _ T’c intet rT 


AT. 
where 


1 aT 
(QAzr +k, A, (are, dei —- —= ) 
) i 1 


we, dz 


The details of the analysis are in reference [1] 
term (d?7’, 
versus X. 


The remaining 
(dz?) was measured from a graph of (dT7,,)/(dz 

All the fluid properties were evaluated at the local film tem- 
perature using [2] which contained data obtained from [3]. 

The friction factor was calculated assuming the fluid was com- 
pressible, that there was an entrance pressure loss equal to 1/2 the 
difference between the total and static pressure at the inlet, and 
that the exit static pressure was equal to the mixing-chamber 
pressure: In addition, the flow area for the rough tube was taken 
as that between the tips of the protrusions 
Fig. 2). 


(shaded portion of 


Results and Discussion 


The dimensionless heat-transfer coefficients (Nu,,) for the 
smooth and rough tubes are graphically presented in Fig. 5 and 
tabulated in Tables 5 and 6. The increase in Nusselt number for 
the rough tube is almost twice that of the smooth tube. Included 
in Fig. 5 is the Colburn equation (4) for smooth round duets in 
turbulent flow, and it is seen to fit the smooth-duct experimental 
data very well in the turbulent-flow range. The dashed curve was 
fitted to the rough-tube data, assuming it should be parallel to 
the smooth curve and that the data for Re 5000 was in the 


transition region. The resulting equation is 


Nu,, = 0.042 Re®* Pr’? (4 


as compared to the Colburn equation for smooth round pipes 
Nu,, = 0.023 Re Pr’? (5 


where the fluid properties are evaluated at the fluid film tem- 
perature. 
W. Nunner [5] 


that the protrusion size, shape, and spacing had the same type of 


assumed (based on his experimental results 
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Fig. 5 Comparison of smooth and rough-duct experimental results. 
(Heat-transfer data.) 
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Fig. 6 Comparison of smooth and rough-duct experimental results. 
(Friction data.) 
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Fig. 7 Comparison of smooth and rough-duct experimental results. 
(jn factor and j,/X ratio.) 








effect on the friction factor as on the heat-transfer coefficient and 
concluded that he could express the Nusselt number as a function 
of Re, Pr, A, and A» alone (where Xo is the friction factor for a 


smooth pipe). On this basis he derived the following equation: 


\/8 Re Pr 


Nu,, = : (6) 
f *Pr-/4Pr A/Xo — 1) 


Equation (6) is plotted in Fig. 5 for the experimentally deter- 
mined friction factor \ and can be seen to compare favorably 
with the experimental results. This is particularly interesting 
since the equation also agreed well with the continuous ring type 
of roughness considered in [5] although various rings and spacings 
were tested. 

The friction data are presented in Fig. 6 and it becomes im- 
mediately apparent that the “smooth tube’’ was not as smooth 
as assumed. When these results were obtained the smooth tube 
was reassembled and additional data taken, but no significant 
The tube was then 
cut open and the surface finish rechecked and again the results 
agreed with the previous measurements. 

Apparently the fine surface scratches, which run in the direc- 
tion of the tube axis, were sufficient to cause the smooth tube 


deviation from the previous data was found. 


to act like a rough tube. Even though these scratches are barely 
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Fig. 8 Smooth duct and guard heater temperature distribution for Re = 
8200 
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Fig. 9 Rough duct and guard heater temperature distributions for Re = 
8400 





Table 3 Heat-transfer results—smooth tube 


X, in Nu, Re 
1.5 3010 17.4 

2.0 2930 17.7 

2.5 2870 16.8 $950 
2.75 2820 17 4905 
3.0 2790 16! 1820 
3.25 It 1770 


Num Re 
5160 
5040 


Weighted mean 
values 


2010 Ww 4980 


Table 4 


Nu, Re 

18.25 

16.87 

16.00 

t 16.14 

3.0 15.73 

Weighted mean 
values 


5100 
5150 
$960 
1910 
{S70 


16.85 5050 


visible to the naked eye and cover a relatively small portion of the 
duct surface, they are apparently the controlling factor. As 
noted previously, the measured surface finish normal to the tube 
axis was 21 to 24 microin., which corresponds to an average peak- 
to-valley height of about 74 to 84 microin. This corresponds to a 
radius-to-protrusion-height ratio of about 150. In addition, there 
were very small traces of braze material in one of the tube corners 


As pointed out to the author in [7] this roughness may explain 


why the results agree with the Colburn equation with a coefficient 


172 


Nu» Re N te Nu, Re 
26.7 8510 12600 21200 
27 8316 12400 20700 
27 8169 12200 20400 

7 8090 12200 20200 
7950 11900 19800 
7860 11800 19700 


27 0 
27.0 
i. $212 


12300 20500 


Heat-transfer results—rough tube 


te Re Nu, te 
S760 15430 27700 
8460 15210 27300 
8230 14900 26600 
S100 14800 26200 
7980 14690 25800 
8380 


15080 26900 


of 0.023 rather than with the lower values of 0.019 to 0.021 now 
accepted for air. 

Fig. 7 shows that the ratio of j-factor to the friction factor 
which is given as constant by the usual Reynolds’ analogy) is 
appreciably higher for the really smooth tube 
Blasius friction factor) than for the rough tubes of this experi- 


(as given by the 
ment In general, this agrees with the previous concept that 
generally the gain in heat transfer due to surface roughness is 


smaller than the loss due to increased friction. The closeness of 
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Table 5 Friction-factor results—smooth tube 


Po, P 2, 
in. Hg in. Hg 
70 41.0 
120.‘ 58.7 0.298 
150 69.¢ 0.298 
181.8 82.¢ 0.297 
93 .: 49 0.297 
70 41 0.283 
48 34.2 0.236 


M, 
Ent. 
0.283 


Table 6 Friction-factor 


W, Pw, P 2, M 
Ibs/hr in. Hg in. Hg Ent. 
173 70.3 38.0 0.213 
187 91 43.3 215 
323 110. 3.5 229 
).057 127.! 53.3 . 225 
7.560 160 x 223 
7 .225 
216 
0.193 


577 202 
11.59 251 


2.282 55.:! 


the smooth-tube results to the rough tube 
roughness of the smooth tube is inefficient 


because the roughness was unintentional 


indicates that the 
(This is primarily 
In other words, it 
appears that a tube surface roughness can be designed which 
would have a greater Nusselt number for the same friction factor 
as the smooth tube (presumably reducing the height of the pro- 
trusions of the rough tube would result in such a tube). 


Conclusions and Recommendations 


a) The equation derived in [5] appears to be sufficiently ac- 
curate for predicting the Nusselt number and is recommended for 
use at the present time. 

(b) Very small tubes (Dy of the order 0.040 in.) will act like 


rough tubes unless extreme care is taken to insure a smooth 
surface. 

c) Appreciable gains in the dimensionless heat-transfer co- 
efficient (Nu,,) can be made by increasing the tube-surface rough- 
ness but generally at the cost of an even greater increase in fric- 
tion factor 
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This very interesting paper illustrates the large effects on heat 
transfer that can result from surface ro 


ighness. 


UO! some con- 
cern to the author in his discussion is the flow which existed in 
the smooth tube and the effects of the relatively minute cireum- 
ferential roughness. A point in question in the case of the smooth 
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M, f r Re 
Exit Fanning Darcy 
0.445 0.0125 0.0494 
0.549 0.0126 0.0503 
0.575 0.0130 0.0519 
0.587 0.0133 0.0531 
0.510 0.0120 0.0481 
0.443 0.0123 0.0491 
0.309 0.0130 0.0520 


12300 
22100 
27500 
33100 
17100 
12300 

7131 


results—rough tube 


M. Es d Re 
Exit Fanning Darcy 
0.370 0.0224 0.0894 
0.428 0.0242 0.0967 
0.482 0.0214 0.0858 
0.500 0.0229 0.0916 
0.525 0.0236 0.0955 
0.550 0.0237 0.0947 
0.555 0.0264 0.1055 
0.300 0.0256 0.1022 
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tube is the state of the boundary-layer development. Even at 
Reynolds numbers as small as 3000, for an appreciable portion 
of the length, fully developed tube flow will not have been 
achieved. Perhaps this can explain a deviation from fully de- 
veloped flow equations for smooth tubes 

The 0.023 constant used in the Colburn correlation does not 
seem to be unusually high, although perhaps it is higher than the 
latest data. Because it is on the high side, it is interesting to 
compare the thickness of the laminar sublayer in turbulent flow 
with the height of the surface roughness in the smooth tube. 
If the surface roughness is greater than the thickness of the 
laminar sublayer, the results probably will be somewhat affected. 
Unfortunately, the computation of the sublayer thickness is 
clouded by the need to estimate the shear stress and the existence 
of entrance conditions. Results show that in some of the test 
runs the height of the roughness may approach the thickness 
of the laminar boundary layer, but it is believed that generally 
the surface-roughness height will be less than the normal thick- 
ness of the laminar sublayer. This is not necessarily inconsistent 
with the relatively high friction factors found using the smooth 
tube since short tubes sometimes give apparent friction factors 
which are higher than would normally be expected. 

Since the foregoing discussion is centered around the results in 
the smooth tube, it illustrates the need for more work in the area 
of marginal roughness; that is, the thresholds of roughness which 
will produce increased heat-transfer characteristics in a tube of 
this general size. 


Myron Tribus‘ 


Just before “Conclusions and Recommendations,”’ the author 
states: “In other words, it appears that a tube-surface roughness 
can be designed which would have a greater Nusselt number for 
the same friction factor as the smooth tube.” 

The writer has always believed that whereas one can obtain 
an increase in heat transfer through the addition of roughness, 
the increase in pressure drop always exceeds the increase in heat 
transfer and therefore, the smooth tube is always the best heat- 
transfer surface if one is comparing heat transfer and pressure 
drop only. The writer would therefore like a clarification of 
what is meant by the foregoing statement. 


’ 
Author's Closure 
While for low Prandtl number fluids in laminar flow the inlet 
4 Professor of Engineering, University of California, Los Angeles, 


Calif. Mem. ASME. 
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length required to obtain fully developed flow will be of the order 
of 50 to 150 tube diameters, shorter lengths are required in tur- 
For example, Deissler (NACA Report 1210) shows 
analytically (for Prandtl number = | and Reynolds numbers from 
10* to 2 & 
developed turbulent value in 5 to 7 diameters 


bulent flow 


10°) that heating rates reach 99 per cent of the fully 
Furthermore, 
Boelter, Young, and Iverson (NACA TN 1451) give experimental 
data for air flowing in a tube with sharp edge entrance conditions 
which indicate that fully developed heating rates will be obtained 
in 15 to 17 diameters for Reynolds numbers as high as 36,000. 
Since local heat-transfer coefficients were determined in the range 
of 35 to 75 diameters, it is unlikely that they were influenced by 
entrance effects 

The friction factors, however, represent average values for the 
entire tube and must therefore be influenced by the entrance 
effects. However, the extent of the entrance effect can be esti 
mated, for the smooth pipe, by considering the effect of the en- 
trance on the average heating rates. Using the experimental 
data of NACA TN 1451 it was found by the writer that after 17 
diameters the average heating rate was approximately 18 per cent 
greater than the local value. Extending these data to a tube with 
an L/D of 115 


effect would cause the average value to be only 2.6 per cent 


corresponding to this experiment) the entrance 


greater than the local value Therefore entrance effects alone 


cannot expl iin the high smooth tube friction factors which were 
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from 33 per cent to 122 per cent greater than truly smooth tube 
While the writer does not completely understand these high val- 
ues of friction factor, it is felt that they may result from phe- 
nomenon which cannot be directly scaled down to the small diam- 
eters encountered in this paper. 

As it was pointed out in the last paragraph of the paper before 
the ‘Conclusions and Recommendations” the gain in heat 
transfer due to surface roughness is smaller than the loss due to 
increased friction.’’ This can be seen more clearly in the revised 
Fig. 7 where both the experimental data of this paper and equa- 
tion (6) from reference [5] show the ratio of the j( = Nu/RePr 
factor to friction factor for the rough ducts to be less than the 
theoretically smooth duct. Thus the data of this paper are in con- 
currence with the literature and Professor Tribus. The state- 
ment in question (see comment by Myron Tribus) is based on the 
assumption that not all types of roughness produce the same heat- 
transfer rates for a given pressure drop. Referring again to the 
revised Fig. 7 we see that equation (6) agreed reasonably well 


with the rough tube data. However, equation (6) indicated ap- 
preciably higher values of heating rates could have been obtained 
for the so-called smooth tube of this experiment if the roughness 
were properly designed. In other words, if one must accept sur- 
face roughness, proper design will produce the maximum heat- 
transfer rate for the given pressure drop. However, this value 


will be less than the theoretically smooth value shown in Fig. 7. 
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PROCEEDINGS, 3rd U. S. NATIONAL CONGRESS 
OF APPLED MECHANICS 


Published 1058 20.00 to members and 

nonmembers.—-Within this 904-page vol 

ume are reports of important studies in the 
fields of dynamics, vibration, elastic waves, elas 
ticity, elastic structures, plasticity, viscoelastic 
flow, fracture, fluid flow, aerodynamics, and heat 
transfer There are also discussions of new con 
cepts and methods, and a large amount of new data 
which will aid in generating new ideas and new 
approaches toward solving other problems of the 
applied mechanics field The information was pre 
sented by 155 well-known engineers at the lune 
1958 congress 


MECHANICAL IMPEDANCE METHODS 
FOR MECHANICAL VIBRATIONS 


Published 1958 $5.50.—Shows how im 

pedance methods apply to lumped and con 

tinuous systems of simple and moderate 

complexity reviews measurement tech 
mques, demonstrates the power of digital com 
puters by comparing the calculated and measured 
characteristics of a highly symmetrical system of 
moderate complexity, gives -measured values of 
typical structures of large size and high complexity, 
discusses the importance of the impedance in in 
fluencing shock and vibration spectra measured in 
field service, indicates how to apply impedance 
methods to the calculation of vibration insolator 
effectiveness, treats impedance of some disordered 
systems, and illustrates how impedance methods 
may be used to find the response to random excita 
ton 


SHOCK AND VIBRATION INSTRUMENTATION 


Published 1956. $5.00.—This book provides 

a fund of useful information on recent ad 

vances and new concepts in the technology 
of shock and vibration Among the subjects dis 
cussed are random vibration, effect of non-rigid 
structures on vibration isolation, hysteresis and 
slip damping. vibration and shock testing. design 
approaches and analog methods. CONTENTS 
Evaluation of Mounts Isolating Nonrigid Machines 
from Nonrigid Foundations. Experimental Study 
of the Effects of Foundation Resilience on Vibration 
Isolating Effect of Material and Slip Damping on 
Kesonance Behavior Shock and Vibration Eno 
vironments Mechanical Design for Random Vi 
bration and Shock Influence of Electrical and 
Motional Impedance on the Control and Perform 
ance of Some Vibration Machines. Shock Test 
ing Machines and Procedures Damaging Poten 
tial of Shock and Vibration. High Speed Comput 
ing Methods for Shock and Vibration Problems 


METALS ENGINEERING: PROCESSES 


Published 1958. $13. 50.—Presents detailed 

data on the various processes by which 

metals are converted into finished products 
For each of the manufacturing methods, there is a 
compilation of the basic physical characteristics to 
be considered and the general advantages and limi 
tations usually enceuntered. Concise, parallel data 
concerning the suitability of various metals for 
each process, and the tolerances on size and surface 
finish obtainable are included 


ENGINEERING TABLES 


Published 1956 £12.00.—A collection of 

tables often wanted by engineers but not 

commonly found in handbooks covering 
Bar Stock and Shafting. Conversion Factors 
Formulas for Stress and Strain Properties of Sec 
tions and Cylinders Bearings Bearing Load 
Analysis. Spur Gears. Helical and Herringbone 
Gears Bevel Gears Worm Gears. Cylindrical 
Fits Standard Tapers Keys and Keyseating 
Bolts Counterbores Serew lhreads Slots 
Serrations and Splines. Nuts. Pins Snaprings 
Washers Wrench Openings Springs Aircraft 
and Mechanical Tubing Pressure Tubes Pipe 
Pipe Threads and Fittings Electric Motors 
Graphical Symbols. Welding. Gaskets Hydrau 
lic Standards and Symbols. O-Rings Packings 
Seals Bibliogtaphy 


THE AMERICAN SOCIETY OF 


THERMODYNAMIC AND TRANSPORT 
PROPERTIES OF GASES, LIQUIDS AND SOLIDS 


Published 1959 $12.50 Provides a wealth 

of information that can be profitably applied 

to the solution of many problems in various 
areas demanding a knowledge of thermal properties 
of materials. Specifically, the book surveys the 
present theoretical and experimental state of the 
science, reports a large amount of data on experi 
mental and theoretical techniques; reviews and 
evaluates the present state of knowledge in the 
specific areas covered, and indicates the gaps of 
knowledge existing in both transport and thermo 
dynamic properties, particularly at high temper 
atures and at high or even modetate pressures 
The information was originally presented at a Sym 
posium sponsored by the Standing Committee on 
Thermophysical Properties of the ASME Heat 
Transfer Division in February, 1959 


PROCEEDINGS OF THE CONFERENCE ON 
LUBRICATION AND WEAR 


Published 1958 23.00—Here is the com 

plete record of the 1957 Conference arranged 

by the Institution of Mechanical Engineers 
with the collaboration of the ASME for the purpose 
of evaluating current lubrication knowledge, 
bringing out the findings of important experiments 
conducted in the United States, Canada, United 
Kingdom, and Europe, and obtaining the views of 
those participating in the discussion. Subjects 
dealt with include hydrodynamic lubrication 
boundary friction, boundary lubrication, bearing 
metals, novel bearing materials, glands and seals, 
solid lubricants, surface treatments, ball and rollex 
bearings, gear lubrication, engine lubrication, mis 
cellaneous lubricants and applications, additives, 
and wear The 34 papers on the latter subject cover 
laws of wear rate; mechanism of wear; testing of 
materials for resistance to wear; surface films; 
effective temperature and environment; wear due 
to abrasion, fretting, scuffing, pitting, wear on 
cutting tools, wear in engines, gears, etc.; influence 
and nature of rubbing surfaces 


THERMODYNAMIC PROPERTIES OF COMPRESSED 
WATER 


Published 1957. $2.50.—Here is a chart 

that fills the need for an accurate and con 

venient presentation of the thermodynamic 
properties of compressed water These properties 
which include temperature, pressure, specific vol 
ume, enthalpy and entropy, are presented graphi 
cally in a large scale chart broken up into fifteen 
plates The temperature range covered is from 32 
to 705 .4°F and the pressure range, from 0 08854 
to 6000 PSIA 


METALS ENGINEERING: DESIGN 


Published 1953. $10.00.—This book dis 

cusses the essential properties which need 

to be evaluated by the design engineer in 
his selection of one material over another. Com 
prising 48 sections and written by 43 well-kuown 
authorities, it deals with the over-all problem of 
selection of material and takes up such specific 
items as high temperature considerations, plasticity, 
residual stresses, vibration, fatigue, shot peening 
cold working, nitriding, flame strengthening, im 
pact, corrosion, non-destructive testing, surface 
finish and mass production, and design theory and 
practice 


METALS PROPERTIES 


Published 1954. $11.00.—Provides, in con 
venient charts and tables, data on a broad 
range of metals in common industrial use 
AISI steels, ASTM steels, cast copper alloys, 
aluminum alloys, tin, magnesium, etc Tabulated 
under each of the more than 500 metals listed is 
such information as the chemical composition of 
the metal; its brittleness, heat treatment and other 
characteristics; its industrial uses; treatment tem 
veratures for forging, annealing, quenching, etc.; 
uch technological properties as recrystallization 
cemperature and hot working temperature, and a 
great deal of other pertinent information to help 
the designer choose the proper metal for each pay 
or product 
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MECHANICAL ENGINEERS 


FLUID METERS, THEIR THEORY 
AND APPLICATION 


Published 1959. $8.00.—This Fifth Edi- 

tion includes information on a number of 

new types of fluid meters and taetering pro- 
cedures that have been developed in recent years 
It is divided into three sections with the first giving 
the classification and nomenclature of fluid meters, 
together with definitions of special terms and other 
general information. The theory of fluid measure- 
ment and the steps taken to develop practical work- 
ing equations from the theoretical relations are set 
forth in the second section Figures and tables 
for use in solving practical fluid measurement 
pro'!ems, along with examples illustrating their 
proper use, are contained in the third section 


FLOW MEASUREMENTS (Chapter 4, Part of 5 instru- 
ments and Apparatus) 


Published 1959. $5.00.—These rules apply 

when the primary element is an orifice, a flow 

nozzle, or a Venturitube. The fluid may be 
compressible or incompressible. The primary ele- 
ment may be installed within a continuous section 
of pipe flowing full, or at the inlet or exit of a plenum 
chamber Specifically this Chapter includes in- 
formation on the construction of these three primary 
elements, the recommended techniques governing 
tests, the necessary equations tor computing rate ot 
flow, examples to illustrate the application of typical 
data, a discussion of tolerances applicable to certain 
of the factors involved in the measurements, and an 
outline of the major advantages and disadvantages 
of various types of primary elements 


77-YEAR INDEX TO ASME TECHNICAL PAPERS 


Published 1957. $20.00 —This ig an index 

to the technical papers published in the 

ASME Transactions, the Journal of Ap 
plied Mechanics, and Mechanical Engineering 
for the past seventy-seven years It is divided 
into two parts with the first covering the papers 
published between 1880 and 1939, and the second 
indexing those published from 1940-1956. Each 
part has a subject and an author index. The 
former lists the titles of the papers chronologically 
under as many headings as are necessary and in 
cludes numerous cross-references to direct attention 
to related material. The alphabetically arranged 
authors’ index gives the headings and the year 
under which each paper is listed in the subjegt index 


VISCOSITY OF LUBRICANTS UNDER PRESSURE 


Published 1954. $3.00.—This publication 

reviews and coordinates twelve experimentai 

investigations, made over a period of thirty 
five years, on 148 lubricants comprising of 25 fatty 
oils, 94 petroleum oils, 17 compounded oils and 12 
other lubricants Data are coordinated by means 
of sixty tables in which the results originally ap 
pearing in diversified units are compared. The 
methods proposed for correlating viscosity-pressure 
characteristics of oils with properties determined 
at atmospheric pressures are reviewed and illus 
trated Experimental work on heavily loaded 
bearings, lubrication calculations, and additional 
techniques for viscosity are covered. Conclusions 
and recommendations are presented. Computa , 
tion of the temperature coefficients of viscesity and 
the method of computing pressure coefficients are 
also given 


DIAPHRAGM CHARACTERISTICS, DESIGN AND 
TERMINOLOGY 


Published 1958. $3.75.—This Manual de- 

fines a diaphragm and its performance char- 

acteristics, describes methods of measuring 
and representing them, shows how they are related 
and used, describes and illustrates the effects of de 
sign details and manufacturing methods on the 
characteristics. In the appendixes, the terms and 
notation used are defined, the equipment and 
methods of testing described, applications of the 
diaphragms classined and equations given 
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